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Создан в целях информирования читательской аудитории о новейших достижениях и перспективах в области 
механики, машиностроения, информатики и вычислительной техники. Издание является форумом для 
сотрудничества российских и иностранных ученых, способствует сближению российского и мирового научно- 
информационного пространства. 


Журнал включен в перечень рецензируемых научных изданий, в котором должны быть опубликованы 
основные научные результаты диссертаций на соискание ученой степени кандидата наук, на соискание 
ученой степени доктора наук (Перечень ВАК) по следующим научным специальностям: 


1.1.7 — Теоретическая механика, динамика машин (технические науки) 

1.1.8 — Механика деформируемого твердого тела (технические, физико-математические науки) 

1.1.9 — Механика жидкости, газа и плазмы (технические науки) 

1.2.2 — Математическое моделирование, численные методы и комплексы программ (технические науки) 

2.3.1 — Системный анализ, управление и обработка информации, статистика (технические науки) 

2.3.3 — Автоматизация и управление технологическими процессами и производствами (технические науки) 

2.3.5 — Математическое и программное обеспечение вычислительных систем, комплексов и компьютерных сетей (технические науки) 
2.3.7 — Компьютерное моделирование и автоматизация проектирования (технические, физико-математические науки) 
2.3.8 —- Информатика и информационные процессы (технические науки) 

2.5.2 — Машиноведение (технические науки) 

2.5.3 — Трение и износ в машинах (технические науки) 

2.5.5 — Технология и оборудование механической и физико-технической обработки (технические науки) 

2.5.6 — Технология машиностроения (технические науки) 

2.5.8 — Сварка, родственные процессы и технологии (технические науки) 

2.5.9 — Методы и приборы контроля и диагностики материалов, изделий, веществ и природной среды (технические науки) 
2.5.10 — Гидравлические машины, вакуумная, компрессорная техника, гидро- и пневмосистемы (технические науки) 
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Наименование Выписка из реестра зарегистрированных средств массовой информации ЭЛ № ФС 77 - 78854 
органа, от 07 августа 2020 г., выдано Федеральной службой по надзору в сфере связи, 
зарегистрировавшего информационных технологий и массовых коммуникаций 
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Аннотация 

Введение. Проектирование систем управления движением мобильных роботов требует построения 
математических моделей. Исследователи неоднократно обращались к этой теме. В частности, опубликованы 
работы, посвященные расчетам мультифизических процессов, моделированию движения различных видов колес 
в тех или иных условиях. Кроме того, рассматривались динамика деформируемых контактирующих тел при 
скольжении, качении и вращении, вопросы автономности и управляемости мобильных роботов. Отметим, что 
динамика и точность позиционирования колесных роботов во многом определяется трением. В литературе не 
представлены исследования динамики робота с дифференциальным приводом с учетом взаимосвязей эффектов 
трения скольжения, верчения и качения на основе теории поликомпонентного трения. Изыскания в этой области 
способны выявить новые динамические эффекты. Основываясь на полученных таким образом данных, можно 
заняться повышением точности позиционирования при построении математических моделей управления. Цель 
представленной работы — исследовать движение автоматического устройства с дифференциальным приводом с 
учетом трех моделей контакта: неголономного, кулоновского трения и поликомпонентного. 

Материалы и методы. В качестве базовой приняли схему двухколесного робота с дифференциальным 
приводом и безотрывным движением по опорной поверхности. Перемещение устройства обеспечивалось 
программным управлением. Динамику описали в виде уравнений Аппеля. Для вычислений использовали 
математические модели, по-разному учитывающие трение. Координация действий механизма формировалась на 
динамическом уровне. Управляющие воздействия — моменты двигателей колес. При визуализации исследуемых 
моделей задействовали встроенные численные методы системы \оЙат Мафетайса («Вольфрам математика») 
с минимальной точностью 10-5. 

Результаты исследования. При построении математической модели определили равенства для угловых 
скоростей колес. Учли наличие площадки контакта и вывели уравнения динамики робота с дифференциальным 
приводом. Элементами системы были проекции силы и момент, показатели верчения платформы, массы, угловые 
ускорения и инерция колес. Показали, как формируются управляющие воздействия в рамках неголономной 
механики. Описали модель двигателей, которые создают момент управления на ведущих колесах. Вывели 
решение как взаимосвязь индуктивности проводников электродвигателей и работы электропитания. Детально 
рассмотрели три модели, описывающие динамику робота с дифференциальным приводом. Первая — 
неголономная. Вторая и третья включали систему уравнений динамики робота с дифференциальным приводом 
для общего случая, при наличии площадки контакта. При этом во второй модели игнорируется время 
переключений в двигателе и задействуется трение Кулона. В третьей ввели параметр, определяющий скорость 
переходных процессов в двигателе, и задействовали разложения Паде. Это модель с поликомпонентным трением. 
Итоги расчетов показали в виде графиков. На них исследуемые модели визуализировали в виде кривых разного 
цвета. Сопоставление графиков показало, в каких случаях после завершения переходных процессов управление 
обеспечит требуемую точность. 

Это модели 1 и 2. В модели 3 программное управление генерирует ошибку в угловой скорости вращения 
платформы. Данную ошибку нельзя спрогнозировать в рамках 1-Й и 2-й модели. Во всех рассмотренных системах 
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скорость скольжения колес в поперечном направлении падает до нуля. Получено и обосновано условие 
безотрывного движения опорного колеса. 

Обсуждение и заключение. Программное управление допустимо в моделях, не учитывающих трение колес при 
одновременном скольжении, верчении и качении (общий случай пространственного движения). Однако важно 
учитывать взаимосвязь данных процессов и поликомпонентного трения. Это необходимо для более точного 
выполнения роботом программных движений. Установлено, что программное управление в модели, 
учитывающей трение верчения и качения, приводит к отклонениям от программных значений угловой скорости 
платформы. Полученные результаты можно использовать при построении системы управления 
с прогнозирующими моделями. 


Ключевые слова: программное управление роботом, управление в неголономной механике, модель трения 
Кулона, поликомпонентное трение, ошибка угловой скорости вращения 
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Введение. В современном мире растет востребованность автоматических мобильных колесных аппаратов, 
активно исследуются возможности координации их работы. Значительный научный и прикладной интерес имеют 
вопросы автономности и управляемости мобильных роботов. В частности, рассматриваются подходы к расчетам 
мультифизических процессов, моделированию движения конкретных видов колес в разных условиях. Некоторые 
изыскания фокусируются на динамике деформируемых контактирующих тел при скольжении, качении и вращении. 
Результаты этих исследований находят применение в робототехнике [1-3]. Управление мобильными колесными 
аппаратами рассматривается на примере электрических самокатов «Сигвей» (Зез\мау) [4]. 

Мобильные механизмы с простыми и эффективными системами управления движением колес нужны 
в первую очередь для работы в узких пространствах [5], для проведения исследований в сложных 
промышленных условиях [6]. 

В [711] показаны возможности точного моделирования динамики колесных аппаратов с учетом скольжения, 
верчения и качения колес. 

При этом в литературе не описаны явления, возникающие при комбинации эффектов трения. Цель 
представленной работы — исследовать влияние поликомпонентного трения на динамику робота 
с дифференциальным приводом. Задачи: построить математическую модель движения робота 


с поликомпонентным трением и смоделировать программные движения. 
Материалы и методы. Рассматривается схема двухколесного робота с дифференциальным 
приводом (рис. 1). Предполагается безотрывное движение колес радиуса К по опорной поверхности. 
Хх] 
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Рис. 1. Кинематическая схема робота: а — платформа; 6 — одно из ведущих колес 
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Центр тяжести платформы С смещен относительно геометрического центра платформы А по оси Х на 
расстояние Й, ‚ по оси У — на расстояние й,. Точки С; и С, — центры масс колес (индекс 1 соответствует 


левому колесу, индекс 2 — правому). 

Вводятся следующие обозначения: т —щ масса платформы робота; с, — момент инерции платформы 
относительно центра масс С; Ох, у — неподвижная система координат; х,, у, определяют положение точки 
Ав О,х,у,; у — угол курса платформы; ф, — углы поворота колес (индекс 1 соответствует левому колесу, 
индекс 2 — правому); АС =АС, =Ё — расстояние между колесами; М,,М, — моменты двигателей, 
приложенные к колесам. 

Предполагается, что на опорное колесо 4 (рис. 1) действует только нормальная реакция опоры. Во время 
движения опорное колесо располагается на оси Х на расстоянии А, от центра платформы (А, < Г). 

Результаты исследования 

Математическая модель динамики робота. При описании движения робота с дифференциальным приводом 


в рамках неголономной механики применяются выражения, которые учитывают ограничения, связанные 
с отсутствием скольжения колес по отношению к поверхности: 


0х, =00х, =0,ох,, =0, (1) 


Ок 


1х 
где Эк Ок,х> Оку ОК, —Щ_ проекции на оси Хи? скоростеи точек К, иК, › являющихся центрами площадок 
контакта колес с опорной поверхностью. 

Отметим, ЧТО ДЛЯ неголономной модели предполагается точечный контакт колес с поверхностью, 


т.е. площадка контакта каждого колеса вырождается в единственную точку контакта. 
Рассмотрим выражения: 


ок. =0-©ОЁ- Ав, 
ок =о+0- К®,, (2) 
Ок б Оку =": 


С их помощью можно получить равенства для угловых скоростей колес: 
о- ОР 5+ ор 
в, = ‚в, = : (3) 
К К 
В общем случае, когда предполагается наличие площадки контакта, получим уравнения динамики робота с 
дифференциальным приводом: 


т, [6-0-9(и +1,0) |+2т, (5-Ои) = А, соза, + А, соза, — Е эт а, — Е, та, 
т, [#+л.0+9(5-л,9) | +2т, (й+05) = Е, эта, + Е, та, — Е. 0$ а — Р,.с08 а,, 


(с. +2. +т, (+ )+2т,Р )+т, [(и+95)л, —й,(6-©и)] = 


4 
Мс +М.. +(5, со5а, — Е зш а, —Ё, с0за +Ё 5 о, ) Е, _. 
Ль®, =М,-М„ -(Е, соза, - Е, зто, К, 
Люб, =М,-М„, (=, сова, -Ё\, зто, ) К. 
Здесь К,Ё, — проекции силы трения на направления скоростей скольжения центров пятен контакта колес; 


Е, 5, —_ Проекции силы трения на направления, перпендикулярные скоростям центров пятен контакта колес; 


М.,М.,— слагаемые момента трения при верчении для двух колес; @ — угловое ускорение верчения 
платформы; &,,6, — угловые ускорения колес, проецируемые на их ось вращения; т, — масса каждого колеса; 
Ук, и Лк. — моменты инерции колеса относительно соответствующих осей. 


Для формирования управляющих воздействий будем использовать модель динамики робота, описанного в 


рамках неголономной механики [12]: 
ь 2Л М+М 
т [5-9 ол,о |+2т =, 
| | К К 
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т, Г(+^.5 +9(5-л,9))л, —^, (-л,0-0^,0) | у 


5) 
р. (М,-М,)г ( 
2 2 2 — 2 1 
(ле. +2. +т (в +®)+2т,Р) +25 кто —— 
Ее также можно представить в виде: 
м, +М м,-м,)Е 
Е = 1 те 2 „Мо — ( 2 ы , (6) 


где Ё, и М. — левые части уравнения (5). 
Модель двигателей на ведущих колесах. Опишем модель двигателей, способных генерировать момент 
управления на ведущих колесах М,,М,. 


Уравнения движения колесных двигателей постоянного тока можно выразить формулами [13]: 


4(1.) 


в в. @) 

М =сж, (^ =) 
Здесь Г, —_ индуктивность проводников электродвигателей; в — сопротивление электрическому току в цепи 
двигателя; и, —_ электропитание для двигателя к 1, —щ ток, проходящий через цепь якоря двигателя; СС —= 


конструкционные постоянные двигателей. 
Предполагается, что двигатели ведущих колес имеют одинаковую характеристику. Из (7) видно, ЧТО 


ре. (8) 


Учитывая (8), в первом уравнении из системы (7) получим дифференциальное уравнение относительно 
момента М, : 
Е... В 
е ыы 
оп о. —с0,. (9) 
2 


ё 


К 
Уравнение (9) разделим на множитель — и получим: 


с2 
Р.С. с сс 
- -М,+М, =-—® 0, - 2 ,. (10) 
[6 К К 
Введем обозначения: 
Ес сс [6 
ПЕК =, К, =. (11) 
с, К у У: 
В таком случае управляющие моменты можно описать: 

ТМ, +М, =К.Ц, -К\о,, (12) 
где К,К, — константы; 1, — параметр, определяющий скорость переходных процессов в электрической части 
двигателей. 

При игнорировании времени переключений в электромеханической системе двигателя: 


М, =К,И, -Ко,. (13) 


Авторы [14] получили похожие формулы. 

Исследуемые модели. Рассмотрим модели, описывающие динамику робота с дифференциальным приводом. 

Модель 1 описывается уравнениями (6) и (12), когда длительность переходных процессов полагается 
исчезающе малой 1 =0: 


: 2/7, Ки, -К КИ, -К 
т [6-л,о-Ол,О ]+2т,5+ Ку о = КОД 10, + 50. > 
1 р; р: к В? В 


2 


т, [(+^.0 +9(5-л,о))л,—й, (5-1,0 -9^,9). +2. 0+ 


14 
+0 (о. +2.Л. +т, (2 +й, ) + 2т, 1?) = (к, ВВ р : _ 
К 
Здесь ®,,®, вычисляются по формулам (3). 
Для описания модели 2 применяются уравнения (4), (13), а также модель трения Кулона: 
Е =-ЛМ,Е, =0,Мс =0. (15) 
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Модель 3 представлена уравнениями (4), (12). Кроме того, учитывается поликомпонентное трение, 
представленное в [12] с помощью разложений Паде: 


21 -К. 1, со [ ма с) 
оо, ЕЛЕЫ - со-[ “) 


Е =-ЛМ№, ° = › (16) 
21 -К. 1, сиз[ А-а) 
А Аа 
0 
Кпо 
Е, =- ЛМ, з (17) 
21 —К, 1, со п-ов 
п 
КГ, т. сон с) ок, Нок, +] 
п 
ешь [ А-а) 
' 27 
2 2, со-[ с) ох, +0 
рт Га 
Ме =-2 1, . (18) 
в - 
А, 08| ——а | |2. +9х, +] 
Ь 
М, = ль". (19) 
| || 


Здесь [,1[,,/, — постоянные величины, которые находятся по формулам в [12]; и — угол, характеризующий 
направление линейных скоростей скольжения относительно точек контакта колес КК; К — 


7 


коэффициент, знак которого зависит от направления качения. 
Для определения программного управления следует указать законы движения робота 5“ (1),О“ (1). По (5) 


можно выразить величины Ё, = у (0, Е = т (г), Мо=М ь (г) в зависимости от времени. Используем 
известные значения Р., М. для вычисления управляющих напряжений (/ и (0, , подаваемых на приводы колес: 


В М К М. 
Ч, ре Пка ко Це = пк -ь Ко ь (20) 
- К (2 7Й К, \ 2 р 


м м 
Движение системы с заданными начальными условиями 0(0)=0,5—,и(0)=0,2—,©(0) = 0,2 
с с с 


смоделировали в математическом пакете «Вольфрам математика» (\МоШат Мафетайса). 


т 1 
Рассмотрим 5“(0=0,5 —, 0“ (1) =0,3—. Для интеграции системы автоматически выбирали шаг 
[6 [6 


интегрирования и задействовали встроенные численные методы У\УоШат Мафетайса с точностью результатов 
не менее 105. 


У(В, м/с 


Рис. 2. Зависимость продольной скорости у. Зеленая кривая — неголономная модель 1, 
красная кривая — модель 2 с трением Кулона, синяя кривая — модель 3 с поликомпонентным трением. 
Синяя и красная кривые совпадают 
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Рис. 3. Зависимость поперечной скорости и. Зеленая кривая — неголономная модель 1, 
красная кривая — модель 2 с трением Кулона, синяя кривая — модель 3 с поликомпонентным трением. 
Зеленая, синяя и красная кривые совпадают 


©(0, 1 
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Рис. 4. Зависимость угловой скорости робота. Зеленая кривая — неголономная модель 1, 


красная кривая — модель 2с трением Кулона, синяя кривая — модель 3 с поликомпонентным трением 


Из рис. 2 и 4 видно, что управление гарантирует точную реализацию программных законов движения после 
завершения переходных процессов в моделях | (неголономная) и2 (с трением Кулона). В модели 3 
(поликомпонентное трение скольжения, качения и верчения) после завершения переходных процессов из-за 
программного управления возникает постоянная ошибка в угловой скорости вращения платформы. Во всех 
рассмотренных моделях скорость скольжения колес в поперечном направлении быстро снижается до нуля. 
Случаи совпадения зависимостей, полученных для разных моделей, показаны на рис. 2 и3. 

Для проверки выполнения условий безотрывного движения опорного колеса запишем теорему изменения 
кинетического момента системы в проекциях на ось 7: 


Лк, (в, +6.) =М,+М,-М, -М, -М, В, +т, эй, -(Е 0$ 0, — 2 зто +, соза, -Ё, уто, ) К. 


ьЬ 
Согласно (4): 
М+М, =-М, А, +т,8й,. 
Из последнего уравнения получаем условие безотрывного движения опорного колеса: 
тай, —-(М,-М,) > 0. 
С учетом (5): 


а. , 
т)0+ Вт [и,Ф+ 1,0] >0. 


Бйр://уезиК-допза.ги 


тей, — К(т, +2т, + 
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Ускорения робота не превышали 1 м/с?, поэтому данное неравенство выполняется в каждый момент времени. 

Обсуждение и заключение. Итоги моделирования позволили сделать ряд выводов. Программное управление, 
основанное на неголономной модели 1, приемлемо для роботов, конструкция которых не учитывает контактное 
трение колес, возникающее при комбинации скольжения, верчения и качения. Однако при формировании 
управляющих воздействий для достижения более точного следования программным законам движения 
необходимо учитывать модель поликомпонентного трения. Полученное решение можно использовать при 
построении системы управления с прогнозирующей моделью (то4е! рге@сйуе сопёго|. Такая разработка будет 
предметом дальнейшего исследования. 
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Анализ точности позиционирования при конструкционных 
погрешностях в установке меканум-колёс мобильной платформы 
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Аннотация 

Введение. Мобильные роботы, способные осуществлять всенаправленное движение, широко применяются 
в различных областях человеческой деятельности. Для обеспечения высокой точности позиционирования 
всенаправленных платформ с меканум-колёсами требуется разработать их детальные математические модели, 
используемые при построении системы управления движением. Из-за непростой конструкции меканум-колбс 
при построении всенаправленных платформ могут возникать различные погрешности, включая погрешность 
установки таких колёс на платформу. Влияние ее на точность движения платформы ранее не исследовалось. 
Целью данной работы является оценка ошибок позиционирования, возникающих из-за наличия 
конструкционных погрешностей в установке меканум-колёс, и анализ влияния указанных погрешностей на 
точность отработки программного движения при использовании управления на кинематическом уровне. 
Материалы и методы. Анализ точности позиционирования основан на математическом моделировании 
кинематики платформы с учетом конструкционных погрешностей в установке меканум-колёс. Для описания 
взаимосвязи между угловыми скоростями вращения колес и скоростями платформы используются условия 
непроскальзывания точек контакта по опорной поверхности. Численные расчеты проведены в математическом 
пакете Ио/ат Мафетайса. 

Результаты исследования. Получена формула для оценки погрешностей псевдоскоростей платформы при 
программном управлении, формируемом на кинематическом уровне. Проведена оценка погрешностей 
скоростей платформы для простых движений. По результатам расчетов показано, что погрешности скоростей 
являются значительными для роботов с меканум-колёсами, функционирующих автономно. 

Обсуждение и заключение. Результаты расчетов продемонстрировали существенное влияние погрешностей установки 
колес на точность позиционирования меканум-платформы и подтвердили необходимость учитывать указанные 
конструкционные погрешности при создании автономных меканум-платформ. Построенная модель кинематики робота 
позволяет прогнозировать ошибки отработки скоростей платформы, возникающие при программном управлении, 
атакже отклонения координат геометрического центра платформы от программного движения. Кинематическая 
модель может быть использована для повышения точности позиционирования за счет формирования управления 
движением платформы, компенсирующего влияние неточностей установки меканум-колёс. 


Ключевые слова: меканум-колесо, всенаправленная платформа, кинематическая модель, конструкционные 
погрешности, программное управление 
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Введение. Мобильные роботы (или всенаправленные платформы), способные осуществлять движение 
в любом направлении, широко применяются в военной, медицинской, образовательной и промышленной 
областях человеческой деятельности [1, 2]. Возможность всенаправленного движения может быть достигнута 
за счет оснащения мобильной платформы омни-колёсами различного вида (меканум-колёса, классические 
омни-колёса). Для обеспечения высокой точности позиционирования всенаправленных платформ требуется 
разработать их детальные математические модели, используемые при построении системы управления 
движением. Необходимость обеспечения высокой точности движения связана с потребностью обеспечить 
автономное функционирование всенаправленных платформ. 

Зачастую при математическом моделировании кинематики меканум-колёс принимаются допущения, 
упрощающие описание движения таких колёс, об этом, например, идет речь в [3—5]. Так, в работах [3, 4] 
рассматривается динамика меканум-платформы с учетом сил сухого трения качения и условий 
непроскальзывания точек контакта, выражающихся неголономными связями. А в статье [5] отмечается, что 
интерес представляет исследование влияния конечных линейных размеров роликов и вибраций, связанных 
с этим эффектом, на динамику меканум-платформы. 

Актуальными направлениями исследований всенаправленных колёс являются моделирование влияния 
конструкции роликов меканум-колёс и оптимизация формы роликов [6-9]. Необходимость оптимизации формы 
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роликов связана с тем, что при некоторых формах роликов омни-колёс могут возникать вертикальные вибрации 
точек контакта при движении всенаправленных платформ (подробнее см. в обзоре патентов [10]). 

Другим фактором, существенно влияющим на точность позиционирования меканум-платформы, является 
проскальзывание точек контакта меканум-колёс с опорной поверхностью. Одним из способов учета 
проскальзывания точек контакта при описании кинематики является использование допущения, что скорость 
скольжения пропорциональна произведению угловой скорости вращения колеса и радиуса меканум-колеса, 
причем коэффициент пропорциональности считается постоянным [11, 12]. Это допущение позволяет свести 
уравнения кинематики к задаче с неголономными связями. Однако результаты экспериментов показывают, что 
указанный коэффициент изменяется при движении платформы [13, 14]. 

Стоит отметить, что из-за непростой конструкции меканум-колёс при построении меканум-платформ могут 
возникать различные погрешности изготовления, в том числе погрешность установки меканум-колёс на платформу, 
малые отклонения формы роликов и др. Влияние погрешностей установки меканум-колес на точность движения 
платформы ранее не исследовалось. В рамках данной работы будет исследовано влияние неточностей установки 
меканум-колёс на точность позиционирования всенаправленной платформы и её кинематику. 

Для достижения поставленной в статье цели, то есть для оценки ошибок позиционирования, возникающих из-за 
конструкционных погрешностей в установке меканум-колёс, и анализа их влияния на точность движения 
всенаправленной меканум-платформы, необходимо построить математическую модель движения такой платформы 
с учетом конструкционных погрешностей в установке меканум-колёс. Полученная математическая модель позволит 
исследовать их влияние на точность позиционирования мобильной платформы. Эта модель может быть применена 
для формирования управления, компенсирующего погрешности позиционирования платформы. 

Материалы и методы. Описание конструкции робота. Как основа для данного исследования 
рассматривается модель робота в виде прямоугольной платформы с четырьмя меканум-колёсами при его 
движении по горизонтальной плоскости (рис. 1). При этом учитываются погрешности, связанные 
с неортогональностью осей вращения колёс и продольной оси платформы. 

Для описания кинематики робота вводится подвижная система координат Сху2 , связанная с платформой 


робота и началом в геометрическом центре платформы. 


При математическом моделировании используются следующие обозначения: Ай —щ половина расстояния 
между осями колёсных пар (при ортогональных осях колёс); / — половина ширины платформы; х, и у, — 
абсцисса и ордината точк и О, крепления привода 1-го колеса к платформе; О, — геометрический центр 1-го 


меканум-колеса; И,,Г. — продольная и поперечная скорости геометрического центра платформы; 5 = [2.0 — 
длина вала привода; © — угловая скорость верчения платформы; о, — угол между осью Сх и осью вращения 


1-го колеса (параллельной вектору ДО, ) (рис. 1). Номера колёс платформы указаны на схеме на рис. 1, 


а эллипсами на колёсах показано расположение контактирующих роликов каждого колеса. 


Рис. 1. Кинематическая схема платформы 


Описание метода построения кинематической модели. Движение платформы будем описывать в рамках 
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неголономной постановки задачи, используя упрощенную модель меканум-колёс: колесо моделируется диском 
радиусом у ‚› а точка контакта может проскальзывать относительно направления, перпендикулярного оси 
контактирующего ролика. При этом в случае отсутствия погрешности в расположении осей колёс существует 
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связь между угловыми скоростями вращения колёс вокруг собственной оси Фф, и скоростями платформы 


Р.Р, [15]: 
, 1 
Ф; = + апб, —(у,-х ‘пб, )© |, | 


где 6, —_ угол между направлением оси контактирующего ролика и плоскостью меканум-колеса. 


Для исследования влияния конструкционных погрешностей необходимо описать зависимость между 
угловыми скоростями колёс ф, и скоростями платформы Г,,Г,,© при учете погрешностей, связанных с 


расположением осей колёс. Параметры конструкции платформы и расположения колёс приведены в таблице 1. 


Таблина 1 
Параметры конструкции платформы 
1 х, я 5, о; а, 
1 [ —45° 90° + До, 5 
2 —/ 45° —90° + Аа, —5 
3 —й [ 45° 90° + Да, 5 
4 —й —/ —45° —90° + Да, —5 
Здесь Ла, — угол, характеризующий  неортогональность осей колёс и продольной оси 


платформы Сх; 4, — проекция вектора Д.О; на собственную ось вращения {-го колеса. 


Результаты исследования. Описывая кинематику меканум-платформы в рамках неголономной механики, 
как это было сделано в работах [3—5], при отсутствии конструкционных погрешностей ( До, =0), определяем 


условие непроскальзывания точки контакта вдоль оси контактирующего ролика: 

Г. ‚ с0$(5, + Ла, )+ И, зип (5, + Ла, ) = 0, (2) 
где Г’. „, т’, —_ проекции скорости точки К; контакта 1-го колеса с опорной поверхностью, которые 
определяются из векторного уравнения: 


И; =И. +0, хСО, +в, хОК,. (3) 


Здесь Г» а“ ’., о] —щ вектор скорости точки К, контакта Г-го колеса с опорной поверхностью в 
проекциях на оси подвижной системы координат Схуг; И; =[И, И, 0] — вектор скорости геометрического 
центра платформы в проекциях на оси Схуг; О, =[0 0 о] — вектор угловой скорости платформы в 
проекциях на оси Сху2; СО, = СО, +О.Б, = [х,-а,зтАа, у, +4, с0$ До, 0] — радиус-вектор, соединяющий 


геометрические центры — платформы С и Гго колеса О, в  проекциях на оси Схур; 


. : . т КИ а т 
®, =[ф.зшАа, Фф,с0зАо, ©] — вектор угловой скорости в проекциях на оси Сху2; ОК, =[0 0 -К] 
радиус-вектор, соединяющий геометрический центр #-го колеса О, и точку контакта ролика К,. 


Переписывая векторное уравнение (3) и используя кососимметрические матрицы для вычисления векторных 
произведений, перепишем уравнение (3) в матричной форме: 


"| [7] [0 -О 9 [х-а,зтда, 0 о ф,созДа, [0 
= ИО 0 л 0|у +4, 008 Да, |+ [Ф) 0 —ф.зшАа, || 0 |. (4) 
0 0 ооо 0 —Ф, со5Да, ф,‚ зшАа, 0 —^ 


После проведения вычислений: 

7, ЕЙ, -О (у, +4, созДа, )- Еф, созДа,, 
Г 5 
к, =У. +0 (х,-4, зтАа, )- Еф, зш Аа... ы 


После переноса параметров платформы из таблицы 1 в уравнения связей (2) с учетом (5) их можно 
представить в виде: 
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Ф, и, -И, - (1+1 |соз Аа, + И, +, + (1-19 ] т Аа, то, 


ф, р. +7, +(1+1О ]соз Ла, : [А 7. -(# Го] та, +0, 
Г. в _ (6) 


ФИ +1, (6+9 воза, = [и 7. +(# о] пла, -.0, 


| 1 1 : 5 
ф, = НЙ И. + (1+1 | сз Да, ай +, - (1-1) о | т Аа, а 
Отметим, что в при отсутствии погрешностей крепления осей двигателей колёс (ла, = 0) уравнения (6) 
принимают вид: 


. меа 1 : 4еа 1 
фи й [И -,-(#+15], фи й и +7, +(1+10 |, 


1 


(7) 
Чей _ 1 ме _ 1 
фа = р +7, -(1+)0], фм = Ее =. +(#+1 а]. 
И совпадают с известными результатами [6]. 
В матричной форме уравнения (6) и (7) можно представить в виде: 

9=Н (Аа, )%, (8) 
9= Ном, (9) 
где п=[И, И, о] — вектор псевдоскоростей платформы; 4=[$, ф, $, $. — вектор угловых 
скоростей колёс; Н (ла, ) — матрица коэффициентов при псевдоскоростях в уравнениях (6); 


аН,=Н(Аа, =0). 
При формировании управления на кинематическом уровне по желаемым значениям псевдоскоростей 
платформы Я“ рассчитываются значения угловых скоростей платформы по формуле. 

4 =Ной, (10) 
где верхний индекс а указывает на желаемые (программные) значения скоростей. Полученные значения 
угловых скоростей подаются на вход системе управления двигателями колёс. 

Однако при наличии неточностей крепления осей двигателей колёс по результатам отработки программных 
значений угловых скоростей колёс будет реализовано движение платформы в виде: 
&=(Н"(ла,)Н(Аа,)} НТ (да,)4* =(Н` (да,)Н (Аа, НТ (до, ) Но". (11) 
По найденному выражению можно оценить погрешность отработки программного движения платформы, 
обусловленную наличием неточностей крепления осей двигателей колёс, по формуле: 


Аи=Я- д“ = |5" (ла,)Н(да,)) " Н` (Аа, ) Ну -в, №, (12) 


где ЕЁ, — единичная матрица порядка 3. 
Поскольку псевдоскорости платформы связаны с координатами геометрического центра в неподвижной 
системе отсчета соотношениями [15], 
Хе ЕЙ, соз\у — И, эту, 
Ус ЕЙ, зш\у + И, со5\у, (13) 
ф=0, 
то, учитывая погрешности псевдоскоростей, можно переписать указанные выражения в виде: 


Хе = № с05 Ду — у“ зт Ау +ДИ, с03 (м + Ау )- АУ, эт (м - Ау), 


Ус = Ха зщ Ам + у“ с0$ у + АИ, т (/ + ду) + АТ. с0$ (и + Ау), (14) 
у=\“ + Ау, 
где Хи =И/ созлу" И ту“, у“ = эшли" +11 соз\у" — законы изменения скоростей геометрического 


центра платформы В неподвижной системе отсчета, соответствующие программному движению; 
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1 
ди =“ (0)+ [© (1)441 — закон изменения угла курса платформы, соответствующий программному 
0 


движению; А\у = Глоба, — погрешность угла курса. 
0 


С учетом введенных обозначений получим формулы для ошибок позиционирования в виде: 
Ахс = с (с05 Ду —1) - ус зтАу + АИ, с0$ (м + ду) АХ. йа (и + ду), 
Ду = эт Ау + у (сз Ау -— 1 +АИ, зщ (м + ду) + АТ. с05 (и + Ау). __ 
В качестве числового примера для оценки погрешностей рассмотрим параметры робота КОКА уопВо [16]: 
р=0.235м, [=0.14м, №=0.05м, 5=0.01м. 


В качестве примера примем погрешности крепления осей двигателей колёс равными 
До = Да, =1°, До, = Да, = 2°. 
Ниже приведем результаты расчета погрешностей скоростей платформы при простых движениях 
платформы. В качестве примеров простых движений платформы рассмотрим: 
— поступательное движение вперед: И, =1м/с, Г. =0м/с, О=Орад/с; 


— поступательное движение влево: И, =Ом/с, ГИ, =1м/с, О=0рад/с; 


— вращение вокруг вертикальной оси, проходящей через геометрический центр платформы: 
„=0м/с, Г. =0м/с, О=1Трад/с. 

Погрешности псевдоскоростей платформы при указанных движениях рассчитаем по формулам (12): 

— движение вперед: 

ДИ, =1.35.10* м/с, ДИ, =-8.73.10° м/с, ДО = 6.8.10 рад/с; 

— боковое движение влево: 

ДИ, =8.72.10° м/с, ДИ, =-3.81-10“ м/с, ЛО =1.46-10“* рад/с; 

— вращение вокруг вертикальной оси, проходящей через геометрический центр платформы: 

ДИ, =-2.49.10* м/с, ДУ, =-4.22.10° м/с, ДО =1.97.10* рад/с. 

Полученные значения погрешностей псевдоскоростей являются значительными для мобильных роботов 
с меканум-колёсами, поскольку приводят к отклонениям по координатам порядка 10 см за 2 минуты движения 
робота. Это делает недопустимым функционирование меканум-платформы в автономном режиме в случае 
отсутствия компенсации указанных погрешностей позиционирования. Отметим, что наличие ненулевой 
погрешности угловой скорости вращения платформы приводит к нелинейному росту ошибки 
позиционирования платформы (формулы (15) при рассмотрении управления на кинематическом уровне. 

Данные результаты позволяют сделать вывод о необходимости учета неточностей крепления осей колёс при 
построении системы управления. Отметим, что, поскольку погрешности крепления осей двигателей колёс 
являются неизвестными величинами, то для их определения требуется разрабатывать методики идентификации 
параметров математической модели движения робота [16]. 

Применение разработанной кинематической модели меканум-платформы, учитывающей неточности 
крепления осей двигателей колёс, даст возможность построить управление на кинематическом уровне 
с большей точностью позиционирования. 

Таким образом, в рамках данной работы представлена модель кинематики мобильной платформы 
с меканум-колёсами при наличии погрешности установки осей двигателей колёс. Получена формула для 
оценки погрешностей псевдоскоростей платформы при программном управлении, формируемом 
на кинематическом уровне. Проведены оценки погрешностей отработки программного движения для простых 
движений меканум-платформы (движение вперед, боковое движение влево и вращение на месте). 

Обсуждение и заключение. Полученная кинематическая модель мобильной платформы с меканум- 
колёсами может быть использована для повышения точности позиционирования за счет формирования 
управления движением платформы, компенсирующего влияние неточностей установки колёс. Для 
использования полученной модели в задачах управления меканум-платформой необходимо разработать 
в будущем методики идентификации параметров математической модели. 

В качестве направления для дальнейших исследований в этой области можно назвать разработку 
математической модели динамики меканум-платформы с учетом погрешностей установки осей двигателей колёс. 
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4 врр-161@уапдех.ги 


Аннотация 


ЕШМ: ЛВОЕЗО 


Введение. При численном решении задач теории упругости в трехмерной постановке методом конечных 
элементов применяются конечные элементы (КЭ) в форме параллелепипедов, призм и тетраэдров. Обычно 
построение матриц жесткости объемных КЭ базируется на принципе изопараметричности, суть которого 
состоит в использовании для аппроксимации геометрии и перемещений полиномов Лагранжа. В расчетной 
практике наибольшее распространение получили так называемые полилинейные изопараметрические КЭ 
с линейным законом аппроксимации перемещений. Главный недостаток данных элементов кроется в эффекте 
«юос_те» («запирания») при моделировании изгибных деформаций. Причем погрешность численного решения 
существенно возрастает в случае, когда конструкция, по сравнению с обычными деформациями, претерпевает 
значительные смещения как жесткое целое. Многолетний опыт решения задач механики деформируемого 
твердого тела методом конечных элементов показал, что существующие объемные КЭ обладают медленной 
сходимостью при моделировании изгибных деформаций пластин и оболочек. Цель настоящего исследования 
состоит в построении на основе метода двойной аппроксимации матриц жесткости полилинейных объемных 
КЭ повышенной точности, позволяющих учитывать жесткие смещения. 

Материалы и методы. Для построения матриц жесткости объемных КЭ применен математический аппарат 
метода двойной аппроксимации, суть которого состоит в раздельном представлении функций распределения 
перемещений и деформаций внутри элемента. Хранение и обработка результирующей системы уравнений 
реализованы в алгоритмических терминах разреженных матриц. Разработка программного обеспечения 
и проведение вычислительных экспериментов осуществлены с использованием 64-х разрядной вычислительной 
платформы Мисгозой У15иа1 Зо 2013 и компилятора пе® РагаПе!| Зию ХЕ 2019 со встроенным текстовым 
редактором Пие!® У15иа1 Еогпгап Сотрозег ХЕ 2019. Визуализация результатов расчетов выполнена с помощью 
дескрипторной графики пакета компьютерной математики МаЙаб. В качестве тестового образца использован 
объемный восьмиузловой КЭ $ОГЛО185 программного комплекса АМЗУ$ Месратса|. 

Результаты исследования. Разработано математическое и программное обеспечение для исследования 
напряженно-деформированного состояния массивных конструкций при различных видах внешнего 
воздействия. На тестовых примерах с известными аналитическими решениями выполнена верификация 
авторизированного пакета прикладных программ. Показано, что построенные КЭ по точности удовлетворяют 
основным требованиям, предъявляемым к конечно-элементному моделированию пространственных задач 
теории упругости. 

Обсуждение и заключение. Проведенное тестирование разработанного математического и программного 
обеспечения показало, что построенные на основе метода двойной аппроксимации конечные элементы успешно 
конкурируют с аналогичными объемными элементами $ОГ/ЛО185 программного комплекса АМЗУ$ Месвашса|. 
Предлагаемые элементы могут быть интегрированы в отечественные импортозамещающие программные 
комплексы, реализующие метод конечных элементов в форме метода перемещений. 
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Гайдсуров П.П. и др. Применение метода двойной аппроксимации для построения матриц жесткости объемных конечных элементов 


Введение. При конечно-элементном моделировании напряженно-деформированного состояния массивных 
тел используются объемные конечные элементы (КЭ) в форме параллелепипедов (гексаэдров), призм 
и тетраэдров, построение матриц жесткости которых выполнено, как правило, по изопараметрической 
технологии [1-5]. Вместе с тем известно, что полилинейные изопараметрические КЭ при использовании 
однослойной схемы неудовлетворительно моделируют изгибные деформации даже при существенном 
сгущении сетки [6, 7]. Суть данной проблемы состоит в эффекте «/юс®ие» («запирания») элемента вследствие 
так называемой деформации «ложного сдвига» [8, 9]. Для «улучшения» изопараметрических КЭ используют 
аппарат несовместных элементов, созданных путем введения дополнительных неузловых степеней свободы или 
вспомогательных аппроксимирующих полиномов [8]. Вместе с тем наиболее эффективным способом решения 
проблемы «заклинивания» КЭ является применение моментной схемы метода конечных элементов, 
теоретические основы которой были разработаны А.С. Сахаровым [7]. В последующем данный подход получил 
название метод двойной аппроксимации (МДА) [6]. Концептуально МДА базируется на раздельном 
представлении функций распределения перемещений и деформаций внутри элемента. Целью настоящего 
исследования является построение на базе МДА и тестирование новых объемных полилинейных КЭ, 
позволяющих моделировать поведение различных конструкций при различных видах внешнего воздействия. 

Материалы и методы. Рассмотрим семейство объемных КЭ, состоящее из восьмиузлового 


и шестиузлового элементов, в глобальных декартовых осях 2,, т=1,2,3 (рис. 1). Геометрию 


т 


и перемещения КЭ представим в следующем виде: 


п. п", 
= (К) р = (К) 
2 =>. = ФСО, х,,х,) Ия =У и ФС х,,х,), 
К=1 


К=1 
ую 


т ›? т 


где 2 — узловые координаты и перемещения; ф,(х,,х.,х.) — «функции формы», представляющие 
собой произведение одномерных полиномов Лагранжа первой степени; х,, х›, х, — локальные в общем 
случае неортогональные координаты КЭ; и, — число узлов элемента. Для базового восьмиузлового элемента 


(рис. 1 а) п, = 8 «функции формы» определяем по формуле: 
1 3 

Фи хх) = Ч+рих,), (1) 
г=1 


здесь р,„, — координаты узлов в локальных осях. Значения Р,‚ задаем в форме матрицы: 


РК 


РЕ ЕЕ Е 


7 


ИА 
Рис. 1. Объемные КЭ: а — восьмиузловой; 6 — шестиузловой 


Зависимость между ковариантными компонентами тензора деформаций в локальном базисе и 
перемещениями в глобальных осях имеет вид [7]: 


1 
5; == ( 2т,] ии: 2: а 
2 
(суммирование по повторяющемуся индексу) 
где 2„,=02„/дх,; и„,=ди„/дх,, 5=Ь]. 
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Связь между вектором деформаций {=} и вектором узловых перемещений {и} представим в матричной 


форме: 
{}=ФИ из, 


где блочная матрица [2] = [2 [Р1,..[В,, | 


(6хи,) 
субматрица 
[В], =[{2°}, (ро, 12°] ‚К =12,...п,. 


(6х3) 
Выражения для вектор-столбцов рассматриваемых КЭ имеют следующий вид [10]: восьмиузловой элемент 
(рис. Та): 


Ри |2 ил + (2и12 Е Рик )хо + (2илз +21 Ре) Хз г 


+(2 „23 +212 Рзк 12 из РЖ +21 Рок Ри) Хх Хз | 


Рэк | т? о +22 Ри )х, на +2 ш,2 Рзк)хз + 


+2 ила +2 12 Рзк 12 из Ри 122 Ри Рзь)х! Х: | 


Рзк |2 из (в +2 из Ри )х а +2 из Рок )х, + 


1 Ц 
(т) р 25 о р 5 я 
ни #23 а Ри) хи Х. | ь (2) 


211 Рок 22 Ри + (2 из Рэь 1 Рок Рук из Ри + 


„Ри Ри) 


211 Рук т2из Ри нь Рак +21 Рок Рк + 2тз Ри + 


+2 3 Ри Рэк)х, 


2 то Рук Т2тз Рок ЕР +22 Рик Рак +2 из Ро ыы 


+2 из Рак Рэь)х 


шестиузловой элемент (рис. 1 6): 


| 
з 
5: 
12 
т 
+ 
— 
м 
15 
=) 
С 
+ 
| 
з 
5: 
м 
| 
— 
< 
| 


2 из Физ + (Рив Физ + Раз ф 3 ре 
1 о Фа Ра Фен ) 


{р М ЕЕ - 
, = а Е ыы 
8 | 21 Фл2 +22 Фу: + (из Фи ил Физ + 


Здесь введены обозначения: 


г _ 92 ое О’я. 
та > бтаВ — 5 
х 
9х. ХЕХ =хз=0 дх.д в ХЕХ =Хз= 
ы 9°2„ 
т,123 
дх, дх, дх 
1 й з жехо=хз=0 
- _9ф 9 а 
а ав 5 Ох ; о, В =Ъ2,3. 
о. © В ХЕХ =Хз=0 


ХЕХ =Х. =0 
Выражения для «функций формы» Ф,(х,х,,х,) шестиузлового КЭ, полученные на основании 


полинома (1) с помощью принципа «вырождения», имеют следующий вид: 


| 
а 


. 1 
Ф, о 


Гайдсуров П.П. и др. Применение метода двойной аппроксимации для построения матриц жесткости объемных конечных элементов 
[ния 


й 1 
Ф: 5 [+ ры) + рых) + рззх,) + 1+ ра) + рых) + рых) 
ь 1 
Фф. Рае ыыы 
Е 1 
ф; Рю Ва) ; 


— 
ф. = [И+рух) + рых.) + рух) + + рых) + рых.) + рых.) 


Формулы (2) и (3) являются основой для построения матриц жесткости рассматриваемых КЭ. 
Соответствующее программное обеспечение разработано на базе вычислительной платформы Мисгозой У15па1 
Эшаю и компилятора П\е® РагаПе| Зо ХЕ с встроенным текстовым редактором Пие® У15йа| Еоггап 
Сотрозег ХЕ. Процессы хранения и обработки глобальной матрицы жесткости реализованы в терминах 
разреженных матриц [11]. Для визуализации результатов расчетов использована дескрипторная графика 
компьютерной системы МайаБ. 

Результаты исследования. Исследование точности и сходимости разработанного конечно-элементного 
алгоритма выполнили на тестовых примерах, имеющих аналитическое решение. В тестовых примерах 
приведены численные решения, полученные с использованием разработанных элементов и аналогичного по 
размерности элемента $0110185 программного комплекса АМЗУ$ Меспашса| [5, 11]. Ниже приведены 
примеры, подобранные таким образом, чтобы в них присутствовало сочетание изгибных деформаций и жестких 
смещений КЭ. 

Пример 1. Разрезное кольцо, жестко закрепленное в одном сечении и загруженное сосредоточенной силой 
на свободном конце. Расчетная схема кольца показана на рис. 2. Исходные данные: средний радиус КЁ = 0,2 м; 
размеры поперечного сечения 0,2х0,2 см; модуль упругости Е= 10" Н/м?; коэффициент Пуассона у = 0,3; 
сосредоточенная сила Р = 10 Н. 

Прогиб в точке приложения силы по теории кривых стержней (точное решение) составляет [7]: 


Ед В? 10.3,14.0,23 
= ЧЕ 


0,00302 м. 


ЕЛ 1.10'.8. 333-108 


Рис. 2. Расчетная схема разрезного кольца 


Результаты сходимости представлены в таблице 1. 
Таблипа 1 
Результаты сходимости для разрезного кольца 


1/6, м/% 
Сетка 
МДА АМЗУ$ 
2х2х32 —0,00170 / 44 —0,000433 / 86 
2х2х64 —0,00292 / 3,3 —0,00125 / 59 
2х2х128 —0,00293 /3 —0,00251 / 17 


В таблице 1 в числителе приведено значение прогиба в знаменателе — относительная погрешность б. 


Механика 


369 


Бйр://уезииК-доп$а.ги 


370 


Афуапсеа Епатеетия Кезеагсй (Коз1оу-оп-Ооп). 2023;23(4):365-375. е155М 2687—1653 


Визуализация перемещений и ‚ , полученная на базе МДА и АМ$У$ на сетке 2х2х64, представлена на рис. 3 


и 4. Отметим, что поле вертикальных перемещений, полученное с помощью АМЗУЗ, не отражает зоны 


с экстремальным значением и =-0,00314 м, показанной на рис. 3 стрелкой. 


и,,м 
-0.003138 
-0.002456 
^^ -0.001793 
-0.001121 
-0.0004484 


Рис. 3. Картина распределения И (МДА) 


и, м 
. 001346 
. 001186 
. 001026 
уе в х а . 865Е-03 
7 . 705Е-03 
. 545Е-03 
. 385-03 
.225Е-03 
. 646-04 
. 956-04 


Рис. 4. Картина распределения И, (АМЗУ$) 


® 
> 


В данном примере особенно отчетливо прослеживается важность учета жестких смещений. 

Пример 2. Квадратная пластина, жестко защемленная по контуру и нагруженная равномерно 
распределенной нагрузкой. Исходные данные: длина стороны а= 1 м; толщина й= 0,01 м; модуль упругости 
Е = 10$ Н/м?; коэффициент Пуассона у= 0,25. 

Точное значение прогиба в центре пластины определяется по формуле [12]: 

4 
94" 


и Б 


2 


ЕЮ 
тде а= 0,00126; Р= 12а) — цилиндрическая жесткость; 4= 0,00888889 Н/м? — интенсивность 
—у 
распределенной нагрузки. Точная величина и, (в метрах) равна коэффициенту %. 


В данном примере рассматривалась и часть пластины с учетом условий симметрии. Результаты сходимости 


в виде графиков И, —5 для однослойной и двухслойной моделей представлены на рис. 5 и 6. Здесь и далее под 


слоями подразумевается разбивка пластины на КЭ по толщине. На этих рисунках значения параметра 5 
1,2,3,4 соответствуют сеткам: 4х4, 8х8, 16х16, 32х32. Приведенные графики отображают результаты 
решения, полученные с помощью АМЗУЗ (линия 1) и МДА (линия 2). Горизонтальная линия, обозначенная 
цифрой 3, соответствует точному решению. 


Гайджуров П.П. и др. Применение метода двойной аппроксимации для построения матриц жесткости объемных конечных элементов 


Рис. 5. Графики И, -5; Рис. 6. Графики И, -5; 


однослойная схема двухслойная схема 


Из рис. 5 следует, что при однослойной схеме разбивки значения относительной погрешности на сетке 
32х32 составляют: $5ОГЛО185 — 16 %; МДА — 10,5 %. При использовании двухслойной схемы (рис. 6) на сетке 
32х32 имеем: 5011185 — 36 %; МДА — 2,8 %. 

Патч-тестирование КЭ выполнено для схемы разбивки 16х16х2 с искажением сетки (рис. 7). Результаты 


патч-теста в виде картин распределения прогибов и. приведены на рис. 8и9. 


Рис. 7. Схема разбивки пластины для патч-теста 
и ,м и ,м 
-0,001263 
су -. 958Е-04 
7 -. 851Е-04 -0,001010 
Е -. 745Е-04 
-. 639Е-04 -0,0007575 
И 
. 532Е-04 
СУ -. 426Е-04 -0,0005050 
| 
-. 319Е-04 
Е В ы -0,0002525 
. 213Е-04 
Г 
-. 106Е-04 
И 0,0 
Рис. 8. Распределение и, ($0Г10185) Рис. 9. Распределение и, (МДА) 


Как видно из рисунков, искажение сетки при использовании 30110185 приводит к более заметной 
асимметрии поля и, ‚ чем при использовании МДА. При этом значение максимального прогиба для КЭ МДА 


и, = 0,001263 м совпадает с точным решением. 


Пример 3. Круглая пластина, жестко защемленная по контуру и нагруженная равномерно распределенной 


нагрузкой. Радиус и толщина пластины: А = 1м; Й= 0,01 м. Механические константы аналогичны данным 


примера 2. 
Точное значение прогиба в центре пластины определяется по формуле [12]: 
_ 1 да‘ 
‘р 


При значении интенсивности распределенной нагрузки 4 = 0,00888889 Н/м? значение и,= 0,01563 м. 
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В процессе тестирования были использованы два варианта дискретизации 4 части пластины (сектора) на КЭ. 
В первом варианте три стороны сектора делились на равное число отрезков. Второй вариант базировался на 
радиальной регулярной схеме разбивки сектора. При этом число элементов вдоль радиуса и циркульной части сектора 
принималось одинаковым. Рассмотренные варианты дискретизации сектора для сетки 32х32 приведены на рис. 10. 


Вариант 1 Вариант 2 


хх 


Рис. 10. Варианты дискретизации сектора круглой пластины на КЭ: 1 — равномерная разбивка трех сторон сектора; 
2 — радиальная регулярная разбивка 


Результаты сходимости в виде графиков и, —5 для дискретизации 4 части пластины по варианту | при 
однослойной и двухслойной схемах разбивки представлены соответственно на рис. 11, 12 и рис. 13, 14. 


1 2 3 4 
|| И .9 1 2 3 45 


-0,002 -0,002 
-0,004 
-0,006 
-0,008 
-0,010 
-0,012 


-0,014 


-0,004 
-0,006 
-0,008 
-0,010 
-0,012 
-0,014 


и, М и, М 
Рис. 11. Графики и, - 5 для варианта 1; Рис. 12. Графики и, - 5 для варианта 2; 


однослойная схема однослойная схема 


На этих рисунках значения параметра 5 1, 2, 3, 4 соответствуют сеткам 4х4, 8х8, 16х16, 32х32. Как и в 
примере 1, линия 1 соответствует решению на базе 30110185, а линия2 — МДА. Горизонтальная линия, 


обозначенная цифрой 3, соответствует точному решению. 

Из приведенных графиков следует, что элемент, построенный по моментной схеме на сетке 32х32х2 
вариант 2, имеет относительную погрешность 4 %. 

Картины визуализации поля распределения вертикальных перемещений и, для МДА и $0185 при 


радиальной схеме разбивки '4 пластины (сетка 32х32х2) показаны соответственно на рис. 15 и 16. 
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Рис. 13. Графики и, - 5 для варианта 1; Рис. 14. Графики и, - 5 для варианта 2; 
двухслойная схема двухслойная схема 
и, м и, М 
-. 004536 
в. МИЯ _ соо 
-0,01202 — -. 003528 
== -. 003024 
-0,009016 __ 00252 
-0,006011 С | -. 002016 
С] -. 001512 
о: 7 -. 001008 
0,0 Г] -. 504Е-03 
. 410Е-06 
Рис. 15. Распределение и, (МДА) Рис. 16. Распределение и, ($0110185) 


Из приведенных рисунков видно, что несмотря на качественное совпадение картин и,, относительные 


погрешности для максимального прогиба составляют: МДА — 3,8 %; $0110185 — 71 %. Такая значительная 
погрешность при использовании 30110185 объясняется тем, что разработчики использовали для 
аппроксимации геометрии и перемещений шестиузлового КЭ «функции формы» аналогичные функциям, 
примененным для восьмиузлового элемента, т. е. без принципа «вырождения» [7]. 

Обсуждение и заключение. Построенные на базе метода двойной аппроксимации матрицы жесткостей 
объемных полилинейный конечных элементов позволяют моделировать напряженно-деформированное 
состояние строительных конструкций произвольной геометрии при различных видах внешнего воздействия. 
Принципиальным отличием предлагаемой концепции от ранее известных конечно-элементных технологий 
является то, что перемещения в данном случае задаются в глобальных координатах, а компоненты тензора 
деформаций определяются в местных, в общем случае, неортогональных осях. 

На тестовых примерах показано, что объемные конечные элементы, построенные по методу двойной 
аппроксимации, обладают устойчивой сходимостью и успешно конкурируют с элементом аналогичного типа 
50110185 вычислительного комплекса АМЗУ$ Меспашса|. 

Разработанное математическое обеспечение может быть внедрено в отечественные импортозамещающие 
программные комплексы, реализующие метод конечных элементов в форме метода перемещений. 
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Преимущества сварки трением штуцеров с конической 
контактирующей формой малого диаметра 


Ю.В. Полетаев, В.В. Щепкин 54 
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Аннотация 

Введение. По свидетельствам многочисленных производственных испытаний, применение ручной дуговой 
сварки (РДС) толстостенных штуцеров малого диаметра (до 80 мм) не обеспечивает получения качественного 
сварного соединения, удовлетворяющего требованиям нормативно-технических документов атомных 
электрических станций (АЭС). Решение указанной проблемы возможно на основе разработки технологии 
сварки с оптимальным тепловложением и применения ее взамен РДС. Существующие технологии сварки 
плавлением не позволяют обеспечить оптимальное регулируемое тепловложение. Однако это может быть 
реализовано при разработке и дальнейшем использовании способа сварки трением (СТ). Поэтому цель данной 
работы заключалась в разработке технологии на основе автоматизированного способа сварки трением, 
позволяющей повысить качество сварных соединений штуцеров малого диаметра энергетического 
оборудования до уровня нормативных требований. 

Материалы и методы. Использовали штуцеры малого диаметра с конусной контактирующей поверхностью из 
низколегированной стали 10ГН2МФА. Экспериментальное исследование выполняли на машине сварки 
трением МСТ-41. Использованы методы неразрушающего и разрушающего контроля качества в соответствии с 
нормативно-технической документацией атомного энергетического машиностроения. 

Результаты исследования. Разработана методика и определены оптимальные размеры конусной 
контактирующей поверхности при сварке. Показано, что оптимальное тепловложение при сварке трением 
достигается при подготовке конусной контактирующей поверхности в диапазоне углов а = 30°—40°. 
Экспериментально отработана методика и параметры режима сварки трением моделей штуцеров малого 
диаметра. В ходе исследований получена и описана циклограмма процесса сварки трением, подтвердившая 
стадийное формирование сварного соединения за счет последовательного включения в стадию нагрева 
кольцевых участков конической поверхности соединяемых деталей. Получены результаты неразрушающего и 
разрушающего контроля, подтвердившие наличие качественного сварного соединения на уровне требований 
нормативно-технических документов АЭС. 

Обсуждение и заключение. Полученные результаты исследований могут быть использованы для разработки 
технологии сварки трением патрубков, а также изделий из углеродистых и низколегированных сталей. 


Ключевые слова: имитатор штуцера, ручная дуговая сварка, сварка трением, коническая поверхность, 
циклограмма, параметры режима, оптимальное тепловложение, механические свойства 
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Введение. Основной проблемой сварки штуцеров малого диаметра является малая протяженность сварных 
швов, что не позволяет полностью контролировать качество процесса ручной дуговой сварки. Подвод 
сварочной аппаратуры к месту сварки также не всегда технологически возможен. 

Вварные и приварные штуцеры соединяются с основной деталью угловым швом способом ручной дуговой 
сварки покрытым электродом. Вварные штуцера приваривают снаружи или изнутри с конструктивным 
непроваром. В любом случае наличие незавариваемых зазоров в соединениях, являющихся концентраторами 
напряжений, снижает работоспособность конструкции, так как может быть причиной возникновения трещин. 
Более совершенным является соединение, не имеющее незаваренный зазор. Иногда с целью получения 
гарантированного провара всей стенки штуцера и исключения возможности образования трещины от непровара 
используют удаляемую после сварки формирующую кольцевую подкладку или вваривают толстостенную 
заготовку патрубка. После сварки ее рассверливают до заданного чертежом внутреннего диаметра (рис. 1). Это 
довольно трудоемко, но оправданно в ответственных конструкциях, например, при изготовлении оборудования 
атомных электростанций [1]. Данный способ изображен на рис. 1, где изображен штуцерный узел до (а) и после 
сварки (6). Здесь: р — внутренний диаметр штуцера; й, Ь — катеты сварного соединения; 5 — толщина стенки 
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штуцера. Согласно нормативному документу, из которого заимствован рисунок, допускаемый размер 


на толщину стенки штуцера должен быть >2 мм; притупление сварного соединения — 4-1 мм; угол разделки 
кромок под сварку — 50°-=5°. 


>2 р 5 


4=1 
МЕ 


® Ай < 


а) 6) 
Рис. 1. Штуцерный узел: а — до сварки; б — после сварки и удаления корня шва 


Для приварки штуцера выполняются следующие технологические операции: точение отверстия в корпусе, 
разделка кромок на штуцере, сборка узла с соблюдением перпендикулярности, прихватка штуцера, 
многопроходная сварка, удаление корневой части сварного шва путем расточки внутреннего диаметра 
до нормативных размеров. 

Известные недостатки РДС не позволяют обеспечить стабильное качество сварного соединения 
и характеризуются большим объемом литого металла шва в ограниченном пространстве. Это приводит 
к перегревам и высоким остаточным напряжениям, что требует проведения дополнительных операций в виде 
термообработки. Стоимость ремонта оборудования в случае поломки высока. Это приводит к временным 
и экономическим затратам [2—6]. 

Анализ результатов литературного обзора и производственного опыта свидетельствуют о технической 
и экономической целесообразности отказа от использования способа РДС для получения сварных 
толстостенных соединений штуцеров и патрубков малого диаметра (до 80 мм) [7, 8]. 

В этой связи было решено исследовать влияние формы разделки кромок и способа сварки трением взамен 
РДС с целью обеспечения оптимально-минимального тепловложения и стабильного качества сварных 
соединений. Сварка трением как раз характеризуется оптимальным тепловоложением (по сравнению 
со способами сварки плавлением), при котором формируется мелкозернистая структура с низким уровнем 
остаточных сварочных напряжений и высоким уровнем механических свойств [9, 10]. 

Цель работы заключается в повышении до уровня нормативных требований качества сварных соединений 
штуцеров малого диаметра энергетического оборудования на основе разработки автоматизированного способа 
сварки трением. 

Материалы и методы. Теоретические и экспериментальные исследования выполнены на сварных соединениях 
из низколегированной стали 10ГН2МФА, используемой в качестве основного конструкционного материала для 
изготовления парогенератора, компенсатора давления и другого ответственного оборудования АЭС?. 

В работе использовали расчетно-экспериментальный метод выбора оптимальной формы соединяемых 
поверхностей при сварке. 

Разработанную технологию сварки трением реализовали при выполнении сварных соединений штуцеров 
на установке МСТ-41. 

Для оценки качества сварных соединений использовали методы и объемы неразрушающего 
(радиографический, ультразвуковой, визуально-измерительный, определение твердости) и разрушающего 
(испытания статической прочности и ударного изгиба, металлографические исследования сварного соединения) 
контроля, применяемые при изготовлении сварных конструкций АЭС. 


! Основные положения по сварке и наплавке узлов и конструкций атомных электростанций, опытных и исследовательских ядерных 
реакторов и установок ОП 1513-72. ОВГ: В@рз://9осз.сп.га/40оситеп 564412851 (дата обращения: 02.08.2023). 

? Федеральные нормы и правила в области использования атомной энергии. Сварка и наплавка оборудования и трубопроводов атомных 
энергетических установок НП-104-18. Приказ Федеральной службы по экологическому, технологическому и атомному надзору № 554 от 
14.11.2018 г. Консультант плюс. ПВГ: рз://зидасе.гаЛау//рика7-гозбеКкБпа4тога-0+-14112018-п-554-о/пр-104-18/ (дата обращения: 
02.08.2023). 
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С учетом технологических возможностей СТ на первом этапе была опробована технология сборки без 
подготовки кромок в виде «классического» таврового соединения (рис. 2 а) [11]. Параметры режима сварки 
были выбраны в соответствии с рекомендациями, изложенными в [11]. 


а) 6) 


Рис. 2. Имитаторы штуцерного узла: а — имитатор таврового соединения; 6 — соединение с предварительной ВЫТочЧкоЙ 


канавки: 41 — диаметр имитатора штуцера; @2— диаметр имитатора корпуса; А— зазор 


Результаты исследования. Контроль качества показал, что сборка штуцера и цилиндрического корпуса 
оборудования обеспечивает получение между корпусом и штуцером зазора А, величина которого меняется в 
зависимости от диаметра корпуса 4 и штуцера 41. Такая форма стыка не обеспечивает равномерного нагрева 
контактирующей поверхности, что приводит к образованию участков пластичного материала и 
проскальзыванию их с нарушением сварного стыка. 

Учитывая отрицательные результаты СТ на первом этапе, на втором этапе была предложена и отработана 
технология сварки трением с предварительной выточкой канавки на корпусе (рис. 2 6). Как показали 
результаты экспериментов, в начальной стадии сварки произошел разогрев кольцевой части, нарушение 
центровки и последующий отрыв нижней части имитатора штуцера, что в конечном итоге привело 
к формированию бесформенного сварного соединения. 

В ходе проведенных экспериментальных исследований было выявлено, что для получения качественного 
сварного соединения основания большой толщины с некоторым радиусом кривизны и толстостенного штуцера 
малого диаметра, полученного способом сварки трением, необходимо: 

— разработать конструкцию соединяемых поверхностей, обеспечивающую надежное соединение и высокие 
механические характеристики сварного соединения штуцеров малого диаметра в состоянии без последующей 
термической обработки или с термообработкой до более низких температур; 

— обеспечить оптимальный нагрев стыкуемых поверхностей за счет применения обоснованных значений 
параметров режима сварки. 

Анализ полученных экспериментальных данных свидетельствует о необходимости разработки конструкции 
соединяемых поверхностей, которые обеспечат более равномерный нагрев и формирование при СТ качественного 
соединения. Одним из вариантов получения такого соединения была выбрана коническая форма поверхностей [12—15]. 

В процессе выбора оптимальной контактирующей поверхности выполнен поверочный расчет прочности 
сварных соединений цилиндрической и конической форм (рис. 3). Расчеты для цилиндрической поверхности 
производили по формуле [16]: 

М, = (с, |.Ё.5, (1) 


где [ор] — предел прочности, МПа; 2 — длина сварного шва, мм; 5 — толщина сварного шва, мм. 
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Рис. 3. Типы сварных соединений: а — цилиндрическая соединяемая поверхность; б — коническая соединяемая 
поверхность: М— нагрузка; 1 — длина сварного шва; а — угол соединяемых поверхностей 


Расчеты ДлЯ конической поверхности производили по формуле: 


м _2]-2-5 


Ш 0, 


(2) 


Можно увидеть, что наименьшая длина сварного соединения и величина разрушающей нагрузки характерна 
для модели с цилиндрической контактной поверхностью. При переходе к конусной поверхности контакта длина 
сварного соединения и величина разрушающей нагрузки М закономерно увеличиваются. Эти характеристики 
напрямую зависят от величины угла (4, оптимальный диапазон которой по результатам расчетов представлен в 


таблице 1. 
Таблица 1 
Влияние угла наклона а на показатели прочности 

Тип соединения Угол наклона, 4,° Нагрузка №, Н 

Цилиндрический = 59 283 

30,0 118 566 

32,5 110 396 

Конический 35,0 101 609 

37,5 97 396 

40,0 94 959 


Следует отметить, что необоснованное изменение величины @а как в большую, так и в меньшую сторону 
будет приводить к снижению прочности (при уменьшении длины сварного шва) или увеличению перегрева и 
величины сварочных деформаций (при увеличении длины сварного шва). 

В этой связи предложено использовать конусную контактирующую поверхность с одной стороны — на 
штуцере, с другой — ответную конусную поверхность на корпусе (рис. 4 а). 


Металл, подвергнувшийся 
термодеформационному 
воздействию 


а) 6) 


Рис. 4. Штуцерный узел: а — форма контактирующей поверхности; б — выполненное сварное соединение 


Данная конструкция позволяет выполнять сварные соединения вне зависимости от диаметра корпуса, так 
как пластичный металл в процессе сварки равномерно заполняет соединяемую поверхность (рис. 4 6). 
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Так как разработанная конструкция обеспечивает равномерное заполнение разделки сварного соединения 
вне зависимости от формы корпуса, было предложено использовать плоскую поверхность корпуса для 
дальнейших экспериментальных исследований. Было выполнено конструирование и изготовление моделей 
штуцеров с разными углами конусной поверхности (рис. 5). 


Рис. 5. Модель штуцерного узла: а — угол конуса штуцера; В — угол конуса корпуса 


Размеры углов сопрягаемых поверхностей стыка были выбраны с учетом результатов расчетов, а также 
энергетической возможности сварочной машины МСТ-41 и ожидаемых размеров сварного соединения [18]. 

При подготовке поверхностей элементов моделей под сварку с углами конуса а = 30,0°; 32,5°; 35,05; 37,5° 
и 40,0° сварные соединения качественно формировались при обосновано выбранных параметрах режима СТ. 
При этом наиболее качественное соединение получалось с углом поверхности 37,5°. 

В ходе экспериментальных исследований была выявлена закономерность при формировании сварного 
соединения, а именно, происходил поэтапный процесс трения, нагрева и пластического деформирования по 
поверхности сварного соединения. Данные отличия от известных способов формирования соединений при 
СТ не позволяют использовать классическую циклограмму процесса сварки [11]. Это потребовало более 
детального изучения последовательности формирования соединения и описания циклограммы процесса 
в соответствии со спецификой (рис. 6) 


п | Е Е 
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— 
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— 
Тм т 
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ИА 
—— 
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Рис. 6. Циклограмма процесса СТ: 
п — угловая скорость; Ё’` — давление нагрева; М — момент трения; № — мощность тепловыделения; 
Т — температура в зоне соединения; Д| — осадка; т!-т7 — фазы цикла 
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В процессе сварки фазы ти, 12, 13 схожи с первоначальными фазами СТ цилиндрических поверхностей. 
В первой фазе происходит притирка поверхности по отдельным микровыступам, но, в отличии от класических 
способов, притирка идет не по всей поверхности, а только на соприкасаемом с конусом участке. Во второй фазе 
происходит увеличение площади контактирующих поверхностей и, в следствие этого, рост температуры. Третья 
фаза характеризуется подъемом температуры и выделением пластифицированного металла (грата). Четвертая фаза 
14 — фаза последовательного нагрева всей конусной поверхности. В данной фазе происходит частичный повтор 
процессов, протекающих во второй и третьей фазах. Отличие заключается в том, что участки, вступающие 
в процесс нагрева, уже частично подогреты. В пятой фазе 15 наступает процесс равномерного нагрева всей 
свариваемой поверхности, близкий к квазистационарному состоянию. Шестая тб и седьмая 17 — фазы торможения 
и проковки, в которых происходит процесс полной остановки вращения и сжатия. 

Реализация указанной циклограммы возможна при использовании разработанной методики расчета 
параметров режима сварки трением [17]. К основным параметрам относятся: время сварки %, давление нагрева 
Р‚, скорость вращения У и давление проковки Ри. Результаты расчетов приведены в таблице 2. 


Таблица 2 
Параметры процесса СТ модели штуцера 
Номер Давление Давление Скорость Время 
Материал 
образца нагрева, МПа проковки, МПа | вращения, об/мин сварки, с 
19-1 0,323 0,539 1200 27 
и 10-2 0,367 0,539 1200 32 
10ГН?МФА - - 
19-3 0,340 0,539 1200 30 


За счет последовательного включения в стадию нагрева кольцевых участков конической поверхности 
соединяемых деталей можно значительно уменьшить величину давления нагрева и увеличить время сварки. 

Предложенные технологические решения позволяют получить качественные соединения при меньших параметрах 
режима, что снижает тепловложение в сварном шве, обеспечивая надежное соединение и минимальный рост зерна. 
Кроме того, при таком способе сварки вступающие в нагрев кольцевые участки в месте трения имеют прослойку из 
пластичного металла из предыдущего слоя, который выступает в роли смазывающего материала. Это снижает 
коэффициент трения и, тем самым, уменьшает требуемую мощность сварочной машины. 

Авторами был выполнен контроль сварных соединений неразрушающими и разрушающими методами. 
Результаты контроля качества и механические свойства сварных соединений приведены в таблице 3. 
Результаты измерения твердости (рис. 7) и микроструктура (рис. 8) выполненных сварных соединений 
удовлетворяют требованиям нормативно-технических документов?. 

Значения твердости во всех образцах указывают на достаточно однородную структуру по всему сечению 
шва. Значений твердости, соответствующих твердости закалочных структур, не выявлено. На рис. 7 
наблюдается небольшое увеличение твердости по линии сплавления. Для выравнивания значений по сечению 
соединения рекомендуется провести термообработку — термофиксацию. 


Таблица 3 
Результаты механических испытаний, рентгенографического и ультразвукового контроля 
выполненных сварных соединений штуцеров 
ы - д а 
= Е В о о в Е 
Е а Я яз Е 5 Е — 8 2= Е 
Ы 5. а = Е в >= Е ВИ 5 
о о 2й= Е = ВИ- во Я ан ы 
з о &> я Е - ы ее 
Е. 5. на = 3 он 55 5 вр 
В ы о ое = ЕВ г во |= 
> 5 о 5 ЕН Е > > эм г © 
т 5. 5. 5. о? я 2 > 
Е Е Р 
10-1 635 595 14—17 70—78 123 251 уд | уд 
10ГН2МФА 10-2 510 425 12—16 67—71 127 195 уд Уд 
10-3 770 650 15 68—70 125 273 уд | уд 


3 Правила контроля металла оборудования и трубопроводов атомных энергетических установок при изготовлении и монтаже. НП-105-18. 
Приказ Федеральной службы по экологическому, технологическому и атомному надзору №554 от 14.11.2018 г. ОВС 
рз://9осз.зеспгз.г/Чосиглеп/прз/%00%90%[00%9Е-105-18/%20%90%00%9Е-105-18.р4Ё (дата обращения: 10.09.2023). 
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Рис. 7. Распределение твердости по ширине сварного соединения из стали 10ГН2МФА 
(ЗТВ — зона термического влияния) 
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Рис. 8. Микроструктура различных зон сварного соединения из стали 10ОГН2МФА, 
выполненного сваркой трением (х500): а — зона термического влияния; б — линия сплавления 


Обсуждение и заключение. Высокое и стабильное качество сварного соединения достигается за счет 
оптимального тепловложения при сварке трением контактирующих поверхностей путем последовательного 
включения в процесс нагрева и трения кольцевых участков конусных сварных соединений. Это способствует 
измельчению структуры металла шва и зоны термического влияния и получению высоких механических 
свойств металла сварного соединения. 

В результате экспериментального исследования установлено положительное влияние основных 
конструктивных и технологических факторов и, прежде всего, конусной формы соединяемых поверхностей на 
качество толстостенных сварных соединений штуцеров малого диаметра (до 80мм). Полученные 
положительные результаты разрушающего и неразрушающего контроля подтверждают требуемое качество и 
перспективность применения сварки трением сварных соединений штуцеров с конусной поверхностью при 
изготовлении оборудования АЭС. 
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Влияние армирования стекловолокном на механические 
свойства полиэстерных композитов 
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Аннотация 

Введение. Стекловолокна существенно улучшают качество композитных материалов, делают их легче, прочнее, 
устойчивее к коррозии и термически стабильнее. В научной и прикладной литературе активно обсуждаются 
сильные и слабые стороны конкретных композитов. При этом недостаточно исследовано влияние соотношения 
волокон и материала матрицы на механические характеристики композитов. Представленная работа призвана 
восполнить этот пробел. Цели исследования — изготовление композиционного материала на полимерной 
основе, армированного стекловолокном, а также изучение влияния весовых соотношений элементов на 
механические характеристики композита. Впервые публикуется отчет о сопоставлении характеристик 
композитов (с различным содержанием волокна) друг с другом и со сталью. 

Материалы и методы. В качестве исходных материалов использовались стекловолокно и полиэстер 
с добавлением медиатора для ускорения процесса формования. Образцы изготавливались вручную и при 
помощи стандартного оборудования испытывались на растяжение, твердость и ударную прочность. Результаты 
обобщали в виде таблиц, визуализировали в виде графиков и обрабатывали методом сравнительного анализа. 
Результаты исследования. Показаны способ создания образцов и методы их испытаний. Изыскания 
позволили установить, что твердость, прочность на разрыв и устойчивость к удару возрастают с увеличением 
процентного содержания случайно распределенного стекловолокна до 50% при 50% ненасыщенного 
полиэстера. В этом случае достигаются максимальные значения прочности на разрыв — 175,4 МПа, 
твердости — 38 НУ и ударопрочности — 1,56 Дж/мм?. Экспериментально доказана нецелесообразность 
превышения доли стекловолокна более чем на 50 %, т. к. механические свойства ухудшаются. Это объясняется, 
в частности, хрупкостью стекла, которая при нарушении пропорций передается всему композиту. Кроме того, 
при чрезмерно высоком объеме армирующих волокон смолы будет недостаточно для качественного скрепления 
элементов, что существенно снизит прочность материала. Большинство механических характеристик композита 
из 50 % полиэстера и 50 % стекловолокна лучше, чем у стали. 

Обсуждение и заключение. Доказано, что свойства композитного материала существенно зависят от содержания 
стекловолокна. Полученный композит сравнили со сталью. Выяснилось, что он обладает лучшими 
механическими характеристиками и меньшим весом. Это позволяет рекомендовать использовать данный 
материал для изготовления корпусов лодок. 


Ключевые слова: армирование стекловолокном, соотношение полиэстера и стекловолокна в композите, 
механические характеристики композита, сравнение композитных и стальных образцов 
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Введение. Полимеры привлекательны как матрицы композитов благодаря сравнительно низкой плотности, 
легкости в обработке и превосходным механическим свойствам. Высокотемпературные смолы используются 
для производства самолетов, ракет, лодок и другой техники. Арматура (стекловолокно) воспринимает 
основную нагрузку, особенно если композит состоит из волокон, рассредоточенных в слабой матрице 
(пример — углеродно-эпоксидный композит). Таким образом, прочность и жесткость материала определяются 
прочностью и жесткостью волокон [1-4]. Основные преимущества композитов, армированных стекловолокном, 
— высокие прочность и модуль упругости. Благодаря низкому удельному весу, хорошему соотношению 
прочности и веса, эти материалы превосходят металлы. Кроме того, выверенное соотношение усталостной 
прочности и веса существенно повышает устойчивость многих композитных ламинатов к усталостным 
повреждениям [5, 6]. 

Авторы [7] доказали, что механические показатели композитов, армированных волокнами, зависят от свойств 
материалов, входящих в их состав (тип, количество, распределение и ориентация волокон, пустоты). Важную роль 
играют также природа межфазных связей и механизмы передачи нагрузки на границе раздела фаз. 

Исследования [8] показали, что прочность на разрыв пластин из стеклопластика изменяется в зависимости 
от условий окружающей среды. В работе [9] рассматривались два типа композита с разным армированием 
стекловолокном: рубленым и смешанным с ненасыщенной полиэфирной смолой (стекловолокно 0/90). 
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Установлено, что композит, армированный рубленым стекловолокном, имел более высокий модуль упругости, 
максимальное напряжение и предел текучести, чем композит, армированный стекловолокном 0/90. 

Авторы научных работ концентрировали внимание на особенностях различных типов стекловолокна 
и материалов матриц. При этом недостаточно исследована зависимость качества композита от соотношения 
объемов волокон и материала матрицы. 

В представленной статье качества материалов сопоставляются по двум направлениям. Первая параллель — 
композиты с разным соотношением волокон и матрицы. Вторая — композит и сталь. Данные публикуются впервые. 

Цель исследования заключалась в изучении влияния содержания стекловолокна на механические свойства 
материала, то есть на его поведение при нагрузке. Оценивались твердость, устойчивость к разрывам и ударам. 

Материалы и методы. Для изготовления образцов в качестве исходных материалов использовали смолу 
(полиэстер) с добавлением веществ, ускоряющих реакцию, и стекловолокно (рис. 1). 


Рис. 1. Стекловолокно 


Полученные образцы испытывали на растяжение, твердость и ударную нагрузку. Аналогичные опыты 
проводились на стальных деталях. Каждый эксперимент повторялся пятикратно. От испытания к испытанию 
увеличивали процентное содержание стекловолокна, а затем брали средние значения результатов. 

В таблице 1 приведены механико-физические свойства используемого стекловолокна. 

Таблица 1 
Механико-физические свойства стекловолокна [10] 


Показатель Значение 
Плотность, кг/м? 2500 
Предел прочности, МПа 2400 
Модуль упругости, ГПа 85 
Удлинение, % 0,01 


В качестве основного материала использовался ненасыщенный полиэстер, который при комнатной 
температуре представляет собой вязкую жидкость. Его механические свойства представлены в таблице 2. 
Таблица 2 
Механические свойства полиэстера [10] 


Показатель Значение 
Плотность, кг/м? 1200 
Предел прочности, МПа 42 
Модуль упругости, ГПа 2,8 
Удлинение, % 3,2 


Вспомогательные материалы. К ненасыщенному полиэстеру добавлялись следующие вещества [11]: 

— кобальтовый жидкий материал для упрочнения (3—4 %); 

— медиатор — стимулирующее вещество, используемое с целью ускорения осаждения (1,5-3 %). 

Измерение плотности стекловолокна. Для определения плотности стекловолокна его погружали в воду 
в стандартном контейнере. Полученное значение р = 2500 кг/м? совпало с эталонным. 

Испытание стекловолокна на растяжение. Средний диаметр волокон выяснили с помощью 
микроскопа — 20,14 мкм. Рассчитав поперечное сечение одного волокна, определили их среднее число в пучке, 
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сравнивая массы волокон и пучка, из которого их взяли. Среднее число волокон в пучке оказалось 
равным 2200. Затем в зависимости от количества волокон рассчитали сечение пучка А = 0,7 мм". 

Эксперименты затруднялись обрывами в местах крепления из-за хрупкости волокон. Давление зажимов 
растягивающего устройства на концы волокон снизили следующим образом. Специальную бумагу с хорошей 
впитывающей способностью залили смолистым материалом, поместили в него волокна и оставили при 
комнатной температуре до полного высыхания. Затем бумагу разрезали и получили образцы из пучков волокон, 


готовых к испытаниям (рис. 2). 


Рис. 2. Образцы волокон на растяжение 


Размеры образцов на растяжение определялись в соответствии с указаниями международных стандартов 
А$ТМ О 2343-95! для таких испытаний. Параметры экспериментов: длина волокон между оправками — 25 см, 
длина оправок на обоих концах — 5см, скорость приложения растягивающей нагрузки — 12,7 м/с, 
лабораторная температура — 23°С. 

Результаты исследования. 

Изготовление образцов. Изготовили образцы с различным массовым содержанием стекловолокна (10%, 
20%, 30 %, 40 %, 50 % и 60 %). Для этого смесь из 5 % полиэстера, кобальта и катализатора готовилась при 
комнатной температуре. Некоторое ее количество равномерно распределяли кистью по внутренней 
поверхности формы. Затем наносили слой со случайно распределенными стекловолокнами. Далее слои волокон 
пропитывались смесью и накладывались друг на друга, пока не получался окончательный образец 


соответствующих пропорций (рис. 3, 4). 


а) 6) 


Рис. 3. Формы образцов: а — изготовление; б — готовые 
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Рис. 4. Форма стандартного образца 


Готовые образцы доставали из формы, очищали и помещали в сухую среду на 24 часа. 

Испытания на растяжение. Испытания на растяжение проводили в соответствии с 130 326 АЗТМО 638?, 
как показано на рис. 5. Образец прочно закреплялся на обоих концах растяжного устройства. К нему 
постепенно прикладывалась растягивающая нагрузка от нуля до 20 Н/мин, пока образец не разрушался. 


Рис. 5. Растяжное устройство 


Процесс растяжения разделялся на три области. В первой происходило скольжение между губками зажима 
и поверхностями образца, поэтому образцы хорошо зажимались. Во второй области из-за упругого поведения 
контактирующих поверхностей материалов наблюдалась прямая зависимость между напряжением 
и деформацией. Третья область — после разрушения образца. Оно не было ни полным, ни внезапным, потому 
что шло равномерно в зависимости от углов ориентации волокон в слоях стекломата. В результате получили 
показатели максимальной прочности на разрыв и упругой деформации образцов. 

Предел прочности при растяжении фей ‚ который выражает напряжение разрушения композиционного 


материала, рассчитывался по соотношению [12]: 


2. а 
ое =—\, 
В.е 
где К, — максимальная сила, которая привела к разрушению образца, Н; В — ширина образца, мм; 


е — толщина образца, ММ. 
На рис. 6 показаны результаты испытаний стекловолокна при растяжении. 
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Напряжение на растяжение, Н/мм? 


Удлинение, % 


Рис. 6. Поведение волокон при растяжении 


По данным рис. 6 рассчитали значение модуля упругости при растяжении. Установили его среднее значение 
Е = 80000 МПа, которое приняли для дальнейших испытаний. 

На рис. 7 и в таблице 3 показаны результаты влияния процентного содержания волокон на прочность при 
испытаниях на разрыв. 


„ 200 
Е 180 
= 160 
& 140 
8 120 
8 100 
а 80 
2 60 
|= \ 

2 40 
Я 20 


Содержание волокон, % 


Рис. 7. Влияние процентного содержания волокон на прочность образцов при испытаниях на разрыв 


Снижение значения прочности на разрыв при содержании волокон 60 % связано с тем, что полиэстера 
становится недостаточно для связывания такой массы волокон. Следовательно, силы сцепления между 
основным и армирующим материалами уменьшаются. 


Таблица 3 
Значения прочности на разрыв в зависимости от процентного содержания волокон для каждого образца 
Содержание волокон в образце, % Значения прочности на разрыв, МПа 

0 50 

10 63,3 
20 101,5 
30 145,7 
40 174,3 
50 175,4 
60 150,2 


Из рис. 8 и таблицы 4 видно, как содержание стекловолокна влияет на упругую деформацию образцов. 
Следует отметить, что с увеличением содержания волокон показатель упругой деформации возрастает 
и достигает максимального значения (2,71%) при 50%. С дальнейшим ростом массы волокон упругая 
деформация уменьшается, что связано с высокой прочностью стекловолокна на разрыв. 
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Рис. 8. Влияние содержания волокон на упругую деформацию 


Таблица 4 
Влияние содержания волокон в образцах на упругую деформацию 
Содержание волокон в образцах, % Значения упругой деформации, % 
10 1,04 
20 1,30 
30 2,32 
40 2,60 
50 2,71 
60 2,10 


Испытания на ударопрочность. В опытах задействовали устройство Шарпи в соответствии с 1$0 1793. 
Использовались образцы с размерами (80х10х4 мм) аналогичные тем, что показаны на рис. 9: 


Рис. 9. Образцы для испытаний ударопрочности 


Прочность рассчитывалась делением значения затраченной энергии на площадь сечения образца: 


Е 
В =—, 
А 
где К — ударопрочность, Дж /мм?; Е — затраченная энергия, Дж; А — площадь поперечного сечения 


образца мм?. 
На рис. 10 и в таблице 5 показано, как содержание волокон влияет на ударопрочность. 


3 150 179-1:2010. Пластмассы. Определение ударной прочности по Шарпи. Часть 1. Неинструментальный метод испытания на удар. ОВГ: 
Брз://уууууу. созит о .га/сайа!о/Оейа5/214=4569214 (дата обращения: 27.09.2023). 
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Рис. 10. Влияние содержания волокон на ударопрочность 


Ударопрочность увеличивается с повышением содержания волокон до 50 %. Это объясняется ростом силы 
сцепления в результате насыщения волокон полиэстером в композиционном материале. Если же вес волокон 
превышает 50 %, ударопрочность снижается. Это обусловлено хрупкостью стекла [7, 13]. 

Таблица 5 


Влияние содержания волокон на ударопрочность 


Содержание волокон в образцах, % Значения ударопрочности, Дж/мм? 
10 0,75 
20 1,27 
30 1,37 
40 1,45 
50 1,56 
60 0,70 


Испытания на твердость. Твердость образцов проверялась на твердомере, а число Виккерса 
рассчитывалось по соотношению: 
Р 
НУ о. 
где: Р —_ приложенная нагрузка, кН; р —_ Диаметр удара, мм. 
Испытания показали рост твердости при увеличении содержания стекловолокна до 50% (рис. Пи 


таблица 6). 
40 
> 35 
ее 
р. 
5 
8 30 
5. 
[22] 
| 
25 
20 
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Рис. 11. Влияние содержания стекловолокна на твердость 
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Таблица 6 
Влияние процентного содержания волокна на значения твердости 
Содержание стекловолокна в образцах, % Значения твердости, НУ 
10 23,0 
20 26,0 
30 33,0 
40 37,2 
50 38,0 
60 32,0 


При увеличении процентного содержания стекловолокна до 50% усиливается твердость композита. 
Волокна и полимерные цепи соединяются в случайном порядке. В результате усиливается сопротивление 
композита проникновению, то есть твердость [7, 13]. Если же содержание волокон превышает 50 %, твердость 
сокращается из-за диспропорции волокон и полиэстера (его в этом случае недостаточно для завершения 
соединения армирующего и основного материалов). 

Отметим также, что в процессе исследования доказана неточность выводов в работе [14]. Рассматривались 
полиэстерные композиты с массой армирующего Е-стекловолокна 15 %, 30%, 45 % и 60%. Исследовалось 
влияние содержания стекловолокна на такие механические свойства, как прочность на разрыв и изгиб, ударная 
вязкость. Твердость композитов оценивалась при помощи твердомера Бринелля. Результаты показали 
значительное улучшение механических свойств композита с увеличением массы стекловолокна, особенно при 
показателе 60%. Однако авторы использовали слишком большой диапазон содержания стекловолокна 
в композите, что осложняет оценку точности полученных результатов. Кроме того, не учитывались процессы 
разрушения волокон, не проверялись их свойства. 

Испытание стального образца. В соответствии со стандартом АЗТМ“ образцы изготовили из листовой 
стали АЗ6. Испытания на растяжение, удар и твердость проводились так же, как для образцов из полиэстера, 
армированных стекловолокном. 

В таблице7 приводятся результаты механических испытаний стальных образцов и их сравнение 
с образцами из композитного материала со стекловолокном 50 % (такой состав композита показал лучшие 
результаты). 

Таблица 7 
Сравнение результатов механических испытаний стальных образцов и образцов 
из композиционного материала с содержанием стекловолокна 50 % 


Образцы из композиционных материалов 
Механические испытания Образцы из стали марки АЗб 
с содержанием стекловолокна 50 % 
Прочность на разрыв, МПа 400 175,4 
Ударопрочность, Дж/мм? 0,61 1,56 
Упругая деформация, % 0,11 2,71 
Твердость, НУ 135,5 38 


Обсуждение и заключение. Итак, прочность на разрыв полиэстерного композита существенно возрастает 
с увеличением массы стекловолокна и достигает максимального значения 175,4 МПа при содержании 
армирующего материала 50%. Незначительное содержание волокон (они постепенно разрушаются) 
обусловливает слабость композита. 

Механические характеристики материала из 50 % полиэстера и 50 % стекловолокна лучше, чем у стали 
(кроме прочности на разрыв). Это позволяет использовать композит для производства лодок, т. к. в этом случае 
прочность на разрыв не основная характеристика [15]. Лодки, изготовленные из ненасыщенного полиэстера, 
армированного стекловолокном, на 75 % легче стальных. Сталь не устойчива к коррозии, что делает лодки из 
композиционного материала более долговечными. 


* АТМ (Атенсап зостейу Юг {езНпе ап4 таена15). Англ. Американское общество испытаний и материалов. Название комплекса стандартов 
и организации в США, издающей нормативы для использования в различных отраслях производства. Стандарты АЗТМ регламентируют 
химический состав, механические, физические и электрические свойства материалов, виды обработки, способы изготовления, методы 
испытаний и тестирований, а также требования к металлопрокату. ОВГ: Б#рз://ашае.га/@тесогу/{апЧагаз/азт.В | (дата обращения: 
27.09.2023). 
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Экспериментально установлено ухудшение качества композита, если содержание армирующих волокон 
превышает 50 %. Указаны две основные причины такого явления. Первая: при чрезмерной массе стекла его 
хрупкость передается всему материалу. Вторая: при непропорционально большом объеме волокон смолы 
недостаточно для их скрепления, что негативно сказывается на прочности. 


Список литературы 

1. Катаг $.М№., Катаг У.С., Китаг У.С., РгабБи М. Ехрегитега| шуезИгайоп оп Месватса| ВеВау!юг о Е- 
О]а$5 апа 5-С1азз Ефег Вешогсе4 \уиН Ро|уезег Кезт. 55 АС ГиетпаНопа! Лоигпа! о} Месйатса! Епетееттв. 
2018;5(5):19—26. № рз://401.0ге/10.14445/23488360/ЛМЕ-У515Р104 

2. ВВоуписК А.К. МесКатшса! Ргорегйез оЁ Ройутегз. Маеяа 5бчепсез & Епотеетив. 2016;1:453—461. ОВГ: 
В рз://уууу\м.е0155.пеИзатр|е-спар!егз/с05/Е6-36-01-03.ра{ (дата обращения: 06.09.2023). 

3. Ншюп М.У., бодеп Р.О., Каддолг А.5. (е4$) РаЙите Стйептоп Е фте-Кетрюгсей-Ро[утег Сотрозйез: Тйе 
Й’ота-ИЛае Гайиге Ехегстзе, 154 е4. Атзег4ат: Е1зеу1етг; 2004. 1269 р. 

4. Вау! Га, Глке Гее. (е4$) ЕБег Катюгсе4 Ро[утег (ЕКР) Сотрозйез юг тГгаягисшге Аррйсапопз. Еосиятя 
оп ипоуаНноп, Тесйпоогу трРетепаноп ап 5ияатабИиу. Ме\у Уогк, МУ: Зрипеег; 2012. 280р. 
Бр://Чо1.оге/10.1007/978-94-007-2357-3 

5. Мецег Н.Е., Соуаег Г.Е. Месвашса| Регогтапсе оЁ Ро[утег Зузетаз: ТБе Веайоп Бебмееп Згасаге ап4 
Ргорегиез. Ргоетезу т Ро[тег 5‹епсе. 2005;30(8-9);915-—938. Вир://4о1.оге/10.1016/.ргозро!утс1.2005.06.009 

6. Гемла М. (е4) НапаБоок о{ Ефег Срептизву. Пиетпанопа! Ефег 5слепсе апа Тесйпоову 5етез, ЗтА е4. Воса 
Ваюп, ЕГ: СВС Ргезз; 2006. 1056 р. № рз://401.0ге/10.1201/9781420015270 

7. Уагоа Сз., М15Ко]с71 М., Вага Г.., Глрос71 ©. Паргоуше фе Месрашса| Ргорегйез ог СЛаз5-Е1ге-ВешЕогсеа 
Ро|уемег Сотрозиез Бу МоФЙсаноп оР Ете биГасе. Маета5 апа Беяеп. 2010;31(1):185-193. 
В рз://Ч01.оге/10.1016/].та4е$.2009.06.034 

8. Арагуе! А., Саге 5., КаКезВ Р.К., Зпей Т., М1зВга В.К. ТепзИе Верау1ог оЁ С]Лаз$ Ефег Кешгсеа Р1азИс$ 
Зи есе4 ю ПШегеп Епупоптета| Соп@@оп$. м@ап Лоитпай ор Епетеетие & Меиепта5 Эсепсех. 
2010;17(6):471-476. УВГ: №&рз://порг.лизерг.гез.ш/6 И геато/123456789/10873/1/13ЕМ$%2017%286%29%20471- 
476.раЁ (дата обращения: 06.09.2023). 

9. Ебабгаа Т. АБаиПав. А З®4у оЁ Вепаше Ргорегаез ог Опзагаеа Роуезег/СЛаз5 ЕБег Вешогсе4 
Сотрозйе$. 41-Майгат Лоитпа! оф 5степсе. 2013;16(3):129—132. Вир://401.0г®/10.22401/ЛМ05$.16.3.18 

10. Оогпеё Г. Сепегашеу зиг [5 таиетаих  сотрозйез. Маме,  Егапсе: ЕСМ; 2008. 
51 р. ОВГ: В@рз://се!.Ва].зс1епсе/ЕМеЛп4ех/4ос14/470296/Шепате/МаСотрозИез.раЁ (дата обращения: 06.09.2023). 

11. МоогзвазШауай Атога Вши Моватта4. Сйагасметганоп апа Ргорегиез ор ше № Ипзашиие4 Ро[уезег 
Кезту юг Сотрозйе АррПсайоп. Сешеог, МУЗ: Зсвоо| оЁ ш4аза| ТесБпо]озу; 2007. 125р. ОВГ: 
В рз://соге.ас.иК/геа4ег/32600400 (дата обращения: 06.09.2023). 

12. УШпа2 Н. Тепз|е Эиепо Тезйпе оЁТЫт $ргау-оп Глпег Ргодис{5 (Т$Т.$) апа ЗБссгее. Тйе Лоигпа! о/ Тйе 
Зоишегп Арлсап тушие ор Митя апа Маайигэу. 2010;110(10):559-569. ОВГ: 
Бр ://уууми\м .зс1е1о.ого.та/рАЁзаииит/у1 10110/01.р4Ё (дата обращения: 06.09.2023). 

13. Агапиае Е.О., Ацапда Р.О., О1Тогаппимо О.О. МесБашса| Ргорегиез оР а Ро|уе$ег Ете СЛаз$ СотрозНе. 
1тегпайопа! Лоигпа][ о} Сотрозйе Миептаё. 2012;2(6):147—151. В рз://401.ог./10.5923/].стаета1$.20120206.06 

14. Е]-У\’ахегуа М.$., Е-ЕМмату М.Т., ГоаМаКаг 5.Н. Мескатса| Ргорегаез ОР СЛазз$ Ефег Вешгсе4 Ро[уезег 
СотрозЦез. — Пиегиайопа! Лоита[ оГ АррПе 5аепсе апа  Епетеетие. — 2017;14(3):121-131. 
Брз://401.0оге/10.6703ЛЛАЗЕ.2017.14(3).121 

15. ББепо: В.А., )иПеч-Вапоп /., Ошши $5. Сотрозйе Мажепа1$ Гог Магте АррПсайоп$: Кеу СБаПепзез юг фе 
Рабиге. ш Боок: №со1а1$ Г... Мео М., МПеПа Е. (е4$) Сотрозйе Миета!5. Гоп4доп: Зрипеег; 2011. С. 69-89. 
Бр://о1.оте/10.1007/978-0-85729-166-0_3 


ВеЕегепсе5 

1. Катаг ЗМ, Катаг У@, Китаг УС, РгабВи М. Ехрегипеща! шуезизайоп оп Месватшса| Веваутог оЁ Е-С]аз$ 
апа 5-СЛазз Ефег КешЕогсе4 уу Ро[уезег Кезш. 55АС Гиетгпапопа! Лоигпа! ор Месйатса! Епетеетпв. 
2018;5(5):19-26. № рз://401.0ге/10.14445/23488360/1УМЕ-\515Р104 

2. ВВомписк АК. МесБашса! Ргорег@ез оЁ Роутегз. Миена! бчепсех & Епетеетие. 2016;1:453—461. ОВГ: 
В рз://уууу\м.е0155.пеИзатр]е-спар!егз/с05/Е6-36-01-03.раР (ассеззеа: 06.09.2023). 

3. Ншюп МУ, бодеп РО, Каддотг АЪ. (е4$) Еайиге Стйетоп Ефте-Кетрюгсе4-Ро[утег Сотрозйез: Тре То"А- 
Ипае ЕаНиге Ехегсчве, 151 е4. Атуегдат: Е1зе\ег; 2004. 1269 р. 


'Антибас И.Р. Влияние армирования стекловолокном на механические свойства полиэстерных композитов 


4. Вам! Ллаш, ГлКе Гее. (е4$) Р\бег Кетрюогсе4 Ро[утег (ЕЕР) Сотрозиез юг тГгаягисшге Аррйсапопз. Госизтя 
оп тпоуаНноп, Тесйпоову ГтрРететаноп ап бияатаБИийу. Мем Уощ, МУ: ЭЗрипеег; 2012. 280р. 
Бр://Чо1.оге/10.1007/978-94-007-2357-3 

5. Мецег НЕ, Соуаег ГЕ. Месрашса| Реогтапсе ог Ройутег Зузетз: Тве Ве!абоп Бебмееп Зиласвте ап4 
Ргорегиез. Ргоэтезу т Ройутег 5сепсе. 2005;30(8-—9);915—938. Вир://401.0ге/10.1016/.ргогро[уп5с1.2005.06.009 

6. Ге\уш М (еа). Нап4Боок о/ Рег Сйепизту. [иетгпапопа! Рег 5счепсе ап4 Тесйпоозу 5етез, Зт4 е4. Воса 
Ваюп, ЕГ: СВС Ргезз; 2006. 1056 р. № рз://401.оге/10.1201/9781420015270 

7. Уагоа Сз, М1зКос7л М, Вайфа Г, Глрос71 @. Паргоуше фе Месватса| Ргорегиез оЁ С]аз5-Е1ге-ВешЕогсеа 
Ро|уемег Сотрозиез Бу МоФЙсайоп оР Е те биГасе. Маета5 апа Беяеп. 2010;31(1):185-193. 
6рз://401.оге/10.1016/].та4ез$.2009.06.034 

8. Арагмуе| А, Саге $, ВаКезВ РК, ЗшеВ [, М!5йга ВК. ТепзИе ВеБауюиг ог С1Лазз Ефег Вешгсеа Р1азИс$ 
Зи] есе4 ю ПШегеп Епупоптета| Соп@@оп$. мфап Лоитпа ор Епетеетие & Миепта5 5сепсех. 
2010;17(6):471-476. УВГ: № рз://порг.лизерг.гез.ш/Б И геато/123456789/10873/1/131ЕМ$%2017%286%29%20471- 
476.раЕ (ассеззеа: 06.09.2023). 

9. Ечабгаа Т АБаиПав. А За4у оЁ Вепате Ргорегаез оЁ Опзаагаеа Ро]уезег/СЛазз ЕБег Кешгсе4 
Сотрозиез. 4/-Майгат Лоигпа! оГ5ачепсе. 2013;16(3):129-132. Бир://401.оге/10.22401/ЛМ05.16.3.18 

10. Оогшеё Г. Сепегашеу зиг [е5 таиетаих  сотрозйез. Маме,  Егапсе: ЕСМ; 2008. 
51 р. ОВГ: Б@рз://се].Ва1.зс1епсе/е/п4ех/4ос14/470296/епате/Ма{СотрозИез.р4Ё (ассеззед: 06.09.2023). 

11. МоогзвазШауай Атога Вши Моватта4. Сйагасметганоп апа Ргорегиез ор ше №» Ипзашгиие4 Ро[уеег 
Кезту юг Сотро5йе АррПсайоп. Сешеог, МУЗ: Зсвоо| оЁР ш4аза| ТесБпо]озу; 2007. 125р. ОВГ: 
Брз://соге.ас.аК/геа4ег/32600400 (ассеззеа: 06.09.2023). 

12. УШла2 Н. Тепз|е Эиепо Тезйпх оЁ ТЬт $ргау-оп Глпег Ргодис5 (ТЗГ.$) апа ЗБосгае. Тйе Лоигпа! о/ Тйе 
Зоишегп Арлсап тушие ор Мття апа Маайигэу. 2010;110(10):559-569. ОВГ: 
Бр ://уууми\м .зс1е1о.оге.ха/рАЁ/]заитт/у1 10010/01.р4Ё (ассеззеа: 06.09.2023). 

13. Агапиае ЕО, А‘апда РО, О]огипи мо ОО. МесВатиса| РгорегНез оЁ а Ройуезмег Ете С]азз СотрозНе. 
1тегпанопа! Лоигпа] о} Сотрозйе Миептаб. 2012;2(6):147—151. В&рз://401.оге./10.5923/].стайема1$.20120206.06 

14. Е]-У/ахегуа М$, Е!-Е!ату М1, ХГоа№Каг ЗН. Месватиса| Ргорегаез ОГ С1азз Ефег ВКешЮгсеа Ройуезег 
Сотрозйез. — Пиегиайопа! Лоита[ оГ АррПеЯ 5аепсе апа  Епетеетте. — 2017;14(3):121-131. 
6 рз://401.оге/10.6703ЛТАЗЕ.2017.14(3).121 

15. ББепо! ВА, РиПечл-Вацоп 7, Ошши 5. Сотшрозие Маепа15 юг Маппе Арр|сайопз: Кеу СБаПепзез$ Рог Фе 
Рабие. ш Боок: №со|а1з Г, Мео М, МПеПа Е. (е4$) Сотрозйе Миета[. Гопдоп: Зритеег; 2011. С. 69-89. 
Бр://Чо1.оте/10.1007/978-0-85729-166-0_3 


Поступила в редакцию 29.09.2023 
Поступила после рецензирования 31.10.2023 
Принята к публикации 22.11.2023 


Об авторе: 

Имад Ризакалла Антибас, кандидат технических наук, доцент кафедры основ конструирования машин 
Донского государственного технического университета (344003, РФ, г. Ростов-на-Дону, пл. Гагарина, 1), 
ЭРП\-код: 7371-0223, Зсори$ ШО, Везеагс О, ОВСТО, пта4.атщураз@тай.га 


Конфликт интересов: автор заявляет об отсутствии конфликта интересов. 


Автор прочитал и одобрил окончательный вариант рукописи. 


Весеуеа 29.09.2023 
Веубеа 31.10.2023 
Ассерд 22.11.2023 


АБош йе Ашфог: 
пад ВымаКаПа Апбра$, Сап4.3с1. (Еп5.), Аззос1айе РгоЕеззог оР ше Еипдатета!$ оР МасЬтегу Оезеп 
Перагитет, Поп Зе Тесбса| ОштуегзИНу (1, Сагагт $4., Козюу-оп-Роп, 344003, ВЕ), ЗРП\-соде: 7371-0223, 


Зсориз О, Везеагс ТО, ОВСТО, ипа4.апбураз(@таЦ.ги 


Соп1сЕ о 1теге5Е зкиетет!: е амфог 40е$ по! Вауе апу соп1с{ оЁ ищегез. 


Тре ашйог ра5 геа4 ап4 арргоуе4 фе па! тапиустри. 


Машиностроение и машиноведение 


397 


Бйр://уезиК-допза.ги 


398 


'А4уапсей Епотеениз Везеагсй (Коз1оу-оп-Ооп). 2023;23(4):398—409. е155М№ 2687—1653 


ИНФОРМАТИКА, ВЫЧИСЛИТЕЛЬНАЯ 
ТЕХНИКА И УПРАВЛЕНИЕ 
МЕГОВМАТТОМ ТЕСНМОГОСУ, СОМРОТЕКВ 
ЭСТЕМСЕ АХО МАМАСЕМЕМТ 


ВУ 


УДК 004.942, 004.02, 519.622, 544.3, 544.2 Научная статья 
6рз://401.оге/10.23947/2687-1653-2023-23-4-398-409 


Автоматизация формирования математического 
описания кинетики для многостадийных химических 
реакций и численное решение прямой задачи 


Н.А. Лысенко Фа, К.Ф. Коледина`® ВБМ: РОРбО 
' Уфимский государственный нефтяной технический университет, г. Уфа, Российская Федерация 
2 Институт нефтехимии и катализа УФИЦ РАН, г. Уфа, Российская Федерация 


Ч шкИка_1узепко_2016@тай.га 


Аннотация 

Введение. Основой для исследования, анализа и математической оптимизации любого химического процесса 
является адекватная математическая модель, учитывающая кинетику объекта. Кинетический анализ 
в химической технологии является важной задачей, поскольку позволяет оптимизировать процессы синтеза 
и прогнозировать их эффективность. Многие химические процессы включают в себя несколько стадийных 
реакций. Для успешного проектирования и оптимизации необходима математическая модель, которая 
описывает каждую стадию. Создание такой модели вручную может быть трудоемким и затратным процессом, 
требующим обработки большого объема информации. Современный уровень автоматизации позволяет 
ускорить получение математического описания кинетики многостадийных реакций. В этом случае значительно 
упрощается работа с данными и уменьшается вероятность совершения ошибок. Полученная математическая 
модель может быть применена для последующего анализа и оптимизации процесса. В работе рассмотрена 
промышленная реакция каталитического риформинга бензина, занимающая важное место в современной схеме 
переработки нефти, поскольку является источником высокооктановых компонентов товарных бензинов 
и индивидуальных ароматических углеводородов. Данный процесс характеризуется участием в нем большого 
числа (до 300) различных углеводородов, изменением числа молей и неизотермичностью. Математическое 
моделирование таких процессов предполагает детализацию стадий до необходимого уровня. Рассмотрена 
детализация до 173 стадий. В такой постановке задачи автоматизация формирования математического описания 
кинетики для каталитического риформинга бензинов ранее не проводилась. Поэтому целью представленной 
работы явилась реализация эффективных численных методов и алгоритмов для автоматизации формирования 
математической модели с учётом кинетики, термодинамики и изменения числа молей. 

Материалы и методы. Математическое описание кинетики многостадийных реакций разрабатывается 
на основе закона действующих масс. Значения кинетических параметров взяты из литературных источников. 
Решение прямой задачи кинетики проводилось с применением следующих алгоритмов: метод Гира, Рунге- 
Кутты 4 порядка и метод зс1ру.одеш языка Руфоп. Концепция автоматизации реализована с помощью 
методологии ШЕЕЙ. Программное обеспечение написано на языке программирования Ру®оп. 

Результаты исследования. Создано новое программное обеспечение для автоматизации процесса 
формирования математической модели с учетом кинетики, термодинамики и учета объема реакционной смеси. 
Приведены результаты работы программы на примере каталитического риформинга бензина. Реализована 
возможность учета в модели промежуточного подогрева смеси в каскаде реакторов. Получены численные 
значения изменения температуры, соответствующие промышленным данным. 

Обсуждение и заключение. Результаты, полученные при моделировании химических превращений в каскаде 
реакторов каталитического риформинга бензина, подтвердили экзотермический характер реакции. 
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Лысенко Н.А. и др. Автоматизация формирования математического описания кинетики для многостадийных химических реакций 


Разработанный программный продукт позволяет отобразить изменения концентраций веществ реакции, а также 
изменение температуры в реакторе и может быть использован в научных исследованиях организаций ДлЯ 
анализа многостадийных каталитических процессов. Результаты моделирования кинетики реакции будут 
использоваться в последующей оптимизации условий проведения процесса на производстве. 


Ключевые слова: кинетический анализ, математическое описание, автоматизация, система дифференциальных 
уравнений, прямая задача, методы численного решения 
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АЪзфтгасе 

ттодиспоп. Те Баз1$ ог гезеагсВ, апа]у51$ ап таетайса| орипихайоп оЁ апу свеписа| ргосезз 15 ап а4еаиае 
тафетайса| то4е| фай {аКез ш®ю ассоип Фе КтеНс$ оЁ Фе об]есё. Ктейс апа|уз1$ 1$ а сва|епее ш сКеписа| 
{есбпо]огу, зшсе й аПо\мз Юг орйпиише зупе$1$ ргосеззез ап рге@сип» феш еЁ йсепсу. Матегои$ сБеплса1 
ргосез5ез шуо]уе зеуега| %асе геасНоп$. Рог зассезз В 4ез1еп ап4 орнпихайоп, а тафетайса| по4е| {Ка{ Чезстез 
еасЬ асе 15 пееде4. Сгеайт» засВ а то4е! папиаПу сап 5е Ите-сопзите ап4 созЙу, зшсе й гедитез ргосеззте а 
атгсе атоци оЁ шЮппаНоп. ТВе тодегп 1еуе| оЁ амютайоп таКез # розз1Ые 1ю ассеегае фе оМашше оРа 
шафетайса| Коппшайоп оЁ фе Кшенс$ оЁ пи$аее геасНоп$. ш 1$ сазе, \уогкте \УИй Дай 15 отеаЙу зпирИйе4, 
апа Фе ргобаб цу оЁРтаКтз пиаКез 15 гедисеа. ТБе гези п табетайса| то4е] сап Бе аррПеа Гог Ви ег апа1у$1$ 
ап4 орипитайоп оЁ Фе ргосезз. ТВе рарег сопз1егз фе ш4из@ла] геасйоп оё саа]уйс геРоплше оЁ сазоПпе, \ась 
оссир1ез ап ппроцап{ р]асе ш фе то4еги зсВете оЁ ой гейпе, зшсе И 15 а зоигсе оЁ 12Ъ-осбапе сотропеп® оЁ 
сопптегс1а! газоПпез ап@ ш@у1иа] аготайс Будгосагфопз. Т1$ ргосез$ 15 сВагасептеа Бу Фе рагистрайоп оЁа 
1атое питбег (ир ю 300) оЁ уапочз Бу4госагФопз, а сБапзе ш Фе питЬег оЁ п1о[ез, ап поп-150егтаШу ш и. 
Машфетайса| тодеНие оЁ засВ ргосеззез шуо[уез аеаШте Ше 5аез 10 Ше гедите4 1еуе|. Те дебате оРпр ю 173 
эбасе$ 15 сопз14егеа. шп 1$ зе пе, ащютайоп оЁ Фе Ююптайоп оЁа тафетайса] оплщайоп оЁ Ктейсз Юг сайауйс 
геРогилте оЁ газоПпе Ваз поё Бееп саге оп беге. ТВегеРоге, Пе ргезете4 ухо апте4 аё ипр!етепип» еЁЙесНуе 
питепса| тей о4$ ап4 а1еогИбил$ Рог ащютайпе фе БаЙ4тя оГа таетайса! по4е] {акте шю ассоипЕ Ктейс$, 
Фегло4упапл1с$, ап сВапеез ш фе пипабег оЁ п10]е$. 

Миепта[5 апа Мешоа$. ТВе тафетайса! югишШаноп оЁ йе Ктейс$ оР ти@$аее геасНоп$ \аз 4еуеюре4 оп фе 
Баз15 оЁ Фе тазз асйоп 1а\у. ТВе Ктейс рагатеегз уа[лез \уеге {аКеп Нот Шегагу зоигсез. Те тес Ктейс$ 
ргоет \уаз зо]уе4 изше а1еогИитз: Фе Сеаг те@о4, Фе Кипзе-Кийа те#о4 оЁ Фе 4 отдег, ап Фе зс1тру.одеш 
тео4 оР Фе РуФоп 1апочаое. ТБе аютаноп сопсер! \аз пар!етещеа изт фе ГРЕЕО тефодо]ову. ТБе зоЙ\аге 
уаз утШеп ш Фе Ру®фоп ргоэтатитте 1апопае. 

Кезий. А пе\ зой\маге \уаз сгеайе4 ю аиютае Фе ргосез$ оЁ Юппише а тафетайса| поде|, аКше ш®ю ассоип Фе 
Кпейс$, Шегто4упаилс$, ап фе уо[ите оЁ Ве геасноп пихее. Те ргостат гезиН$ \уеге ргезеще4 Бу Фе ехатре оЁ 
саа[уйс гео оЁ газопе. ТЬе тоде! ппретеще4 Фе розу оР4акКте шю ассоипе фе пиеппе аще Беайт» оЁ Фе 
шаге ш Фе геасюг сазсаде. Митепса| уа[лез оР{етрегаиге сВапее$ сотгезропате ю шизила| даа \уеге оМатед. 
Оё5си5яоп апа Сопсшизтоп. ТБе геза5 оМатей @гоиеВ тодеНпе сВепуса| тапзогтаНопз ш фе сазса4е о# газоПпе 
саёайуйс геРогпте теасюг$ сопЯгте4 Ве ехофегиис пабаге оЁ 1е геасНоп. Те деуеоре4 зой\уаге ргодисе ргоу14ез 
Ч1зр1ауте свапоез ш фе сопсептанопз оРтеас4ат5, аз \е[ аз {етрегаиге уапайоп$ т е геасюг, ап4 # сап Бе изе4 
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ш зс1тепЯЯс гезеагсЬ огоап1тайотз Гог Фе апа[уз1$ оР пи засе сайа|уйс ргосеззез. ТБе гези5 оЁ Фе геасйоп Ктейс$ 
то4деНие \\Ш Бе изе ш пе заб зедиет орйпутаНоп оЁ е ргосез$ соп@1оп$ ш ргодасНоп. 


Кеуумога$: Кшейс апа[уз$15, табетайса! ЮгишШайоп, ашютайоп, 4егепйа| едиайоп зубет, тес рго ет, 
сотршайопа[ пефо4$ 


Аскпо\едоетепт5. Те аифог$ мо ШКе №0 ФапК Фе геуе\мег$ Юг уаша Ме соплтеп$ па соты оше4 ю Фе 
пиргоуетепЕ оЁ йе агас1е. 


Гипато нюгшайоп. ТВе гезеагсв \уаз 4допе уИйш Фе Нате\могк оЁ Фе Пеадегз$Мр рго]есё “Оеуе!ортепЕ оЁ а 
сотршег зузет Юг апайухштяе фе КшеНс$ оЁ сВеписа| ргосе$$е$ ап Феш шиа@сгИиена орйпитаноп”, О Звае 
Рего!еит Тесппою1са1 Ошхуег$Иу. 


Ког сНайоп. Гузепко МА, Ко]е4та КЕ. Аиютайоп оЁ фе Еогтаноп оРа Мафетайса| Еогишайоп оЁ Ктейс$ Юг 
Мивз$асе СВепуса| Кеасйоп$ ап Миатегтса| Зо[аноп ю а Онесе Рго ет. Адуапсей Епэтеетия Кезеагсй (Ко5юу- 
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Введение. Изучение кинетики многоступенчатых химических процессов имеет теоретическое 
и практическое значение, поскольку позволяет раскрыть механизм реакций и осуществить контроль. При этом 
становится возможным ускорять целевые или замедлять побочные химические реакции. Кинетический анализ, 
будучи важной задачей химической технологии, способствует оптимизации процессов синтеза 
и прогнозированию их эффективности [1-3]. Многие химические процессы являются многостадийными 
реакциями. Для их эффективного проектирования и оптимизации необходимо иметь математическую модель, 
описывающую все стадии [4,5]. Ручное формирование такой модели может быть очень трудоемким 
и затратным процессом, так как подразумевает под собой обработку большого объема информации [6]. 

Объектом исследования является промышленный процесс каталитического риформинга бензина. Целевой 
продукт реакции — риформат — позволяет получать основную массу товарного бензина. Остальные 
составляющие являются продуктами изомеризации, каталитического крекинга и алкилирования. Большинство 
нефтеперерабатывающих заводов в России и мире используют установки каталитического риформинга [7, 8]. 
Каталитический риформинг бензина позволяет получать ароматические углеводороды (ксилол, толуол, бензол), 
которые затем могут быть использованы в нефтехимической отрасли [9, 10]. Однако для товарного бензина 
присутствие в продукте ароматических углеводородов ограничивается по экологическим критериям. Поэтому 
детализированная сложная кинетика, показывающая постадийное образование компонент, позволит управлять 
условиями процесса для контроля выхода продукта и ограничения содержания ароматических углеводородов. 
Процесс каталитического риформинга бензина проводится в реакторном блоке, состоящем из трех 
последовательных реакторов, между которыми установлены печи для нагрева смеси. Это является 
необходимостью, так как данная реакция — эндотермическая. При математическом описании реакторного 
блока необходимо, наряду с изменениями концентраций компонент от условного времени контакта 
катализатора и смеси реакции, учитывать изменение температуры и объема реакционной смеси. 

Автоматизация моделирования таких сложных процессов позволяет значительно ускорить получение 
математического описания кинетики многостадийных реакций, снизить вероятность ошибок и упростить 
процесс работы с данными [11]. В целом автоматизация формирования математического описания для 
многостадийных химических реакций является важным инструментом для эффективного управления 
процессами в различных отраслях промышленности, а также для достижения более точной и эффективной 
реализации химических реакций в лаборатории. 

Таким образом, целью исследования являлась разработка программного продукта, помогающего 
автоматизировать формирование математического описания многостадийных химических реакций с учётом 
кинетики, термодинамики и изменения объёма реакционной смеси, а также получение численного решения 
составленной математической модели. 

Для достижения цели исследования в данной работе решались следующие задачи: 

— разработка программного модуля формирования математического описания кинетики с использованием 
языка программирования Руфоп; 

— разработка программного модуля численного решения прямой задачи согласно составленному 
математическому описанию с использованием языка программирования Ру®оп; 

— получение результатов работы программного продукта для реакции каталитического риформинга бензина. 

Материалы и методы. Математически указанный процесс может быть описан системой обыкновенных 
нелинейных дифференциальных уравнений с начальными данными, которая и является математической 


Лысенко Н.А. и др. Автоматизация формирования математического описания кинетики для многостадийных химических реакций 


моделью. Учитывая, что каталитический риформинг бензина является неизотермической реакцией 
с изменением объема реакционной смеси, модель многоступенчатой химической кинетики может быть 
рассмотрена в следующем контексте: 


ау _\ 
—= хи, [=Ъ....М; 1 
— 2. и, 1 (1) 
40 _ 4» о 
= -;0(0)=О; 2 
г 20-0 (2) 
М М 
иж 9-5». (3) 
м © м © 
Ея Е; 
К, = К, хехр А ; в, =, хех т | (4) 
М 
а У АН (Г) 
=== $ Т(0)=Г; (5) 
ат С ст 
ри и ) 
уй 
н,(Т)=АН, (298) [С„(Т)ат; (6) 
298 
С„(Г)=а, +В Т+сТ +а. 1"; (7) 


с начальными условиями при т= 0, у, (0) = у’, Т = 763К, О =1.014; те [о,г | ПИ. 


В данной формуле реакционной кинетики имеется ряд переменных и коэффициентов, которые 


определяют концентрацию реагентов, время контакта и скорость реакции. Конкретно О иу, (концентрации 


реагентов реакции) могут быть выражены в моль/л или в долях, а переменная т, показывающая условное 
время контакта реакционной смеси с катализатором, в кгхмин/моль или кг катализатора. В формуле также 
задействованы показатели: количество стадий М (173 стадии) и количество веществ № (38 веществ), а 


Также стехиометрическая матрица У; . (. 


в = это все отрицательные элементы матрицы У; на пересечении 


7-ой строки и /-го столбца, а В, — все положительные элементы матрицы У, на пересечении {-ой строки и 
7-го столбца. Скорость /-ой стадии реакции у, может быть выражена в 1/мин или в моль/(кгхмин), а 
константы скоростей стадий А.К, — в 1/мин. Кроме этого, в формуле указаны предэкспоненциальные 
множители в. м выраженные в 1/мин, энергии активации прямой и обратной реакции Е,, Е), 
выраженные в ккал/моль, газовая постоянная К и температура Т, измеряемая в Кельвинах. 


Продолжительность реакции 1* определяется временем, которое требуется для её завершения в минутах, а 


энтальпия образования 1-го компонента АН, (Т) измеряется в Дж/моль. Удельную теплоемкость С’, (Т) 7-го 


компонента, измеряемую в Дж/(мольхК), можно выразить через коэффициенты температурной 
зависимости теплоемкости данного компонента, которые обозначены буквами а,Б,с,,а, [12]. 1 — 


порядковый номер вещества, участвующего в реакции. 

При моделировании химических превращений каталитического риформинга бензина необходимо 
учитывать болышое число индивидуальных углеводородов. Поэтому в работе применяются групповые 
компоненты в модели по выделению типов углеводородов и числу атомов углерода в структуре молекулы: 
А; — ароматические углеводороды; АСР, — пятичленные нафтены; АСН; — шестичленные нафтены; 1Р, — 
изомеры парафинов; пР; — нормальные парафины. Здесь / — количество атомов углерода в структуре 
молекулы. В таблице | приведены усредненные значения коэффициентов температурной зависимости 
теплоемкости групповых компонент, которые будут применяться в расчетах. 
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Таблица 1 
Коэффициенты температурной зависимости групповых компонент 
| ре а, Ь, с; а, 
углеводородов 
1 ПР! 34,942 —0,03996 0,19184 —0,153 
2 пПР› 28,146 0,043447 0,18946 —0,1908 
3 ПР: 28,277 0,116 0,19597 —0,2327 
4 пР4 20,56 0,2815 —0,013143 —0,09457 
5 ТР. 6,772 0,34147 —0,10271 —0,03685 
6 ПР; 26,671 0,32324 0,04282 —0,1664 
7 1Р —0,881 0,47498 —0,24797 0,06751 
8 пРе 25,924 0,41927 —0,012491 —0,1592 
9 ТР —7,123 0,58327 —0,30338 0,06802 
10 ПР7 26,984 0,50387 —0,04748 —0,1684 
И 1Рт 19,245 0,55072 —0,14055 —0,08248 
12 ПР 29,053 0,58016 —0,057103 —0,1955 
13 1Рз —3,367 0,75824 —0,38216 0,05736 
14 пПРо 29,687 0,66821 —0,096492 —0,2001 
15 1Рэ 68,581 0,44754 0,31908 —0,5118 
16 пРо 31,78 0,74489 —0,10945 —0,2267 
17 Ро —46,17 1,108 —0,70316 0,1787 
18 пПРи 125,21 0,31401 0,79137 —0,9141 
19 Ри —8,791 1,0548 —0,5778 0,1192 
20 АСНб 13,783 0,20742 0,53682 —0,6301 
21 АСНУ? 4,296 0,42716 0,21058 —0,3999 
22 АСНв —51,866 0,78827 —0,35255 —0,006855 
23 АСНо —120,89 1,2728 —1,0794 0,4035 
24 АСНуо 90,421 0,23264 0,94595 —1,057 
25 АСНи —65,48 1,1809 —0,68379 0,1458 
26 АСР —9,939 0,42528 0,012521 —0,1886 
27 АСР7 —28,514 0,58607 —0,094379 —0,1644 
28 АСР& —69,713 0.92602 —0,62526 0,2049 
29 АСРо —116,73 1,3097 —1,2439 0,5292 
30 АСР —67,341 1,0922 —0,704 0,1906 
31 АСРи —68,23 1,187 —0,7575 0,2018 
32 Аб —31,368 0,4746 —0,31137 0,08524 
33 А? —24,097 0,52187 —0,29827 0,06122 
34 Аз —17,36 0,5647 —0,26293 0,01122 
35 Ао —10,933 0,64349 —0,27829 —0,01443 
36 Аю —24,187 0,79716 —0,48265 0,1341 
37 Ап —26,717 0,91044 —0,53876 0,1203 
38 Н› 29,07 —0,0836 0,1064 0,5752 


Уравнения (1), (3) и (4) описывают изменения концентраций реагентов в процессе, а в системе (2) 
используются мольные доли для описания изменения мольного расхода смеси. В (5) рассматривается 
зависимость изменения температуры реакции от термодинамических характеристик веществ. В (6), (7) 
приведена зависимость для расчета энтальпии образования 1-го компонента и удельной теплоемкости. 

В данной модели прямая задача (решение системы дифференциальных уравнений) является как жесткой, так 
и нежесткой. Метод Гира, используемый для решения модели каталитического риформинга, основан на 
дифференцировании назад [13-15]. Для сравнения результатов и выявления оптимального метода 
использовались явный одношаговый метод численного интегрирования Рунге-Кутты четвертого порядка 
и явный метод численного интегрирования зс1ру.о4еш() языка программирования Руоп. 
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Концепцию автоматизации формирования математического описания кинетики для многостадийных 
химических реакций и численного метода решений проиллюстрируем с помощью методологии ШЕЕО. 

На рис. | представлена контекстная диаграмма процесса автоматизации формирования математического 
описания кинетики для многостадийных химических реакций и решения численным методом. 


Метод формирования 
математического 
описания 


Метод Рунге-Кутты 
4-го порядка 


Метод 
зсру.одешк) Метод Гира 


Стехиометрическая 
а Автоматизация формирования 
пр математического описания кинетики Численное решение 
к я 
для многостадийных химических прямой задачи 
Экспериментальные реакций и численный метод решения 

данные > 

НИ 5 


Разработанный 
программный 
комплекс 


Рис. 1. Контекстная диаграмма «Автоматизация формирования математического описания кинетики для 
многостадийных химических реакций и численный метод решения» 


Входные параметры — стехиометрическая матрица реакции и экспериментальные данные. В их число 
входят начальные данные по концентрации веществ, температуре, предэкспоненциальным множителям, 
энергиям активации. Результатом решения задачи является численное решение прямой задачи. В качестве 
элементов управления выступают численные методы для сравнения результатов, а также метод формирования 
математического описания (СОНДУ). В качестве механизма выступает разработанный программный комплекс 
автоматизации математического описания кинетики многостадийных реакций и решения прямой задачи. На 
рис. 2 приведена декомпозиция процесса автоматизации формирования математического описания кинетики 
для многостадийных химических реакций и решения численным методом на этапы: 

— формирование математического описания; 

— решение прямой задачи. 


Метод Рунге-Кутты 
Метод формирования Метод Гира 4-го порядка 
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Рис. 2. Этапы решения задачи «Автоматизация формирования математического 
описания кинетики для многостадийных химических реакций и численный метод решения» 
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Возможности разработанного программного комплекса представлены описанными ниже особенностями. 

При помощи программы можно создать кинетическую модель многостадийной химической реакции 
и представить её решение методами, описанными ранее, используя заданную стехиометрическую матрицу. 
Если реакция неизотермическая и происходит изменение объема реакционной смеси, к системе 
дифференциальных уравнений добавляются два дополнительных уравнения для описания этих изменений. 

Математическая модель каталитического риформинга бензина представляет собой систему 
дифференциальных уравнений (СДУ) вида (1)-(7), состоящую из 40 дифференциальных уравнений по числу 
групповых компонент, температурной зависимости и уравнения изменения молей. 

Основной функционал программы представляет собой: 

— главное окно приложения с выводом математического описания по стехиометрической матрице, 
введенной пользователем; 

— окно с визуализацией решения: для каждого численного метода (Руфоп-одешть, Рунге-Кутта 4-го порядка, 
метод Гира) выводится время работы, массив времени { и решение системы дифференциальных уравнений 
в каждый момент времени, которое представляет собой графики изменения концентрации веществ, изменения 
температуры реакции или объёма реакционной смеси. 

Результаты исследования. Реакция каталитического риформинга протекает в трех реакторах, между 
которыми располагается нагревательная печь. Начальная температура первого реактора равна 766 К (493 °С). 
Первые 9,6 условных часа реакция протекает в первом реакторе, температура падает. После этого происходит 
нагревание в печи между двух реакторов до 763 К (490 °С) и химический процесс переходит во второй реактор. 
Во втором реакторе процесс протекает в интервале времени от 9,6 до 32,3 условных часа, температура 
понижается. Прежде чем перейти к 3 реактору, смесь нагревается во второй промежуточной печи до 768 К 
(495 °С) и протекает в этом реакторе в интервале времени от 32,3 до 60 условных часа. 

На рис. 3—5 изображены кинетические кривые групп углеводородов, участвующие в реакции и полученные с 
помощью метода Гира. Скачкообразные переходы на графиках свидетельствуют о том, что смесь переходит из 
одного реактора в другой с промежуточным нагревом. 
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Рис. 3. Кинетические кривые каталитического риформинга бензина 
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Кривые изменения концентраций групповых компонент приведены в заисимости от времени контакта 
реакционной смеси с катализатором (рис. 3—5). Например, символ пР1 обозначает нормальный парафин с одним 
атомом углерода, а символ АЭ — ароматический углеводород с девятью атомами углерода. На рис. б и 7 
представлены соответственно кривые изменения температуры и объема реакционной смеси реакции 


часы 


Рис. 5. Кинетические кривые каталитического риформинга бензина 


каталитического риформинга бензина. 
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Рис. 6. Кривая изменения температуры 
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Рис. 7. Кривая изменения объема реакционной смеси 


Обсуждение и заключение. В представленной работе разработан программный модуль формирования 
математического описания кинетики для сложных многостадийных реакций. Это позволяет анализировать сложные, 
в том числе промышленные процессы, рассматривать разные возможные механизмы за приемлемое время. 

Разработанный программный модуль численного решения прямой кинетической задачи позволяет 
отслеживать динамическое изменение концентраций компонент и температуры процесса в случае его 
неизотермичности. Помимо представленных графиков изменения концентраций компонент, произведены 
замеры времени работы каждого из трех реализованных методов. По методу зсру.одеш() решение прямой 
задачи длилось 80,5 секунды, по методу Рунге-Кутты 4-го порядка — 111,8 секунд. Наиболее оптимальное 
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время, равное 39,8 секунды, показал метод Гира. Это подтверждает, что метод Гира оптимизирован для работы 
с жесткими системами дифференциальных уравнений. 

Реализованный подход применен для моделирования промышленного процесса каталитического 
риформинга бензина. Автоматизация формирования математического описания кинетики и реализация 
численного решения прямой задачи позволили в приемлемые сроки разработать кинетическую модель 
детализированной схемы химических превращений и получить концентрационные и температурные профили 
в реакторах процесса, соответствующие промышленным данным. 


Список литературы 

1. Нарышкин Д.Г. Компьютерная математика в курсе "Физическая химия". Современные информационные 
технологии и ИТ-образование. 2012;(8):430-440. ОВГ: Б&рз://субе|йепшка.ги/агис]е/п/котпрушегпауа-тайетайКа- 
у-Кигзе-Я1свезКауа-Битуа/\1е\уег (дата обращения: 15.08.2023). 

2. О! Нап, Хлап-Та! ХВой, Х1ао-О! Не, Нопг-Вше Л. Меспап15т ап Ктейс$ ог фе Аегоыс Ох1даНоп оЁ 
Веп7у! А1сово| ю Вепха!ерву4е Саутеа Ъу СоБаЁ Рогрцупп ш а Метбгапе Мисгосваппе!| Кеасюг. Сйеписа 
Епетеетте степсе. 2021;245:116847. В рз://401.0ге/10.1016/.сез.2021.116847 

3. ОКеогове Мапа. А ССМ-Вазеа Мод\аг апа Нубла Ктейс Моде! ю Зипиае Фе Тгурюрвап 5упез1$ ша 
Ееа-Ваюп В!огеасюг Озше Моеа Е. Сой Се$. Сотршету апа Сйетса! Епетеетия. 2021;153:107450. 
6 рз://401.оге/10.1016/.сотрсветепз.2021.107450 

4. Файзуллин М.Р., Балаев А.В. Автоматизированная система исследования кинетики сложных химических 


реакций. Вестник Башкирского университета. 2008;13(3—1):835-839. 

5. Хамидуллина 3.А., Исмагилова А.С., Спивак С.И. Анализ информативности кинетических параметров 
сложных химических реакций. Вестник Тверского государственного университета. Серия: химия. 
2020;39(1):70-80. В рз://401.ог®/10.26456/у{спет2020.1.9 

6. Ко!едта К.Е., СибаудяШа 1.М., Коедт 5.М№. Мафетайса! Модейп> ап Сотрщайопа! Арес оЁ Мий- 
Сгиепа ОрНпитаНоп оЁР ше СопаНюп$ оЁ а Габогаюгу СайаЙуйс КеасНоп. М№Митейса![ Апа/уя5 ап Аррйсанопз. 
2022;15:104—111. В ф://Ч01.оге/10.1134/$1995423922020033 

7. ХашияШи В.7., Коедта К.Е., АКБтеюу А.Е., СибаиШт Т.М. Ктейс$ оЁ Фе Саайуйс КеюЮгишя оЁ 
СазоПпе. Атейсу апа Сао!у5. 2017;58(3):279-289. ВИрз://401.0г2/10.1134/$0023158417030132 

8. Зайнуллин Р.З., Коледина К.Ф., Ахметов А.Ф., Губайдуллин И.М. Возможные пути модернизации 
реакторного блока каталитического риформинга на основе кинетической модели. Еестотс баепийс Лоигпа! ог 
ОИ апа Са5 Визтезз. 2018;(6):78-97. 

9. Вгаззеиг @.Р., Гасоб О.Т. Модейия о} Аттозрйетс Сйепияту. Сатби4ее: Сатбиаее Ошуегзиу Ргезз; 2017. 
606 р. Би рз://401.ог2/10.1017/9781316544754 

10. СибаудаШи Т., Коеата К., бауШпа Г. Мафетайса! МодеНпе оЁ ш4асйоп Рено оЁ Фе Оейп$ 
Нуагоавттайоп Веасйоп Фу ПизобиуааттилисЬюнае Саутеа узи Ср›7гСЬ. Епетеетие Лоигпай. 
2014;18(1):13—24. В&рз://401.0г2/10.4186/е].2014.18.1.13 

11. Мезба т У.Р., Роу! У.С., ЗоБо]еуа 1.У. Ме#о4 апа А!гогийт юг Вауезап Езитайоп оРКтейс Рагатеег$ 
ог Свеписа| Ргосеззез у Ешсшайп» Шш4ереп4еп{ Уапаез. Тйеогейса! Коипаапопх ор Српетса! Епатеетив. 
2009;43(6):861—868. В&р://401.ог2/10.1134/50040579509060037 

12. Сафиуллина Л.Ф., Коледина К.Ф., Губайдуллин И.М., Зайнуллин Р.3. Исследование математической 


модели каталитического риформинга бензина методами анализа чувствительности. Вычислительные методы и 
программирование. 2020;21(4):440—451. В рз://401.0ге/10.26089/МитМе(.у2 11435 

13. Мигрву Т.А. РВуяса| Огвапюе Свепизву. Вейяет Лоигпай ор Отгатшс Сйетияйу. 2010;(6):1025. 
66 р5://401.0г2/10.3762/5]0с.6.116 

14. ХаупяШо В.7., Коедта К.Е., СифаудиШпт 1.М., АКБтеюу А.Е., Ко]едт $.М. Ктейс Моде! ог Са уйс 
СазоПпе Веюгите м Сопзегайоп Рог СБапзез п фе Веасйоп Уоите апа ТВегло4дупапс Рагатецегз. Атейсу 
апа Сша[руя5. 2020;61:613—622. ВИрз://401.0г2/10.1134/$002315842004014Х 

15. ЗайяШпа Г..Е., СибаудиШт Т.М., Отуапбаеу В.М., Мизша А.Е. Сотрщайопа! Азрес{ оЁ ЗпарИЙЯсаНоп оЁ 
Мафетайса| Мо4е!5 оГ СБеписа! Веасноп Зузет1з. Лоигпа/ ор Рйу5с5: СотГегепсе 5етеу. 2019;1368:042022. 
Вр://401.0г2/10.1088/1742-6596/1368/4/042022 


Информатика, вычислительная техника и управление 


407 


В р://уезиК-опза.га 


408 


'А4уапсей Епотеениз Везеагсй (Коз1оу-оп-Ооп). 2023;23(4):398—409. е155 № 2687—1653 


ВеЕегепсе5 

1. МагузКш ОС. Сотршег Мафетайс$ ш Фе Сочгзе “РБуз1са! СБепизку”. Модети т/огтайоп Тесйпоовлез 
апа [Т-Едисаноип. 2012;(8):430—440. (ш Влз$.) 

2. ОГ Нап, Хлап-Та1 Ивои, Хлао-О! Не, Нопз-Вте Л. Месвап1зт апа Клмейс$ о Фе Аегоыс Ох14айоп оЁ Веп7у1 
А[соВо] +0 Вепха!4екуде Саёа]уте4 Бу Софа Рогрвупп ш а Метбгапе Мпсгосваппе! Кеасог. Сйеписа! Епетеетия 
Эстепсе. 2021;245:116847. В рз://401.0го/10.1016/.сез.2021.116847 

3. ОБеогеве Мапа. А ССМ-Вазеа Модшаг апа Нубиа Ктейс Моде! ю Зииае Фе Тгурюрвап Зуп®е$1$ ша 
Еед-ВажЬ В!огеасюг Озше МодШеа Е. Сой СеШз. Сотршету ап Сйетса[ Епетеетия. 2021;153:107450. 
6рз://401.оге/10.1016/.сотрсветепз.2021.107450 

4. ЕаухаШп МК, Вааеу АУ. Сотшрщег-А14е4 Зумет Юг Кшейс За41ез оГ Сошр! ех СБеписа| Веасйопз. 
Вийейп о ВБИ. 2008;13(3-1):835-839. 

5. КВапи4иШпа 7.А, [зтаПоуа АЗ, ЗрууаК 51. Апа1уз15 шоппайуйу оЁ Кшейс Рагатеегз оЁ Сотр!ех СКеписа1 
КеасНоп$. Нега о{ТубО. бетех “Српетяу”. 2020;39(1):70-80. Врз://401.0ге/10.26456/у1сВет2020.1.9 

6. Ко]едта КЕ, Сибауди п М, Коеаш ЗМ. Мафетайса! Мо4дейпе апа СотршаНопа| Азресё ог Мш@- 
Стцепа Орипитайоп оЁ Фе Сопа!юоп$ оЁ а Габогаюгу Сауйс Веасйоп. Митейса! Апа/у515 апа АррИсапопз. 
2022;15:104-111. В р://Ч01.оге/10.1134/$1995423922020033 

7. ГашаяШи В7, Коедта КЕ, АКЬтеюу АЕ, СибаиШт 1М. Ктейс$ оЁ ше Саауйс Кеюгиите ог СазоПпе. 
Ктейсх апа Сшайряю. 2017;58(3):279—289. В@рз://401.0тг2/10.1134/$0023158417030132 

8. ГашиШи В7, Коедта КЕ, АКбтеюу АЕ, СибааШт 1М. Розе У/ауз оЁ Оретадте Веасюг Опй оЁ 
Санс КеЮгитз Вазе4 оп Кшейс Моде!. Еестотс 5слепИйс Лоигпа[ оГ ОЙ апа Саз Визтезз. 2018;(6):78-—97. 

9. Вгаззеиг СР, Лласоб О. Мо4ейпе ор Автозрйетс Сйепияту. Сатбиаее: СатбиАее Ошуегзиу Ргез$; 2017. 
606 р. 6 рз://401.ог2/10.1017/9781316544754 

10. Сибау4аиШо 1, Коедта К, бауШпа Г. Мафетайса! Модейпе оЁР шаисНоп Репо@ оЁ Фе Оейп$ 
Нуагоаттайоп Веасйоп Бу ПизобиУааттилисЬюнае Саевутеа узи Ср›7гСЬ. Епетеетие Лоитпай. 
2014;18(1):13—24. №&рз://401.0г2/10.4186/е].2014.18.1.13 

11. Мезба т УР, Роз! УС, боБоеуа ГУ. Ме#о4 апа А1вогибт г Вауезап Езитайоп оЁ Ктейс Рагатаег оЁ 
Среплса! Ргосеззез у Ешсаайпе ш4ерепдет Уапа ез. Тйеогейсай Гоип4айоп$ ор Сйеписа! Епатееттв. 
2009;43(6):861-868. №@р://401.оге/10.1134/$0040579509060037 

12. ЗайаШпа ГЕ, Ко]едта КЕ, СаБбаудиШп ТМ, ХаупаШп ВЯ. Зва4у оЁ а Мафетайса! Моде оЁ СазоПпе 
Сауйс Веюгите Бу Зепяйуйу Апа!у515 Мефод$. Митенса! Мето45 апа Ргоетатпиив. 2020;21(4):440—451. 
№ рз://401.0г2/10.26089/МитМе.у2 11435 (ш Ви$$.) 

13. Мирву ЛА. РБузса! Огватюе СБепизйу. ВеЙяет Лоитпай оГ Отгзатс Сйепизйу. 2010;6):1025. 
6 р5://401.0т2/10.3762/5]0с.6.116 

14. ГаупаШа ВИ, Коедта КЕ, Сибау4и т ТМ, АКАтеюу АР, Ко!едт ЗМ. Ктенс Моде! оё Сауйс СазоПпе 
Веюгиш» у Сопз1еганоп г СКапзез ш фе Кеасйоп Уоите апа Твептодупапис Рагатеегз. Клиейсх апа 
Сашуяю. 2020;61:613—622. В рз://401.ог2/10.1134/$002315842004014Х 

15. ЗайаШпа ГЕ, Сибауа Шт М, Отуапфаеу КМ, Мизша АЕ. Сотшршайопа| Азресё5 оЁ ЗппрИЙЯсайоп оЁ 
МафетаНса| Мо4де!$ оЁ СБеписа| КеасНоп Зузетз. Лоигпа[ оГ Рйуяс5: Сопфегепсе 5етеу. 2019;1368:042022. 
Бр://Чо1.оте/10.1088/1742-6596/1368/4/042022 


Поступила в редакцию 07.10.2023 
Поступила после рецензирования 02.11.2023 
Принята к публикации 16.11.2023 


Об авторах: 

Никита Андреевич Лысенко, магистрант кафедры информационные технологии и прикладная математика 
Уфимского государственного нефтяного технического университета, (450064, РФ, г. Уфа, ул. Космонавтов, 1), 
ЭРГ\-код: 5938-7991, ОВСТО, Везеагсвег О, икИКа_[узепко_2016@тай.га 


Камила Феликсовна Коледина, — доктор физико-математических наук, профессор — кафедры 
информационные технологии и прикладная математика Уфимского государственного нефтяного технического 
университета, (450064, РФ, г.Уфа, ул. Космонавтов, 1), старший научный сотрудник лаборатории 
математической химии Института нефтехимии и катализа УФИЦ РАН (450075, РФ, г. Уфа, 
проспект Октября, 141), ЭРП\-код: 2449-0255, ОВСТО, ВезеагсНег О, Зсориз ЮО, Ко[едтаката@та!.ги 


Лысенко Н.А. и др. Автоматизация формирования математического описания кинетики для многостадийных химических реакций 


Заявленный вклад авторов: 

Н.А. Лысенко — разработка программного обеспечения, подготовка текста, формулировка выводов. 

К.Ф. Коледина — научное руководство, предоставление исходных данных, корректировка выводов, 
доработка текста статьи. 


Конфликт интересов: авторы заявляют об отсутствии конфликта интересов. 


Все авторы прочитали и одобрили окончательный вариант рукописи. 


Весеуеа 07.10.2023 
Веубед 02.11.2023 
Ассеред 15.11.2023 


'АБои! ше Аифогу: 

МКИКа А. ГузепКо, ип4егогадиае збва4епе оЁ йе шЮгтайоп ТесБпоюгу апа АррНеа Мафетайс$ Оерагитетв, 
О За Ремо]еит Тесппоое1са! Ошуетзйу (1, Козтопамюу 51., ЧЁ, 450064, ВЕ), ЗРПУ-соде: 5938-7991, 
ОБВСТО, Везеагспег О, иКИКа_[узепко_2016@та!.га 


Капшйа Е. Коедта, Ог.5с1. (РБуз.-Ма.), Ргоеззог оЁ пе шоппайоп ТесБпо]огу ап4 АррПе@ Маетайс$ 
Рерагалепь, Ча Зе Рего!еит ТесБпоо?1са! ОшуегзИу (1, Козтопауюу 56., О, 450064, ВЕ), Зешог ВезеагсВег, 
Габогафогу ог Мафетайса! СБетилзу, шзНвие ог Редосвепизту апа Са{а1у515, ОЕВС КА$З (141, рг. ОКбуабтгуа, О, 
450075, ВЕ), ЗРПУ-соде: 2449-0255, ОВСТО, Везеагсвег О, Зсориз ЮО, Кое4такатПа(@тай.ги 


С ратей сополшотяр: 
МА ГузепКо: зоЙ\уаге 4еуе!ортепь, {ехё ргерагайоп, огиайоп оЁ сопса10п5. 
КЕ Коеаша: асадепис а4у1$ те, ргоу15юп оЁР и а1 даа, соггесНоп оЁ сопсаз1оп$, геу1$1юп оЁ Пе 1ехё оРе агас[е. 


Соп/ПсЕ о 1теге5! зкиетеп!: 1е алфогз 40 по Вауе апу сопЯ1сЕ оЁ пуеге$4. 


АП ашфогу йауе геа4 ап арргоуе4 ше ппа тапизстри. 


Информатика, вычислительная техника и управление 


409 


Бйр://уезииК-допза.ги 


410 


А4уапсей Епэтеетпия Кезеагсй (Ко51оу-оп-Ооп). 2023;23(4):410—421. е155М№ 2687—1653 


ИНФОРМАТИКА, ВЫЧИСЛИТЕЛЬНАЯ 
ТЕХНИКА И УПРАВЛЕНИЕ 
МГОВМАТТОМ ТЕСНМОГОСУ, СОМРОТЕВ 
ЭСТЕМСЕ АХО МАМАСЕМЕМТ 


м Спеск тог ирда{е$ 
ВУ 


УДК 004.58 Научная статья 
6 рз://401.оге/10.23947/2687-1653-2023-23-4-410-421 


Оптимальная 20 -расстановка виртуальных объектов в физическом 
пространстве для приложений дополненной реальности 


М.В. Алпатова® ‚ Ю.В. Рудяк®@ ЕОМ: ОРЕЕМС 


Московский политехнический университет, г. Москва, Российская Федерация 


4 т.у.араюуа@уапдех.га 


Аннотация 


Введение. Научные и прикладные работы о размещении виртуальных объектов в реальном пространстве чаще 
всего фокусируются на вопросах интерактивности, интеграции реальности и виртуальности, физических 
свойствах виртуальных элементов. Однако недостаточно проработана задача одновременно свободного 
и оптимального размещения объектов с учетом их размеров и окружающей зоны комфортности вокруг них. 
В литературе можно найти описание схожей задачи — об упаковке в прямоугольный контейнер. В нашем 
случае цель не ограничивается максимально плотным размещением. Следует учесть два условия: жесткие 
размеры объектов (их запрещено нарушать) и дополнительные области — зоны комфортности (их 
нежелательно занимать). Цель работы — создание и реализация такого 2)-алгоритма размещения объектов 
в физическом пространстве, который будет учитывать обозначенные выше ограничения. 

Материалы и методы. Используя аппарат численных методов, авторы задействовали созданный ранее 
1-алгоритм размещения объектов. Расчеты основываются на системе линейных уравнений. В одномерном 
случае оптимальное размещение виртуальных объектов сводится к задаче, не зависящей от вида функции 
комфортности. Элементы такой системы — размеры объектов, дистанции между ними, а также расстояния 
до края области встраивания, зоны комфортности. Предлагаемый 2)-алгоритм оптимальной расстановки 
виртуальных объектов реализовали в виде программного кода на языке СЯ с использованием известного 
игрового движка Иийу. Решение тестировали на гаджетах в режиме пиковой нагрузки для 5, 10, 15, 20, 25, 35, 
40, 45 и 50 объектов. Для опытов задействовали 1,8 тыс. устройств. Проанализировали около 77 тыс. событий. 
Чтобы исключить нерепрезентативные значения, каждый расчет повторяли 10 раз, и для каждого значения 
провели 2-оценку. Аномальные (больше 3 и меньше -—3) исключили. 

Результаты исследования. В работе создан алгоритм 2О)-расстановки, который реализует заполнение 
прямоугольной области виртуальными объектами. У каждого из них есть размер и еще одна характеристика — 
зона комфортности. Авторы составили блок-схему реализации данного алгоритма в заданной двумерной 
левосторонней системе координат. Показано, в частности, на каком этапе объекты сортируются по длине, когда 
формируются их партии и выполняются расстановки по двум осям. Первая — горизонтальная, вторая 
направлена вперед от пользователя (это вектор глубины, или фронтальное измерение). Алгоритм 
1-размещения для сформированного ряда позволяет оптимально расположить объекты вдоль оси Х на основе 
рассчитанного коэффициента комфортности К. Выполнены расчеты и составлены схемы с целью достичь 
определенных показателей комфортности. Для каждого объекта первой линии смещение по оси 7 от края 
плоскости определяется так, чтобы комфортность спереди равнялась комфортности по Х. Начиная со 2-го ряда 
для вычисления отступа проверяется наличие потенциальных соседей, которые находятся на ряд выше и имеют 
общие участки по Х` с обрабатываемым объектом. Каждый элемент строки устанавливается по оси 7 так, чтобы 
его комфортность сверху была максимальной из односторонних горизонтальных комфортностей в данной 
и предыдущей строках. Принцип расчета координаты 7 для объекта строки представлен в виде блок-схемы. 
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Исходными данными для реализации этого алгоритма были 7 объектов с 14 разными размерами и 28 зонами 
комфортности. После программной реализации работу описанного 2)-алгоритма проверили на практике — 
в мобильном приложении дополненной реальности. Записали аналитические данные пользовательских сессий. 
Рассчитали среднее время выполнения. Возникшую в ходе работы гипотезу о квадратичной зависимости 
проверили на персональном компьютере. С этой целью провели аналогичный эксперимент для диапазона [10— 
10000] объектов. Гипотеза подтвердилась. Алгоритму можно присвоить сложность О(и?). Для сравнения 
скорости вычисления задействовали 10 самых популярных моделей пользовательских устройств. Результаты 
представили в виде диаграммы. Минимальное зарегистрированное время выполнения — 0,093 мс, 
максимальное — 0,146 мс. Расчеты показали высокую эффективность двумерного алгоритма. Дополнительно 
визуализировали схемы расстановки для разного количества и параметров объектов. 

Обсуждение и заключение. Предлагаемый алгоритм двумерного размещения позволяет работать с набором 
виртуальных объектов с разными размерами и зонами комфортности. Показаны достаточно высокие 
производительность и стабильность. В среднем алгоритм реализуется за доли миллисекунды даже при больших 
партиях объектов. Возможные будущие направления работы: 

— расширение подхода для построения ЗО моделей и алгоритмов; 

— включение в алгоритм вращения объектов для большей гибкости их расположения и лучшего использования 
пространства. 

Итоги работы могут представлять интерес для инженеров и дизайнеров интерфейсов. В перспективе следует 
изучить пользовательский опыт и возможности включения дополнительных ограничений на позиционирование. 


Ключевые слова: виртуальные объекты в физическом пространстве, виртуальные объекты в дополненной 
реальности, комфортное размещение виртуальных объектов 
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Введение. Известная задача приложений дополненной реальности (апотеще4 геащу, АВ) — размещение 
виртуальных объектов в реальном физическом пространстве. Ряд исследований [1-3] фокусируются на 
интерактивности, слиянии реальной и виртуальной сфер. Изучаются также вопросы физических свойств 
виртуальных объектов. 

Недостаточно проработана задача такого размещения объектов, которое одновременно было бы свободным 
и оптимальным, учитывало не только геометрические размеры объектов, но и зоны комфортности вокруг них. 
В [6] отмечена внешняя схожесть этой задачи с известной задачей об упаковке в прямоугольный контейнер [5]. 
Однако есть существенная разница. В исследуемом случае недостаточно обеспечить максимально плотную 
упаковку. При размещении виртуальных объектов нужно учесть не только их жесткие размеры, которые нельзя 
нарушать, но и дополнительные области — зоны комфортности. Их занимать нежелательно. Эти 
дополнительные области позволяют приблизиться к объекту и выполнить с ним какие-то действия. 
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Целями представленной работы являются создание, реализация и оценка 2)-алгоритма оптимального 
размещения виртуальных объектов в физическом пространстве. 

Материалы и методы. Ранее [4] авторы определили, что такое комфортное размещение виртуальных 
объектов, ввели понятие функции комфортности А(х). Она монотонно возрастает от 0 до 1 при 0<х<1 и равна 1 
при х>1, где х=АХ/О, Х — расстояние от края объекта до ближайшего препятствия, р — размер зоны 
комфортности. Для каждого измерения у объекта есть две односторонние зоны комфортности О- и Д-. 
Комфортность с каждой стороны вычисляется отдельно. 

В [6] показано, что задача одномерного размещения п виртуальных объектов в свободной области 
пространства протяженностью Г сводится к системе линейных уравнений, не зависящих от вида функции 
комфортности К(х): 


хо хх) 

ро ро’ 

х0 хе 

5 -`Во ( 

50 С 

„ 2-У(х0-+9 

хе _ Ум) 

50) 50) 
Здесь Х® — расстояние первого объекта от левого края области встраивания; Х®, {=2,3,...п — расстояние 
между объектами с номерами Ги (1—1); Б® и ро — левая и правая зоны комфортности соответственно; /® — 
размер объекта; Б® = 9 =р®+ р, 1=1,2,..., (п-1). 


Матрица системы (1) сильно разрежена, поэтому можно избежать использования не самых быстрых 
универсальных методов и достаточно просто найти решение. Например, в первых (и - 1) равенствах можно 


в каждом Г-м уравнении выразить Х““ через Х®, затем подставить это в последнее уравнение и получить 
> 1 
линейное уравнение относительно Х“. После этого от первого уравнения к (и-1)-му последовательно 


2 Э КА 
определяются значения хх... хо. Авторы реализовали данный 1Л-алгоритм размещения объектов 


с помощью численных методов. Он показал высокую скорость и эффективность. Этот 1Ш-алгоритм стал 
основой для научных изысканий, описанных в представленной работе. Авторы предложили 2)-алгоритм 
оптимальной расстановки виртуальных объектов. Его реализовали в виде программного кода на языке С# 
с использованием игрового движка Иийу, который широко используется для создания мобильных приложений 
дополненной реальности [8]. В ходе эксперимента на устройстве пользователя, работающего в режиме пиковой 
нагрузки, запускался алгоритм для 5, 10, 15, 20, 25, 35, 40, 45 и 50 объектов. Замеряли время, затраченное на 
вычисления (в миллисекундах). Каждый расчет повторяли 10 раз во избежание аномальных значений. Всего 
в эксперименте участвовали 1,8 тыс. устройств, с них собрали и проанализировали около 77 тыс. событий 
с результатами. 

Фильтрация аномалий с помощью 2-оценки позволила выявить значения, которые можно определить как 
выбросы [9]. Авторы стандартизировали данные и вычислили 2-оценку для каждого значения. Аномальными 
сочли те, у которых 2-оценка оказалась больше3 или меньше-3. Их исключили и сосредоточились 
на типичных и репрезентативных данных. 


Результаты исследования. Итак, на плоскости Х, 2 есть прямоугольная область с шириной Г, и длиной Г.. 
Ее следует заполнить некоторым количеством виртуальных объектов. Каждый объект помимо своих размеров 
19 и 19 (где { — номер объекта) характеризуется еще зонами комфортности 20,09, р0,ро . На рис. 1 


представлена блок-схема описываемого алгоритма. 
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Установка * Размеры области 


параметров * Список объектов 


Сортировка 
объектов по 
длине 7 


Уведомление о 
завершении 
расчетов 


Формирование 
партии 
(Баёсь) 


Бас != поП Нет 


Да 


у 


Расстановка по Х 


Расстановка по У 


Рис. 1. Верхнеуровневый алгоритм 2)-расстановки 


В данной работе используется левосторонняя система координат: ось Х обозначает горизонтальный вектор 
размещения, ось Х направлена вперед от пользователя и является вектором глубины или фронтальным 
измерением [7]. Объекты размещаются в горизонтальной плоскости (например, на полу, столешнице и пр.), 
поэтому вертикальная ось У не рассматривается, а все иллюстрации предполагают вид сверху на получившуюся 
композицию. 

В левом нижнем углу каждого виртуального объекта обозначим опорную точку (р1уо. Она будет 
определять конечные координаты размещения объекта. Далее набор объектов сортируется в порядке убывания 


по 19 +50 + ро — габаритной длине объектов с учетом зон комфортности, расположенных спереди и сзади 


в соответствии с рассматриваемой осью координат. После сортировки объекты, занимающие больше 
пространства по глубине, располагаются дальше от пользователя (или выше для схемы, представленной 
в горизонтальной проекции). Блок-схема формирования партии из объектов представлена на рис. 2. 


Загрузка нерасставленных 
объектов, пока сумма их 
широт не превысит размер 
плоскости 


Направление текущей 


Нет Инвертирование порядка 
строки слева направо? 


объектов в строке 


Смена флага 

Инкремент индекса челночного хода Создание объекта 
текущей строки для следующей новой строки 

строки 


Расчет Ш 
оптимального 


положения 


Рис. 2. Алгоритм формирования ряда из объектов 


Затем из упорядоченного множества объектов формируются отдельные строки или партии (Баёс|). В каждую 
строку включается минимальное количество объектов, сумма ширины которых вместе с зонами комфортности 
превосходит ширину заполняемой области Г. Если при добавлении следующего объекта вычисленная 


занимаемая ширина превысит горизонтальный размер доступного пространства внутри прямоугольника, партия 
считается завершенной. Таким образом, каждая партия может быть размещена в виде линии в пределах 
области. Ряды чередуются слева направо и справа налево для смешивания крупных и мелких объектов. Это 
можно назвать «челночным ходом» (рис. 3). 
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[ния 


2 


Рис. 3. Схема «челночного хода» на координатной плоскости 


Для сформированного ряда применяется алгоритм 1Ш-размещения, который позволяет оптимально 
расположить ряд объектов ВДОЛЬ оси Х на основе рассчитанного коэффициента комфортности К, одинакового 
для всех объектов ряда. Для каждого объекта первой линии смещение по оси И от края плоскости определяется 
таким образом, чтобы комфортность спереди равнялась комфортности по Я: 

(1) () 70 (2) 

2, рол ОИ 

т Р.И, (2) 

ро ро ро Во 


2+ х- 2+ х- 


Начиная со 2-го ряда для вычисления отступа 7+ проверяется наличие потенциальных соседей сверху из 
вышележащего ряда, имеющих с текущим обрабатываемым объектом общие участки по координате Х. При 
этом каждый объект размещаемой строки устанавливается по оси Й таким образом, чтобы его комфортность 
сверху была максимальной из односторонних горизонтальных комфортностей в данной и предыдущей строках. 
Алгоритм расчета координат 7 представлен на рис. 4. 
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Рис. 4. Принцип расчета координаты 7 для объекта строки 


Из блок-схемы видно, что отступ вычисляется аналогично для всех трех возможных сценариев: 

— отступ от края плоскости для установки объекта; 

— отступ от соседа сверху; 

— при отсутствии соседа в предыдущей строке берется самый близкий по 7 объект из всех предыдущих строк. 
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Новые партии формируются до тех пор, пока всем объектам не присвоят ряд или пока не закончится 
физическое пространство. 

Исходные данные к описанному алгоритму систематизированы в таблице 1. Объекты рассортированы по их 
вертикальным размерам. 


Таблица 1 
Исходные данные к примеру работы описанного алгоритма 
Номер объекта Размер Зона комфортности 
1 (7; 4) (4; 3; 3; 5) 
2 (10; 4) (4; 3; 5;3) 
3 (5; 4) (2; 3; 5,25; 2,25) 
4 (12;4) (2,8; 3,3; 3,8; 4,65) 
ь (5; 4) (2,5; 4; 1,75; 5,8) 
6 (18; 4) (3,2; 2,4; 5; 2,25) 
7 (11;4) (2,1; 2,45; 2; 14,2) 


Таблица 1 соответствует рис. 5, демонстрирующему размещение объектов. 


Рис. 5. Схема расположения объектов для рассматриваемого примера. Вид сверху 


Серый прямоугольник — это свободная площадь, на которой расставляются объекты. Черные области — 
объекты в их габаритных размерах. Цветные, с наложениями, — предустановленные зоны комфортности для 
каждого объекта. Как отмечалось выше, зоны комфортности объектов могут частично перекрываться 
и выходить за пределы свободного пространства, но применяемый подход позволяет обеспечить баланс, чтобы 
в равной мере уменьшались комфортности различных объектов. 

Напомним, что в данную модель входит понятие функции комфортности /(х). От нее не зависит полученная 
система уравнений (1), а значит и оптимальное размещение объектов. А вот значения односторонних 
комфортностей объектов определяются и их размещением, и видом функции А(х). Приведем два примера: 


— в случае линейной функции А(х) односторонние комфортности объектов верхнего ряда рис. 5 оказываются 
равными 0,6; 


Бйр://уезииК-допза.ги 


2 
— для зависимости д/]- (х _ 1) — значение 0,9. 


После программной реализации описанного 2)-алгоритма его внедрили в мобильное приложение 
с дополненной реальностью и записали аналитические данные пользовательских сессий. Сгруппировали 
информацию по каждой уникальной модели устройства, по каждому количеству расставляемых объектов 
в диапазоне [10—50]. Рассчитали среднее время выполнения (рис. 6). 
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Рис. 6. Время выполнения алгоритма на пользовательских устройствах: ® — реальные данные; а = 0,00498, 


р = 0,10593 — линейная аппроксимация; а = 0,00003, = 0,00317, с = 0,12743 — квадратичная аппроксимация 


На указанном диапазоне невозможно сразу определить точную сложность вычисления алгоритма. Как видно 
из рис. 6, данные близки к линейной и квадратичной аппроксимации. С помощью метода остаточной суммы 
квадратов рассчитана количественная оценка обеих аппроксимаций и получены значения 0,00035 и 0,00021 
соответственно [10]. 

Для подтверждения гипотезы о квадратичной зависимости на персональном компьютере запустили 
аналогичный эксперимент для диапазона [10-10000] объектов. Рис. 7 подтверждает гипотезу, значит алгоритму 
можно присвоить сложность О(п?). Впрочем, ситуация с 10000 объектами, скорее, теоретическая, и на практике 
в дополненной реальности пользователь вряд ли будет работать более чем с 1-2 десятками объектов. 
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Рис. 7. Время выполнения алгоритма для большого количества объектов: ® — реальные данные; а = 0,1365, Е 
Ь = 45,455 — линейная аппроксимация; а = 0,00001, Б = 0,0103, с =-—1,6178 — квадратичная аппроксимация Е 
=] 
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моделей и по результатам вычислений построили диаграмму (рис. 8). 
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Рис. 8. Сравнение скорости вычисления на разных устройствах для 10 объектов 


Минимальное зарегистрированное время выполнения — 0,093 мс, максимальное — 0,146 мс. Для данной 
выборки корреляция между моделью устройства и временем выполнения — всего 0,177. Это довольно низкое 
значение. 

Расчеты показали высокую эффективность двумерного алгоритма. Это обусловлено разбиением всей 
совокупности объектов на отдельные строки и использованием в каждой строке быстрого одномерного 
алгоритма, описанного в начале статьи. Этот алгоритм позволяет решить задачу с меньшими затратами 
ресурсов и времени. 

На рис. 9 представлены дополнительные визуализации схем расстановки для разного количества объектов и 
их параметров. 


а) 6) в) 
Рис. 9. Примеры расстановки разного количества объектов: 
а — 4 объекта; 6 — 9 объектов; в — 21 объект 


Итак, предлагаемый алгоритм двумерного размещения позволяет работать [6 набором виртуальных объектов. 
Каждый из них характеризуется определенными размерами и зонами комфортности. Алгоритм предназначен 
для оптимального расположения таких объектов в прямоугольном физическом пространстве, которое 
воссоздает окружение пользователя. 

Обсуждение и заключение. Для ранее описанной одномерной задачи предложен и реализован 
в программном коде эффективный алгоритм оптимального размещения виртуальных объектов. Для двумерного 
случая построена модель оптимального размещения. На ее основе разработан эффективный алгоритм, 
реализованный в программном коде. 

Проанализирована производительность алгоритма, которая измеряется временем выполнения. Установлено, 
что нет значимой корреляции между временем выполнения и моделью устройства. Среднее время — доли 
миллисекунды Даже Для больших партий объектов. Это говорит об относительно стабильной 
производительности алгоритма в различных условиях. Однако нельзя исключить возможное влияние других 
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непроверенных факторов. Для полной характеристики зависимостей производительности нужны дальнейшие 
исследования. 

Обозначим несколько направлений будущей работы для повышения производительности и адаптивности 
алгоритма. Это, во-первых, совершенствование решения, которое позволит использовать полные размеры 
объектов, по возможности без пересечений комфортных зон. Их размеры будут уменьшаться только в случае 
крайней необходимости. Это должно обеспечить оптимальное использование пространства и дальнейшее 
улучшение интерактивного опыта пользователя. 

Во-вторых, включение вращения элементов в алгоритм повысит гибкость в их расположении и даст 
возможность лучше использовать пространство. К тому же это позволит алгоритму учитывать больше разных 
объектов, тем самым расширяя его практическое применение. 

Заметим, что в приложении с дополненной реальностью расстановка пользователем трех и более 
объектов — достаточно трудоемкая задача. Она требует значительного времени и существенно затрудняет работу. 

Предложенная технология размещения объектов в 20 позволяет оптимально позиционировать виртуальный 
контент для улучшения пользовательского опыта в дополненной реальности. Автоматизированный расчет 
расположения сокращает ручную работу и дает возможность сосредоточиться на полноценном взаимодействии с АК. 

Научные изыскания в данном направлении ведут к бесшовной интеграции виртуального и физического 
пространств. Результаты могут применить на практике разработчики очков дополненной реальности, которые 
сталкиваются с проблемой быстрого интеллектуального позиционирования объектов в зависимости 
от пользовательских предпочтений и контекста. Кроме того, итоги работы могут представлять интерес для: 

— инженеров, работающих над умными домашними системами; 

— дизайнеров интерфейсов, для которых важно эффективно и гармонично сочетать элементы управления 
и информацию на экране. 

Введение нечеткой логики позволяет разработчикам адаптировать решение под конкретные задачи 
и сценарии использования. В будущих исследованиях целесообразно изучить и задействовать пользовательский 
опыт, включать дополнительные ограничения на позиционирование объектов. 
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У\Уе-приложение для математического моделирования 
нестационарного течения нефти в пористой среде 
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Институт нефтехимии и катализа Уфимского федерального центра Российской академии наук, 

г. Уфа, Российская Федерация 
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Аннотация 

Введение. Основу для анализа и интерпретации данных гидродинамических исследований составляют 
теоретические модели и вычислительные алгоритмы. Но, несмотря на востребованность данной тематики, 
многие вопросы, связанные с нестационарными течениями жидкостей в нефтяных пластах, все еще требуют 
решений. Поэтому математические постановки задач, связанных с учетом нестационарного течения жидкости, 
разработка эффективных численных методов и алгоритмов, их решение с применением современных 
уеБ-технологий являются актуальными. Цель данной работы — разработка \еб-приложения для 
математического моделирования процесса фильтрации жидкости в однопоровых коллекторах при проведении 
гидродинамического исследования на добывающей скважине. 

Материалы и методы. Для решения поставленной задачи применены методы механики сплошной среды 
и вычислительной математики. Представлена модель течения нефти в однопоровом коллекторе. При разработке 
приложения использованы языки программирования Руфоп и Лауабсире. Результаты расчетов хранятся 
в реляционной базе данных, реализованной средствами РозёотеЗ ОГ. 

Результаты исследования. Разработано новое \еБ-приложение для моделирования процесса фильтрации 
нефти в однопоровых коллекторах, которое применимо для изучения флюидодинамических процессов и может 
быть использовано для прогноза дебитов, добычи и расчета оптимальных режимов работы скважин. 
Обсуждение и заключение. Разработанное у\еб-приложение позволяет строить поля давления и температуры 
в пласте около работающей и остановленной добывающей скважины и на различных расстояниях от нее. 
Данная информация дает возможность оперативно оценивать длительность гидродинамических исследований, 
а также регулировать работу скважин. Приложение может быть развернуто в существующей сетевой 
инфраструктуре, пользоваться всем функционалом, подключаясь к удаленному серверу. Оно оптимизировано 
для использования на различных платформах и имеет широкие перспективы дальнейшего развития. 


Ключевые слова: у\еБ-приложение, математическая модель, гидродинамические исследования, база данных, 
давление, температура 
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Введение. Российская Федерация является одной из ведущих стран мира по объемам производства нефти. 
Каждый год на месторождениях страны добывается свыше 500 миллионов тонн нефти, подавляющая часть 
которой в сыром виде экспортируется в страны азиатско-тихоокеанского региона, Европы и ближнего 
зарубежья. Развитие нефтяной промышленности в нашей стране оказывает влияние на другие отрасли и сферы 
деятельности, поэтому решение задач по эффективному освоению и разработке нефтяных месторождений было 
и остается перспективным. Основой современных научных технологий нефтяной разработки является 
комплексное изучение свойств продуктивных коллекторов, содержащихся в них веществ и происходящих 
сложных процессов при течении флюида. 

Главная цель при добыче нефти — увеличение темпа отбора флюида из продуктивного пласта 
и поддержание пластового давления. Планирование работ по эффективному освоению и оптимальной 
эксплуатации нефтяной скважины носит комплексный характер, связанный с гидродинамическими процессами, 
происходящими при течении жидкости в системе «скважина — пласт» [1]. Для изучения продуктивных свойств 
пластов и жидкостей применяются гидродинамические исследования скважин (ГДИС). Средствами ГДИС 
проводятся замер совокупности данных коллектора с использованием глубинных или устьевых приборов, 
а также последующая обработка замеров и интерпретация полученных данных (продуктивность, 
фильтрационные свойства жидкостей, тип коллектора и т. д.) [2]. 

При исследовании коллекторов с использованием ГДИС применяются теоретические модели, основанные 
на классических уравнениях и законах, а также вычислительные алгоритмы для моделирования процессов 
фильтрации флюидов. Однако, несмотря на востребованность данной тематики, многие вопросы, связанные 
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с математическими и вычислительными аспектами моделирования нестационарных потоков в коллекторах, 
требуют дальнейшей разработки и развития [3, 4]. 

Для анализа и интерпретации результатов ГДИС существует ряд отечественных («Автон», «Гидрозонд») 
и зарубежных (ОГСА, ГедаРю\у, Е!о\/У191юп, МАЗТ, барыг МГ.) программ [5-7]. При всех имеющихся 
преимуществах они обладают некоторыми недостатками, такими как невозможность проведения полного 
спектра расчетов, просмотра и изменения параметров исходных моделей. Все представленные программные 
продукты являются настольными, коммерческими, сложными в эксплуатации, многие обладают 
перегруженным интерфейсом. В условиях, когда необходимо оперативно или в режиме реального времени 
провести расчет, большие коммерческие симуляторы оказываются малоэффективными. Поэтому для 
выполнения схожих работ предлагается использовать \еб-приложение, обладающее множеством преимуществ. 

В отличие от настольных, \еб-приложения независимы от операционной системы и вычислительных 
мощностей персональных компьютеров. Программный код \еб-приложения пишется единожды для 
конкретной платформы, где оно будет развернуто. \У’еБ-приложения выделяются простой реализацией, 
легкостью применения, высокой производительностью. Их основными характеристиками являются 
доступность, надежность, безопасность, масштабируемость, гибкость, кросс-платформенность. Для работы 
с \уеБ-приложением от персонального компьютера пользователя необходимо находиться в одной сети 
с сервером, на котором оно развернуто. Такой сетью может быть локальная сеть предприятия или глобальная 
сеть Интернет, если речь идет о веб-сайтах. По сравнению с настольными, у \еб-приложений можно выделить 
один недостаток — время отклика, зависящее от удаленности персонального компьютера от сервера. 

Логика работы \еб-приложения отличается от настольного, которое запускается и выполняет код 
на компьютере пользователя. Поэтому настольные приложения отличаются более расширенным и отзывчивым 
пользовательским интерфейсом и позволяют реализовывать более сложную бизнес-логику работы. Однако 
с развитием уеб-технологий это достоинство с каждым годом нивелируется, и У’еб-приложения все чаще 
бывают способны реализовывать сложные пользовательские сценарии. Таким образом, в части реализации 
производственных задач \еб-приложения становятся более конкурентными. 

В связи с этим математические постановки задач с учетом нестационарного течения жидкости во всех его 
элементах (пласт, стволовые трубы, проточные каналы), разработка эффективных численных методов 
и алгоритмов, последующие их решения с применением современных технологий \еБ-разработки являются 
актуальными. Использование полученных численных алгоритмов позволяет количественно описывать 
поведение давления и температуры в нефтяных скважинах. Алгоритм и \еБ-приложение предоставят 
пользователям, специалистам в области анализа и интерпретации гидродинамических исследований, удобный 
и легко адаптируемый к их практическим потребностям инструмент для расчета сложных процессов, 
происходящих при разработке пористых коллекторов. 

Материалы и методы. Рассматривается добывающая скважина, на которой проводится гидродинамическое 
исследование методом кривой восстановления давления. Скважина эксплуатирует терригенный коллектор, 
который относится к однопоровым. Пласт является неограниченным, воздействие соседних скважин 
отсутствует. Течение флюида в пласте — горизонтальное, перетоки между пропластками отсутствуют. 
Скважина работает в установившемся режиме с постоянным забойным давлением. На границе контура 
поддерживаются постоянные давление и температура. Для проведения исследования измерительный прибор 
опускается на глубину верхних дыр перфорации. Расчет динамики давления и температуры проводится 
в пространственно-одномерном случае. 

При численном моделировании важно учитывать все параметры пласта и флюида, чтобы достичь наиболее 
качественного описания процесса массопереноса нефти в однопоровом коллекторе. Рассмотрим 
математическую модель, описывающую процесс фильтрации нефти в однопоровом коллекторе. Основой для 
модели служат уравнения пьезопроводности и теплопроводности. Математическая модель в радиальных 
координатах представляет собой систему дифференциальных уравнений в частных производных: 
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0Т _1 [2 о[бРньее + содР 
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гдеф — пористость пласта; с, — сжимаемость породы (1/Па); Е — проницаемость пласта (м?); и — вязкость 


жидкости (Па с); Р — пластовое давление (Па); г — расстояние до стенки скважины (м); # — время (с); с — 
теплоемкость нефти (Дж/кг К); р — плотность нефти (кг/м?); Т — пластовая температура (К); ХА, — 
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теплопроводность пористой среды (Вт/м3-К); и — скорость конвективного переноса тепла в пористой 
среде (м/с); = — коэффициент Джоуля-Томсона (КЛЛа); п — коэффициент  адиабатического 
расширения (К/Ла). 


Первое уравнение описывает процесс изменения давления в пласте, второе — процесс теплопереноса жидкости. 
Скорость конвективного переноса тепла в пористой среде определяется выражением 


КОР 
Е: (2) 
и д" 
Система уравнений (1) задается в пространственно-временном интервале: 
7, <Г< К, 8) 
0154, 
где г, — радиус скважины; А — расстояние до границы контура исследования (м); Ё/ — время исследования (С). 


Исходя из поставленных выше условий для скважины и рассматриваемой области для системы (1) 
начальные и граничные условия задаются следующим образом: 
1=0:РЕБ,Т=Т,г, <Г< К, 


ГЕГ, :РЕБ,Т=Г,,1>0, (4) 
г в и 0,#> 0. 
дг ди 


После определения начальных и граничных условий получена модель процесса фильтрации нефти в 
однопоровом коллекторе (1)—(4), характеризующая перераспределение давления и температуры в пласте. 

Первоначальным этапом решения системы является приведение его к дискретному виду. Для этого 
используется метод конечных разностей [8]. Строится равномерная пространственно-временная сетка: 
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М-1 
у (5) 
С, = |. = ит, и =0,1,...,М,т= т 
где й — шаг сетки по пространству; т — шаг сетки по времени; № — число узлов по пространству; М — 
число узлов по времени. 
После замены дифференциальных выражений разностными аналогами получается система уравнений: 
1 п-1 п-1 
г Рог, |Р" +" Р 
ТОРЫ ЧИ, 1 й о 
в" —Р" К1 О 5 о = 
Ф<с — 5 
т $ и р? 
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Система (6) представляет собой систему линейных алгебраических уравнений (СЛАУ), которая сводится к 
трехточечному уравнению общего вида: 


Аим Ву + Су" = ЕД=Ь,...М,п > 0. (7) 
Коэффициенты равны: 
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Р т й т _ 
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Система (7) решается методом скалярной прогонки [9, 10]. Задача является нестационарной, поэтому метод 
применяется для определения пространственного распределения давления и температуры пласта на каждом 
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временном слое [11, 12]. Задача решается последовательно, на первом этапе решается уравнение для давления, 
на втором — для температуры. 

Для автоматизации решения задачи разработано \уеб-приложение. В общем виде его структура представлена 
на рис. 1. В ней выделяются пять основных компонентов: входные информационные потоки, выходные 
информационные потоки, методы обработки информации, технические средства обработки информации, база 
данных [13, 14]. 


Методы обработки Технические 
информации средства обработки 
информации 


Входные Выходные 
информационные \\еБ-приложение информационные 
потоки потоки 


База данных 


Рис. 1. Структура приложения в общем виде 


Входными информационными потоками являются данные о параметрах модели, а выходными — графики, 
отражающие временные или пространственные зависимости забойного давления и температуры. Под методами 
обработки информации понимаются численные методы и методы интерпретации гидродинамических 
исследований. Техническими средствами обработки информации являются сервер, на котором размещено 
приложение, и персональный компьютер пользователя. Для хранения входных параметров и результатов 
расчетов используется база данных. 

Конечный вид программного обеспечения зависит от средств, выбранных для его разработки. Средствами 
разработки являются парадигмы и языки программирования, среды и технологии для разработки. Выбор 
инструментария влияет на способ реализации задачи, критерии качества программного обеспечения, 
представление данных и тип графического интерфейса пользователя. Для разработки приложения в качестве 
языка программирования для серверной части выбран РуФоп с использованием микрофреймворка ЕИазК, для 
клиентской части — Лауабсирё с Уче.]з [15, 16]. В качестве СУБД для проектирования и разработки базы 
данных используется РозетеЗ ОГ. [17]. 

Результаты исследования. Разработанное программное обеспечение представляет собой кросс- 
платформенное клиент-серверное \/еб-приложение. Хранение данных и вычисления происходят на \е$- 
сервере, с которым пользователь взаимодействует через браузер (поддерживаются все современные интернет- 
браузеры). Клиентская часть может использоваться на широком спектре платформ (ПК, ноутбуки, смартфоны) 
и операционных систем (\/т4о\з, МасО$, Апаго1а, 105). 

Большую часть главного окна \еб-приложения занимает средство построения графиков температуры и 
давления. Графики строятся после произведения расчета путем выбора соответствующего пункта меню. Также 
имеется возможность построения графиков по ранее произведенным расчетам, для этого выбирается 
необходимый расчет в левой части окна. Графики строятся по пространственным и временным координатам. 
Для изменения вида графика выбираются соответственные пункты в верхней части окна. 

Все произведенные расчеты хранятся в разработанной реляционной базе данных Но\са|с46, структура 
которой представлена на рис. 2. 
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Рис. 2. Структура базы данных Йо\’са]са6 


В общей сложности база данных содержит пять таблиц, которые находятся в третьей нормальной 
форме [18]. Все таблицы имеют первичные ключи. Главная таблица, са1сшайоп, хранит название расчета и дату 
его произведения. Две зависимые от нее таблицы, га4аз и Ите, хранят информацию о пространственно- 


временных координатах. Таблицы 


соответственно. 


фетрегаиге и ргезбиге хранят значения температуры 


и давления 


В разработанном приложении проведен ряд вычислительных экспериментов, позволяющих оценить 
результат его работы. Для проведения моделирования использованы значения из таблицы 1. 


Таблипа 1 


Начальные и граничные условия, параметры расчета 


Параметр Значение Единица измерения 
Радиус скважины, Г, 0,102 м 
Радиус исследования, В 100 м 
Время исследования, 1 24 Ч 
Начальное давление, Р. 25 МПа 
Конечное давление, Рк 30 МПа 
Проницаемость, К 1е-14 м? 
Пористость 0,2 
Вязкость нефти, | 1,10е-3 Па-с 
Сжимаемость нефти 1,58е-9 1ИПа 
Теплоемкость, с, 1,4 МДж/ мз-К 
Теплопроводность, А. 0,6 Вт/м-К 
Коэффициент Джоуля-Томсона, = 3,94е-7 К/ЛЛа 
Коэффициент адиабатического расширения, 7 3,15е-7 К/Па 
Начальная температура, Т, 363,15 К 
Количество точек по времени, М 1000 
Количество точек по пространству, М 1000 


С использованием данных из таблицы 1 построены графики распределения давления и температуры по 


пространству, представленные соответственно на рис. Зи4. 
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Рис. 3. Распределение давления по пространству в конечный момент времени 
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Рис. 4. Распределение температуры по пространству в конечный момент времени 


Как следует из данных на рис.3 и 4, после остановки скважины на гидродинамическое исследование 
наблюдается восстановление давления и температуры, причем чем дальше от стенок скважины, тем больше 
величина восстановления. На рис. 5 и 6 показано, как именно распределены давление и температура в 
различные моменты времени (6, 12, 24 ч.) после остановки скважины. Эта информация полезна специалистам 
по анализу и интерпретации гидродинамических исследований, поскольку позволяет оценить динамику 
давления и температуры и принять решение о продолжении или остановке исследования, а также проведении 
дополнительных мероприятий. 
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Рис. 5. Распределение давления по пространству в различные моменты времени 
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Рис. 6. Распределение температуры по пространству в различные моменты времени 


На основе произведенных расчетов установлено, что наибольшее влияние на поведение давления 
и температуры оказывает проницаемость пласта. На рис. 7 и 8 изображены графики динамики давления 
и температуры в зависимости от различных значений проницаемости вблизи от стенок скважины 
(низкопроницаемый коллектор — 1мД, среднепроницаемый — 10 мД, высокопроницаемый — 100 мД 
(1Д = 10-12м?)). Как видно на графиках, при увеличении проницаемости коллектора скорость восстановления 
давления и температуры увеличивается. На рис.7 изображен график для реальных промысловых 
данных (проницаемость коллектора — 80 мД). 
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Рис. 8. Динамика температуры для разных проницаемостей 
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На основе динамики давления можно определить время проведения исследования (рис. 7). Например, для 
высокопроницаемого коллектора не имеет смысла проводить исследование более чем -—6 часов, поскольку 
после этого времени давление практически не меняется, а простой скважины является неэффективным, 
поскольку идут потери добычи. 

Обсуждение и заключение. Разработанное \еБ-приложение может быть использовано для анализа 
и интерпретации данных гидродинамических исследований. \!еБ-приложение обеспечивает удобный 
пользовательский интерфейс и снижает затраты на поддержку. Оно может быть развернуто в существующей 
сетевой инфраструктуре или подключено к удаленному серверу для обеспечения полной функциональности. 
Применение приложения позволяет качественно и количественно оценивать поведение давления и температуры 
в однопоровых коллекторах. А хранение результатов расчетов в базе данных дает возможность проводить 
анализ с использованием результатов множества исследований, на основе которого можно составлять сводный 
отчет и формировать рекомендации по дальнейшему использованию скважины или месторождения в целом. 
Используемые модели являются устойчивыми и сходящимися. Графики строятся с помощью написанного \еф- 
приложения, и поэтому оно может использоваться на производстве. 

Приложение имеет широкие перспективы развития, например, можно реализовать другие модели течения 
жидкости или предусмотреть наличие нескольких фаз. Также на момент написания данной работы остро встал 
вопрос об отказе от зарубежных коммерческих продуктов в пользу отечественных симуляторов. Поэтому 
разработка подобных продуктов является перспективной. 
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Аннотация 

Введение. Обратные задачи представляют собой специфический тип задач, где изучаются последствия явлений 
с целью определения их причин. Они широко используются в научных исследованиях, особенно тех, что имеют 
дело с большими объемами экспериментальных данных. В представленном исследовании рассмотрены обратные 
задачи в машиностроении и диагностике конструкций. Эти области требуют точных методов для выявления в 
различных материалах внутренних дефектов, которые могут иметь критические значения для обеспечения 
безопасности и эффективности использования технических конструкций. Несмотря на множество имеющихся 
методов дефектоскопии существует потребность в инновационных разработках, способных обеспечить ее более 
высокую точность и эффективность. В данном исследовании объединены различные научные методы и 
технологии, оно открывает новые перспективы в неразрушающем контроле для обнаружения внутренних 
дефектов в различных материалах и структурах. Его цель — развитие и внедрение методов неразрушающего 
контроля на основе нейросетевого аппарата для повышения точности идентификации дефектов, а также 
разработка нейросетевой модели и оценка ее эффективности для усовершенствования процесса ультразвуковой 
визуализации внутренних дефектов в твердых материалах. В связи с этим задача, которую предстоит решить для 
достижения поставленной цели, заключается в создании надежного инструмента для точной визуализации 
размеров, форм, местоположения и ориентации внутренних дефектов в различных материалах. 

Материалы и методы. Применяется методика определения геометрических параметров дефектов в материалах 
с использованием неразрушающего контроля. Также используется метод, объединяющий моделирование 
распространения ультразвуковых волн в акустической среде и технологии искусственных нейронных сетей. Он 
выявляет нелинейные связи между геометрическими характеристиками дефектов и амплитудно-частотными и 
амплитудно-временными данными, полученными при анализе сигналов. Искусственные нейронные сети 
представляют собой модель, которая может обучаться на примерах, что позволяет эффективно решать задачи, 
которые сложно выразить в традиционных формах. В исследовании используется метод конечных разностей во 
временной области. Он применяется для идентификации и визуализации внутренних дефектов в материалах с 
использованием ультразвукового неразрушающего контроля и сверточных генеративных нейронных сетей. 
Результаты исследования. Разработана сверточная нейронная сеть для визуализации внутренних дефектов 
с использованием техник ультразвукового неразрушающего контроля. Эта нейронная сеть успешно определяет 
размер дефектов, их местоположение, форму и ориентацию с высокой точностью и надежностью. 

Обсуждение и заключение. Авторы подчеркивают ключевое влияние размера дефекта на точность 
ультразвуковой визуализации в различных сценариях. Проведенная валидация модели для трех различных 
случаев дефектов с разными механическими параметрами показала, что для успешной визуализации дефектов 
длина волны ультразвукового импульса должна быть в десятки раз меньше размера дефекта. При анализе влияния 
размера дефектов на точность работы нейронной сети выявлено, что ошибка визуализации увеличивается для 
дефектов меньшего размера. 
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Установлено также, что относительная скорость звука в материалах оказывает большее влияние на точность 
метода, чем относительная плотность материала. На основании полученных авторами результатов можно 
утверждать, что разработанные методики и технические решения имеют большое значение ДЛЯ будущих 
исследований в области дефектоскопии, обладают весомым потенциалом для научных и практических сфер 
применения. 


Ключевые слова: ультразвуковой неразрушающий контроль, дефекты, ультразвуковой отклик, сверточные 
нейронные сети 
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Введение. Обратные задачи в рамках научного исследования — это такие задачи, в которых осуществляется 
анализ последствий некоторых явлений или процессов с целью формулирования гипотез и выводов относительно 
их причинности. Обратные задачи встречаются в различных областях научных исследований, особенно в 
случаях, когда имеются большие объемы экспериментальных данных, позволяющие выявить характеристики и 
особенности процесса, который породил эти данные. Таким образом, обратные задачи предоставляют ученым 
возможность определить параметры процессов, недоступных при непосредственном наблюдении. 

В настоящее время существует установившаяся методология решения обратных задач, которая основывается 
на анализе реакции системы на определенные входные воздействия. Современные подходы к решению обратных 
задач подразделяют их на различные категории, что облегчает их формализацию и адаптацию в зависимости от 
конкретной ситуации. 

В машиностроении широко применяется методика решения обратных задач для диагностики механических 
систем. Изучение динамической реакции таких систем, включая анализ переходных характеристик, позволяет 
получить детальную и точную информацию о состоянии системы, что, в свою очередь, дает возможность 
определить ее производительность и качество. 

В области диагностики конструкций и изделий растет потребность в оценке текущего состояния технических 
систем с целью выявления и устранения потенциальных опасностей, которые могут привести 
к производственным убыткам или несчастным случаям. Такие оценки могут быть даны лишь 
высококвалифицированными специалистами в области реинжиниринга. 

С учетом этих аспектов была поставлена задача разработать и в дальнейшем оценить возможности 
применения нейросетевой модели для решения прикладных проблем ультразвуковой визуализации внутренних 
дефектов в твердых телах для различных материалов. 

Существует большое количество научных работ, посвященных данной теме. Так, в [1] авторы сфокусировали 
свое внимание на анализе композитных материалов, содержащих дефекты. Наличие дефектов в таких композитах 
может существенно влиять на их эксплуатационные характеристики и поэтому представляет особый интерес для 
исследователей. Для поиска дефектов каждый из исследуемых объектов подвергается ультразвуковому сканиро- 
ванию, в ходе которого фиксируются колебания и определяются их характеристики. Следующим этапом анализа 
является классификация полученных характеристик с использованием искусственных нейронных сетей. Этот ме- 
тод позволяет получить изображение области дефекта с высокой точностью. Такая визуализация отражает как 
характер, так и степень распространения дефекта в исследуемом материале. Важно отметить, что данная мето- 
дика обеспечивает четкую идентификацию зоны повреждения, а полученное изображение соответствует визу- 
альным наблюдениям. 

В работе [2] представлен метод акустической эмиссии в контексте оценки композитных материалов, пред- 
ставляющий собой перспективное направление исследований. Особое внимание уделяется идентификации аку- 
стических сигналов, связанных с микроструктурой волокон, в условиях высокого уровня шума и с использова- 
нием пленочных датчиков. Для исследований были подготовлены многослойные композиты на основе армиро- 
ванного стекловолокна, изготовленные различными методами. Эксперименты включали в себя нагружение ком- 
позитных образцов в статических условиях, в то время как РУОЕ-пленка была установлена на их поверхности. 
Акустические сигналы были аккуратно записаны и после этого были классифицированы при помощи искусствен- 
ной нейронной сети. Результаты данного исследования говорят о возможности при помощи нейронных сетей 
успешно классифицировать различные механизмы разрушения в композитных материалах. 

В исследовании [3] сделан анализ инцидентов обрушения, произошедших при раскопках в слоях горных по- 
род на гидроэлектростанции Вудунде в Китае. Его авторы представили усовершенствованную модель искус- 
ственной нейронной сети (ИНС) для прогнозирования вероятности и глубины таких обрушений. Однако следует 
сказать, что одиночный параметр не способен обеспечить абсолютной точности предсказания глубины обруше- 
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ний. В этом исследовании был предложен интегрированный метод интеллектуального прогнозирования, осно- 
ванный на применении искусственных нейронных сетей. Для его разработки проведены аналитическое и стати- 
стическое исследования, что позволило определить шесть ключевых входных параметров: глубина покрытия, 
отношение основных и малых напряжений, индекс геологической силы, метод выемки грунта, сила поддержки и 
ориентация породы. 

Модель ИНС была обучена с использованием данных, собранных в процессе раскопок на гидроэлектростан- 
ции Вудунде, среди них было 45 образцов для обучения и 6 тестовых образцов для верификации. Для оптимиза- 
ции параметров ИНС, таких как структурные характеристики и начальные веса, был применен генетический ал- 
горитм. После завершения этапа обучения обученная нейронная сеть была успешно применена для прогнозиро- 
вания глубины обрушений на нескольких участках раскопок. Результаты прогнозов демонстрировали высокое 
соответствие данным, полученным в ходе наблюдений на этих участках. 

В статье [4] представлен инновационный метод для детектирования дефектов в стержневых конструкциях. 
Он существенно отличается от других и состоит из двух этапов. Первый этап фокусируется на точной локализа- 
ции дефекта в стержне. Для этого используется анализ частот и форм колебаний, связанных с исследуемым объ- 
ектом. Результаты анализа форм колебаний играют решающую роль в решении задачи первого этапа. На втором 
этапе, после успешной локализации дефекта, создается аналитическая или конечно-элементная модель конструк- 
ции, учитывающая выявленные дефекты. С использованием этой модели проводится оценка степени поврежде- 
ния стержневой конструкции. Применение данного метода проиллюстрировано на примере вычисления парамет- 
ров дефекта в упругом консольном стержне. Этот расчет основан на данных экспериментальных наблюдений, 
полученных в конечно-элементном комплексе. Представленный алгоритм предоставляет новый инструмент для 
более точной и эффективной идентификации дефектов в стержневых конструкциях, что имеет большое значение 
для научных исследований и инженерной практики. Результаты моделирования, проведенного в рамках данного 
исследования, свидетельствуют о снижении ошибки прогнозирования осаждения. Это указывает на улучшение 
функциональной аппроксимации и, следовательно, на отличные обучающие характеристики этого метода, по 
сравнению с традиционными нейронными сетями с сигмоидными или другими функциями активации. 

В исследовании [5] предлагается инновационный метод диагностики дефектов в механической коробке пере- 
дач, основанный на анализе сигналов вибрации, полученных в ходе экспериментальных испытаний. Основное 
внимание уделяется раннему выявлению аномалий в работе зубчатых передач. Для эффективного анализа виб- 
рационных сигналов используется обработка данных, ориентированная на частотно-временную область. Одним 
из ключевых методов является применение непрерывного вейвлет-преобразования, что позволяет извлечь стати- 
стически значимые характеристики из сигналов. С целью оптимизации процесса обработки сигналов внедряется 
метод выделения особенностей, который способствует определению наиболее важных параметров, что, в свою 
очередь, снижает вычислительную нагрузку и ускоряет процесс анализа. Для выбора оптимальных входных при- 
знаков и дальнейшего снижения вычислительной сложности исследования применяются генетические алго- 
ритмы. Этот интеллектуальный подход способствует повышению эффективности процесса анализа данных. Сле- 
дующим важным этапом в рамках исследования является оценка эффективности предлагаемого метода. Для этой 
цели использованы передовые технологии нейронных сетей обратного распространения. Полученные результаты 
и производительность классификаторов тщательно оценивались на разнообразных вибрационных сигналах, по- 
лученных в ходе экспериментальных испытаний при различных условиях работы механических зубчатых пере- 
дач коробки передач. 

Исследования систем неразрушающего контроля и поиска дефектов в твёрдых телах зачастую затрагивают об- 
ласть моделирования соответствующих процессов [6]. Основное внимание направлено на решение обратных гео- 
метрических задач, тесно связанных с теорией упругости, для их дальнейшего применения при анализе плоской 
прямоугольной области. Задача заключается в восстановлении характеристик круговых полостей и трещин. 

Для эффективного решения обратных задач исследователи опираются на информацию, полученную о первых 
четырёх собственных резонансных частотах. Эти параметры представляют собой дополнительный набор данных, 
значительно обогащающий информацию о системе. Методика, применяемая для решения обратных задач, бази- 
руется на минимизации функционала невязки. Суть этой методики заключается в сведении к минимуму расхож- 
дения между измеренными входными данными и значениями, полученными в результате численного решения 
прямых задач. При этом предварительно задаются параметры дефектов в твёрдом теле. 

Таким образом, исследования в области неразрушающего контроля дефектов в твёрдых телах ориентированы 
на разработку и применение математических методов и моделирования, что позволяет точно определять и 
анализировать дефекты, скрытые внутри материала. Эти исследования нашли применение в области 
безопасности и контроля качества технических систем. Авторы исследуют возможности применения 
нейросетевых технологий в решении обратных задач механики, в частности задачи двумерной визуализации 
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внутренних дефектов. Данные методы получили широкое распространение в медицине, гражданской инженерии, 
неразрушающем контроле и других областях. 

В работе [7] предложена система для выявления трещин в металлических конструкциях и определения их глу- 
бины на основе анализа двухмерных изображений. Основной целью данной работы было создание доступной и 
удобной системы контроля, более экономически выгодной альтернативы дорогостоящим измерительным устрой- 
ствам. В рамках данного исследования была разработана специфическая стратегия обучения, а также рассмот- 
рено несколько архитектур нейронных сетей, способствующих более точному выявлению трещин и оценке их 
глубины. 

Процесс обучения нейронной сети включал в себя предоставление сети характеристик двухмерных профилей 
трещин, а также данных о максимальной глубине залегания трещины, измеренной с использованием лазерного 
микроскопа. Это позволило создать модель, которая способна точно определять глубину дефектов в материале. 
Следует отметить, что средняя ошибка в работе данной нейронной сети составляет всего около 18 % на тестовых 
образцах. Этот результат значительно превосходит итоги предыдущих работ, которые основывались на разнооб- 
разных стратегиях обучения и методиках анализа. Таким образом, данное исследование значительно повысило 
качество определения глубины трещин в металлических конструкциях. Его ценность заключается в разработке 
новых стратегии обучения и инструмента для более точной оценки глубины трещин. 

В исследовании [8] подчеркивается важность наличия размеченных данных, содержащих изображения, на ко- 
торых отмечены дефекты. Одним из новейших предложений в данной области является представленная модель 
глубокого трансферного обучения, предназначенная для извлечения характерных признаков внутренних дефек- 
тов на рентгеновских изображениях, полученных в результате инспекции авиационных композитных материа- 
лов. Эта модель представляет собой инструмент для автоматического анализа таких изображений, который поз- 
воляет оптимизировать процесс обнаружения дефектов и повысить эффективность неразрушающего контроля. 

В ходе исследования была произведена проверка эффективности предложенной модели на задаче обнаруже- 
ния дефектов-включений на рентгеновских изображениях. Анализ результатов экспериментов показал, что дан- 
ная модель способна достигать высокой точности (96 % по метрике Е1). Такие обнадеживающие результаты сви- 
детельствуют о потенциале данного подхода и его способности давать удовлетворительные результаты в задачах 
по идентификации дефектов. Это открывает новые горизонты в области неразрушающего контроля и позволяет 
более эффективно обеспечивать безопасность и надежность композитных материалов в авиационной индустрии. 

В работе [9] показан метод, который предлагает реконструкцию и визуализацию внутренних аномалий 
в форме трехмерных репрезентаций. Этот подход использует экономически и временно выгодную методику, из- 
вестную как импульсная термография. В данном контексте представлен новаторский метод, который способен 
быстро и точно оценить как глубину, так и толщину внутренних аномалий с помощью односторонней инспекции. 
Целесообразность и эффективность данного метода успешно демонстрировались в ходе исследования композит- 
ных и стальных образцов, содержащих полузакрытые воздушные зазоры. Результаты исследования подтвердили, 
что при предложенном методе ошибка в оценке общего объема трехмерных аномалий составляет всего лишь 
10 % для композитных образцов. Таким образом, данный инновационный метод предоставляет новые возможно- 
сти для неразрушающего контроля и визуализации внутренних дефектов в различных материалах и образцах, 
обеспечивая при этом высокую точность и экономическую эффективность. 

В работе [10] рассмотрены фундаментальные факторы, приводящие к отказам в работе твердотопливных ра- 
кетных двигателей. Один из ключевых аспектов, оказывающих критическое воздействие на целостность двига- 
телей, заключается в явлении отслаивания на границе раздела между ракетным топливом, изоляцией и защитным 
слоем. Современные методы диагностики и оценки структурной целостности ракетных двигателей, как правило, 
ограничиваются визуальным осмотром и анализом их компонентов. Данный подход имеет свои ограничения и 
требует совершенствования. В этом контексте представлен новаторский алгоритм, предназначенный для выяв- 
ления поверхностных дефектов обшивки, способных нарушить связь между твердым топливом и изоляцией ра- 
кетного двигателя. Для оптимизации работы нейронной сети проведен ряд тестов, в которых менялись различные 
параметры сети. Эти исследования позволяют достичь точности классификации в 99,08 % для класса исправных 
элементов, 90,66 % — для класса посторонних объектов и 99,48 % — для класса, связанного с дефектами. Полу- 
ченные результаты свидетельствуют о высокой эффективности предложенного метода и его потенциале для при- 
менения в области диагностики твердотопливных ракетных двигателей. 

Исследование, проведенное и опубликованное автором в статье [11], нацелено на изучение потенциала мето- 
дов глубокого обучения в контексте электромагнитной инверсии. Суть данного подхода заключается в примене- 
нии глубоких нейронных сетей, основанных на сверточной архитектуре, для анализа и обработки данных, полу- 
ченных в рамках данной задачи. Важным аспектом этого исследования является обучение с использованием об- 
ширных наборов синтетических данных, полученных при помощи трехмерного моделирования. Глубокие 
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нейронные сети здесь обучаются на множестве синтетических данных, представляющих собой информацию 
о распределении электромагнитных полей внутри объектов и сред. Эти данные предоставляют ценную информа- 
цию об их свойствах и составах, что делает метод электромагнитной инверсии более точным и эффективным. 
Эффективность предложенного метода была успешно продемонстрирована на различных моделях, которые 
имеют высокое практическое значение. К примеру, при мониторинге электромагнитного поля в области скопле- 
ния углекислого газа, находящегося под землей, инверсионный метод позволяет анализировать и контролировать 
данное поле с использованием источника на поверхности. 

Ранее авторы представленного сегодня исследования уже рассматривали задачи, сочетающие современные 
методы глубокого машинного обучения и хорошо зарекомендовавшие себя классические подходы к идентифи- 
кации дефектов [12—14]. В нынешней же работе акцент сделан на нейросетевую модель, которая эффективно 
выполняет двумерную акустическую визуализацию внутренних аномалий. Этот проект выведет на новый уро- 
вень область неразрушающего ультразвукового контроля, предоставив возможность провести системный анализ 
с использованием численных экспериментов. 

Цель представленной сегодня исследования — разработка и оценка эффективности нейросетевой модели для 
улучшения ультразвуковой визуализации внутренних дефектов в твердых материалах. Для выявления нелиней- 
ных связей между геометрическими характеристиками дефектов и сигналами ультразвукового контроля авторы 
предлагают использовать методику, основанную на решении обратной задачи механики твердого тела и приме- 
нении искусственных нейронных сетей. Основная задача при этом заключается в создании надежного инстру- 
мента для точной визуализации размеров, форм, местоположения и ориентации внутренних дефектов в различ- 
ных материалах. Численные эксперименты аппроксимируют процесс ультразвукового зондирования образцов, в 
которых присутствуют внутренние дефекты с различными геометрическими формами и механическими харак- 
теристиками. Эти эксперименты служат исходными данными для обучения нейронной сети и последующей про- 
верки ее на правильность функционирования и точность. Таким образом, можно утверждать, что данное иссле- 
дование охватывает широкий спектр научных методов и технологий, включая разработку нейронных сетей, мо- 
делирование ультразвуковых волн и анализ акустических данных. 

Материалы и методы. В данной работе применяются многослойные искусственные нейронные сети прямого 
распространения (Рее4-Юг\ууагА пеига| пебмогк$, ЕЕММ) и свёрточные ИНС (сопуо@опа| пепга| пебмогк, СММ). 
Для определения наиболее рациональной архитектуры сетей применяется аналитическое программное обеспече- 
ние. ИНС обучаются алгоритмами КРгор (устойчивое обратное распространение) и Адат (метод адаптивной 
оценки моментов). 

Алгоритм ЕВРгор основан на методе градиентного спуска. В отличие от стандартного алгоритма обратного 
распространения ошибки, КРтор для подстройки весовых коэффициентов использует только знаки частных про- 
изводных. Алгоритм использует обучение по эпохам, а коррекция весов происходит после обработки всего объ- 
ема обучающей выборки. Весовые значения обновляются следующим образом: 


: 720) 
а 0 
() ЕЙ дЕ” В 06, 
А 
а - 0, до, Е“ 
А ‚ до = +40, <0 ‚0<п <1=1, 
› ЕО дЕУ " 0%, 
—_ д , до, ВЕ $ 
у 0%, 
:. (+) — 0 (1) 
после чего происходит подстроика весов: ©; =, Е До, в 


Адат (метод адаптивной оценки моментов) использует сгруппированные средние как градиентов, так и вто- 
рых моментов градиентов. Если установлены параметры у’®, а функция потерь 9, где {отражает индекс текущей 
итерации, пересчёт параметра алгоритмом задаётся следующим образом: 


те) < В, т + +(1 _В, ут 


уе) < Ву + +(1-,) (ул), 
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где = является малой добавкой, используемой для предотвращения деления на 0, а Вл и В> являются 
коэффициентами забывания для градиентов и вторых моментов градиентов соответственно; | — базовая 
скорость обучения. Возведение в квадрат и квадратный корень вычисляются поэлементно. 

Изображения, обрабатываемые свёрточной ИНС, можно представить как тензор размера: 


Чит(7) = (#,6,4), 
где Й и и^ — размеры двумерного изображения; 4 — количество каналов изображения. Ядро свёртки К имеет 
следующие размеры: 
(К) = (4), 
где /— сторона квадратного фильтра свёртки. 


Операцию свёртки можно определить следующим образом. Объектом свёртки является двумерное изображе- 
ние Г. Ядро свертки К имеет размеры й Хх у’. Результатом свертки будет: 


о 


(1 "К — за Х 1-1. 


= = 
На этом операторе строится сверточный слой, который состоит из определенного набора ядер К › соответ- 


ствующих им смещений : 


по 4 


сопу (1, К), =@ + УК ЕЕ 


= Л К= 
где 4 — индекс канала изображения; Ь — составляющая смещения; с — активационная функция сверточного слоя. 

Обучаемыми параметрами данного слоя являются фильтр К и смещение р. 

В области машинного обучения существует несколько техник регуляризации. Для предотвращения переобу- 
чения сложных моделей ИНС с большим количеством параметров применяется метод регуляризации Тихо- 
нова (714е гезтезяоп, или Г2), в машинном обучении также называемый уе! 4есау. Коэффициентом регуля- 
ризации является ^, который управляет минимизацией нормы относительно минимизации потерь на обучающем 


наборе данных. Таким образом, для каждого веса у происходит прибавление к целевой функции Г, ( у, у) слага- 


^, хе ы 
емого: В |6|[= Ув. В рамках применяемого фреймворка Кегаз регуляризация применяется для каждого 


1=1 
свёрточного слоя ИНС-модели. 
Начальная инициализация весов проводится методом Хащег. Во время начальной установки весовых зна- 
чений применяется вероятностное распределение, как равномерное, так и нормальное с дисперсией 


Уаг(И/)=2/(п„ +"), где п, ип 


„„ — Количества нейронов в предыдущем и последующем слоях соответ- 
ственно. 
Дисперсия выходных значений линейного нейрона (без смещений) имеет вид: 


Пт 


У, Ух = Уу, (у, (х) = пу, (у, (У)у, (х ). 
1=1 


1=1 

Подразумевается, что веса и входные значения не коррелируют и имеют нулевое матожидание. Из этого 
получаем дисперсию вероятностного распределения, на основе которого проходит первичная инициализа- 
ция весов, и обеспечивающую сохранение дисперсии входных данных после прохождения через слой. 

В данной работе слои батч-нормализации (Горе и 5<езейу) активно применяются в большинстве разрабо- 
танных ИНС-моделей с использованием фреймворка Кегаз. Метод осуществляет нормализацию входных 
данных таким образом, что их матожидание равно нулю, а дисперсия — единице. При этом нормализация 
выполняется на каждом слое для каждого пакета данных: 


1 т 2 1 т 2 
Вх,....х„ьив = —У лов — —У (> и») . 
т т 
На основе этих характеристик преобразуется активационная функция: 


^ х- Ив ^ 
м = ———,5 > 0, у, = ух, В, 
ноев 
где В и у — параметры, которые возможно оптимизировать методом градиентного спуска на основе обучаю- 


щего набора. 
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Слои Огорои! применяются в ИНС-моделях в задаче идентификации дефектов на дорожном покрытии, кото- 
рые служат для предотвращения переобучения. При использовании данного средства регуляризации предотвра- 
щается взаимоадаптация нейронов на этапе обучения ИНС. Данный метод блокирует (отключает) нейроны во 


время обучения с вероятностью. Учитывая, что п(х)=х"/-+Ь — линейная проекция входного 4, — мерного 

вектора х на 4, — мерное пространство выходных значений и а(й) — активационная функция, применение 

Дгорош к данной проекции на этапе обучения можно представить как модифицированную функцию активации 
Л (п) =В0 а(№), 

где р=(х Вен а — 4л-мерный вектор случайных величин Х,, распределенных по закону Бернулли. Х,, имеет 


следующее распределение вероятностей: 

р, К =1 
К, — ) 

И О 


где К — все возможные выходные значения. 
Применение Огорои! к 1-му нейрону выглядит следующим образом: 


4, 
4. а У <, +Ь |,Х. = 
О =Ха Ух, +Ь = 2 ь о 
в 0, Хх, =0 


где Р(Х!=0)=р. 


Нейронные сети, которые описаны в данной работе, принадлежат к сетям прямого распространения. Входная 
информация составляется из данных численных экспериментов и, проходя через сеть, формирует на выходе 
значения, которые оцениваются метрикой, и на основе этих данных происходит корректировка внутренних 
значений весов нейронной сети. Применяется итерационный подход, который обеспечивает достижение 
удовлетворительного уровня точности работы нейросети. 

Методология идентификации дефектов. В данной работе рассматривается методика определения геометри- 
ческих параметров дефектов на основе данных, полученных в процессе воздействия на исследуемый объект. Эта 
задача входит в область неразрушающего контроля свойств и параметров объектов, ориентированных на обрат- 
ную задачу механики твердого тела. Для решения данной задачи применяется системный подход с ключевой 
ролью искусственных нейронных сетей. Суть данного подхода заключается в выявлении нелинейных зависимо- 
стей между различными геометрическими характеристиками дефектов, такими как их размеры, глубина распо- 
ложения, угол поворота и тип, и данными, полученными в процессе анализа. 

В качестве таких данных могут выступать амплитудно-частотные и амплитудно-временные характеристики 
сигналов, полученных в ответ на воздействие. Этот метод позволяет эффективно выявлять и оценивать пара- 
метры дефектов и обеспечивает точное и надежное контролирование состояния объектов. Искусственные 
нейронные сети — это вычислительная модель, используемая в машинном обучении, информатике и других ис- 
следовательских дисциплинах, которая основана на болыпом наборе связанных простых единиц, называемых 
искусственными нейронами, подобных аксонам в биологическом мозге. Связи между нейронами передают сиг- 
нал активации разной величины. Если объединённые входящие сигналы достаточно сильны, нейрон активиру- 
ется и сигнал перемещается к другим нейронам, связанным с ним. Такие системы можно обучить с помощью 
заранее приготовленных примеров, а не явно программировать. Этот метод дает значительное преимущество в 
тех областях, где решения и функции трудно выразить в традиционных формах. Нейронные сети, а также другие 
методы машинного обучения используются для решения различных задач, которые трудно решить с помощью 
традиционного программирования на основе таких правил, как компьютерное зрение и распознавание речи. ИНС 
обычно используются для моделирования сложных связей между входами и выходами с целью нахождения за- 
кономерностей в экспериментальных или иных данных (ава пе) [15, 16]. 

Разделение технической задачи на ряд подзадач, рассматриваемых отдельно друг от друга и решаемых с 
применением существующих средств и подходов, позволяет говорить о системном подходе к рассматриваемой 
задаче идентификации. Сочетание этих задач образует целую систему, которая решает производственные 
проблемы комплексно. 

Предлагается инновационный метод идентификации и визуализации внутренних дефектов с применением уль- 
тразвукового неразрушающего контроля (УЗНК) и сверточной генеративной нейросетевой архитектуры. Сигнал 
УЗНК, зарегистрированный на поверхности объекта контроля, пропускается через обученную нейронную сеть. В 
результате обработки сетью генерируется высококачественное изображение, представляющее собой тепловую 
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карту, демонстрирующую внутренний дефект объекта. Для иллюстрации эффективности предложенного метода и 
его перспектив в промышленных приложениях в качестве объекта исследования была выбрана металлическая 
стальная полоса. Внутри этой полосы предполагается наличие дефектов, которые могут имитировать как полости, 
так и жесткие включения. Также учитывается и возможность отсутствия дефектов. Следует отметить, что характе- 
ристики дефектов, такие как форма, размер и ориентация, могут значительно варьироваться, что представляет до- 
полнительную сложность для задачи визуализации. Методология данного подхода состоит в проведении серии чис- 
ленных экспериментов, на основе результатов которых обучается глубокая нейронная сеть. Для формирования обу- 
чающего набора данных выполняется систематическое изменение механических параметров, геометрических ха- 
рактеристик, а также положения внутренних дефектов. Особое внимание уделяется моделированию процесса рас- 
пространения акустической ультразвуковой волны внутри материала. 

Этот инновационный подход к идентификации и визуализации внутренних дефектов предоставляет возмож- 
ность более точного и эффективного контроля за качеством материалов и изделий в различных отраслях про- 
мышленности. Его применение может способствовать повышению надежности и безопасности конечных про- 
дуктов, что является важным аспектом в современной индустрии. На основе собранных данных строится и обу- 
чается оптимальная структура нейросетевой модели. В результате нейросетевая модель может реконструировать 
форму, размер и положение дефекта из необработанного ультразвукового сигнала, снятого на поверхности по- 
лосы. На рис. 1 показана блок-схема исследования. 


ЕОТЬ Ультразвуковые сигналы 


- я Обучение 
на верхней и нижней поверхностях 


численные эксперименты 


нейронной сети 
| ыыы | |- | | _ 
& + п > Пе | та й 
ражения внутренних де- |. 
фектов 


Рис. 1. Схема проведения серии численных экспериментов с использованием метода конечных разностей во 


временной области (ЕОТО). Фрагмент полосы в области дефекта и распространение ультразвуковой волны (слева). 
Построение набора данных, состоящего из ультразвуковых сигналов и соответствующих бинарных изображений дефекта 
(в центре). Построение архитектуры и обучение нейронной модели (справа) 


Метод конечных разностей во временной области. Он был предложен Кейном Йи [17]. Решение дифферен- 
циальных уравнений в контексте моделирования сред с дисперсными и нелинейными свойствами выявляет уни- 
кальный класс методов — сеточные методы. Эти методы нашли широкое применение в разнообразных областях, 
начиная с геофизических задач и заканчивая задачами в оптическом диапазоне. 

Один из эффективных методов исследования физических явлений в различных средах — метод конечных 
разностей во временной области. Он представляет собой ценный инструмент при решении акустических задач, а 
также может быть успешно применен в упругих средах. 

Основным уравнением данной акустической модели в плоской постановке является: 


ор _ п ду, , 0, ду, _ 1др 9», _ 1 др. 
Гей 0х 0] рах’& р ду 


Происходит дискретизация скорости и акустического давления частиц моделируемого объекта в простран- 


ственной сетке. Путем последовательных вычислений значений этих параметров можно анализировать распро- 
странение звукового поля в течение определенного временного интервала. Этот метод обеспечивает возможность 
изучения различных физических процессов и явлений, предоставляя ценные данные о поведении среды с разно- 
образными характеристиками. 

Метод конечных разностей во временной области хоть и был изначально разработан для решения акустиче- 
ских задач, обладает широким спектром применения и в других областях физики. В частности, он находит свое 
применение при моделировании распространения звука в различных средах, таких как воздух или жидкость. Та- 
кая универсальность метода делает его неотъемлемой частью научных и инженерных исследований в области 
физики и акустики. 
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Для проверки эффективности выбранной методики была разработана эталонная модель в расчетном па- 
кете СОМ$ОГ, предназначенная для анализа распространения акустической ультразвуковой волны. В рамках 
данного исследования были использованы два основных численных метода: метод конечных элементов (МКЭ) и 
метод конечных разностей (ЕОТО). Модель представляет собой структуру в форме квадрата, изготовленную из 
стали, с присутствием внутреннего отверстия. 

Значительные отличия в форме полученных сигналов объясняются различными методами задания источника 
ультразвуковых колебаний. В случае применения метода конечных элементов точки источника располагаются 
на окружности внутреннего отверстия, в то время как в методе конечных разностей во временной области источ- 
ник представлен точкой, находящейся в узле сетки. 

Однако следует отметить, что ограниченные размеры расчетной сетки не позволяют моделировать 
распространение акустических волн за пределами данной области. Для учета этого факта применяются 
специальные граничные условия, такие как поглощающие слои Мура или РМГ. (Реесё МайсВе4 Гауегз). Эти 
граничные условия значительно уменьшают отражение акустических волн от границы расчетной области, что 
позволяет симулировать прохождение ультразвуковой волны внутрь структуры с высокой степенью точности. 

Модель неразрушающего контроля. В качестве модели неразрушающего контроля используется стальная по- 
лоса, содержащая дефект. Еб размер — 100х15 мм. Поглощающие слои установлены на левой и правой границах 
полосы. Верхняя и нижняя части прутка связаны с областями, которым заданы механические параметры воздуха. 
Дефекты представлены в виде выпуклых многоугольников. Разнообразие внутренних дефектов задается путем 
изменения количества сторон многоугольника, их длины, расстояния между вершинами и центром дефекта. Фи- 
зические параметры дефектов изменяются относительно толщины полосы. Положение дефекта варьируется от 
0,3 до 0,7, длина сторон — от 0,1 до 0,5. Область ввода импульса датчика расположена на верхней поверхности 
полосы. Колебания считываются как на верхней, так и на нижней поверхности. Таким образом, моделируется 
эхо- и теневой методы неразрушающего ультразвукового контроля. 

В данном исследовании анализируется входной сигнал, представленный как фиксированное множество дис- 
кретных значений, зависящих от времени проведения эксперимента. Экспериментальное время было строго 
определено с учетом требования, чтобы зондирующий импульс, выпущенный на инициирующей стороне и 
направленный на противоположную поверхность тестируемого бруса, имел достаточное время для пройденного 
расстояния и возврата в точку источника. 

Следует отметить, что использованная частота зондирующего импульса составляет 10 МГц, что обеспечивает 
высокую разрешающую способность экспериментальной системы. Это позволяет более точно изучать свойства 
бруса и его структурные характеристики, основываясь на анализе отраженного сигнала. 

На рис. 2 показана схема проводимого численного эксперимента. Дефект расположен внутри стальной по- 
лосы. Источник ультразвукового сигнала обозначен красным цветом, точки считывания сигнала — зелеными 
метками. Сигнал также считывается в точках его излучения. 


Полоса Точки инициации Дефект 
и снятия сигнала 


Точки снятия сигнала 


Рис. 2. Схема численного эксперимента 


На рис. 3 представлены колебания ультразвуковых сигналов, полученных в указанных точках на верхней по- 
верхности полосы (справа) и на нижней поверхности полосы (слева). Это имитирует расположение датчиков на 
противоположных поверхностях объекта контроля, данная модель предполагает, что датчики могут переме- 
щаться вдоль полосы в поисках дефекта синхронно. 
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Рис. 3. Ультразвуковые сигналы на верхней и нижней поверхностях полосы, полученные в результате численного 
моделирования 


Модель нейронной сети. Сверточные нейронные сети (СММ) — это специальный нейросетевой инструмент 
для обработки дискретных данных (двумерные изображения, одномерные сигналы) [18]. В настоящее время это 
один из наиболее быстро развивающихся и перспективных инструментов глубокого обучения [19-21]. 

Сверточные нейронные сети эффективно применяются в различных областях, включая анализ видео и обра- 
ботку временных рядов. Они позволяют обрабатывать временные ряды как одномерные массивы данных и вы- 
делять важные признаки в них с помощью операции свертки. 

Для визуализации дефектов в объектах контроля используется архитектура модели, включающая слои одно- 
мерной свертки (Сопу1О) и слои понижающей дискретизации (МахРоо|пз), для извлечения признаков из уль- 
тразвуковых сигналов. Затем данные передаются в полносвязный слой для анализа в контексте задачи визуали- 
зации дефектов. Вторая часть сети генерирует изображения, отображающие характеристики дефектов. 

Для создания окончательного изображения дефекта данные из полносвязного слоя обрабатываются и пода- 
ются на двухмерный слой. Затем с использованием обучаемых сверточных слоев (Сопу2) Тгапзрозе) формиру- 
ется финальное изображение. 

На рис. 4 показана архитектура сверточной нейронной сети с указанием размеров входных данных и количе- 
ства сверточных фильтров для каждого слоя. 


| Изображение 


1512 504х32 168х64 56х128 18х256 6х512 3072 1024 32х32х|1 64х64х8 128х128х1 дефекта 


Входной сигнал 


== Входной слой 77% Одномерный НЕ Изменение размера Е? Сопу2) Тгапзрозе + КеГ О 
= Сопу1 + ВеГ О + МахРооНп Г] епзе + ОгороцЕ ет Сопу2) Тгапзрозе + Ве О ЕЕ Сопу2) Тгапзрозе + Это 


Рис. 4. Структура модели сверточной нейронной сети 


Обучение и валидация модели. В рамках комплекса численных экспериментов были рассмотрены три случая. 
В случае А дефект моделирует полость, а параметры материала дефекта соответствуют параметрам воздуха. В 
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случае В дефект моделирует жесткое включение и параметры дефекта имеют параметры материала, который яв- 
ляется более жестким, чем сталь. В данном случае это карбид вольфрама. В случае С варьируются механические 
параметры дефектов, которые соответствуют параметрам газов, жидкостей и твердых материалов. 

Для случаев А и В было решено в общей сложности 50 000 задач для различных геометрических параметров 
дефекта. Из них 35 000 экземпляров используются для обучения, 10 000 — для тестирования и 5 000 — для ва- 
лидации. В случае С решается 60 000 задач с одинаковым процентным соотношением тестовых и валидационных 
образцов. 

В процессе оценки производительности модели ключевыми метриками являются минимальная ошибка обу- 
чения и сокращение разрыва между ошибкой обучения и тестирования. Эффективность модели машинного обу- 
чения подвергается анализу с использованием различных методов регуляризации, таких как регуляризация Ти- 
хонова, также известная как Г.2-регуляризация. Инициализация весов, проводимая, например, методом Хаутег, 
представляет собой значимый этап обучения нейронных сетей. 

Одним из средств, применяемых для уменьшения внутреннего ковариационного сдвига и стабилизации обуче- 
ния нейронных сетей, является метод пакетной нормализации. С целью повышения стабильности и обобщающей 
способности модели применяются разнообразные приемы регуляризации, такие как ранняя остановка и Огороцй. 

Необходимо отметить, что в данной исследовательской работе не применяется аугментация данных ввиду 
обширного объема обучающей выборки. Также не проводится балансировка набора данных, что, тем не менее, 
имеет свои преимущества в связи с нормальным распределением дефектов по размерам. 

Следует подчеркнуть, что обучение нейронной сети представляет собой задачу оптимизации, где выбор соот- 
ветствующей функции потерь играет важнейшую роль. Например, в задачах компьютерного зрения и семанти- 
ческой сегментации коэффициент схожести Жаккара может оказаться весьма значимым. В данной работе рас- 
сматривается применение коэффициента схожести Жаккара в контексте задачи обнаружения и сегментации объ- 
ектов, что позволяет достичь более точных и устойчивых результатов при анализе изображений: 

_ [АоВ| _ [Ао8| 


|Аов| |А+|В|-|АСВ| 


юи 


В данном контексте функция потерь выражается как обратное значение коэффициента сходства Жаккара (1— 
То) и представляет собой метрику, которая оценивает различие между двумя экземплярами. Метрика Жаккара 
известна в области анализа данных и машинного обучения как индикатор пересечения и сходства множеств. 

Для оценки эффективности модели нейронной сети в процессе валидации данная метрика применяется как 
один из критериев. Она позволяет определить степень соответствия между выходом модели и эталонными дан- 
ными, что важно для анализа и улучшения работы нейронных сетей в задачах обработки информации: 

т 2х6 
50 +1 

В настоящем исследовании использован один из наиболее распространенных алгоритмов оптимизации, име- 
нуемый Адат. Для решения поставленных задач машинного обучения и обучения нейронных сетей применяются 
одни их наиболее признанных и широко используемых библиотек — ТепзогЕ1о\\ и Кегаз. Эти программные про- 
дукты выступают в настоящее время как стандартные инструменты для создания и обучения нейронных моделей 
и включают в себя передовые технологии в данной области. 

Для достижения высокой точности работы моделей было проведено обучение в течение 30 эпох в каждом из 
трех рассмотренных случаев: А, Ви С. 

Нейросетевой подход для решения обратных задач уже зарекомендовал себя с хорошей стороны [22—24]. С 
углублением в область машинного обучения и последующим развитием новых методологий открываются пере- 
довые инструменты для анализа данных и решения традиционных задач в механике и дефектоскопии. Одним из 
значительных научных достижений в данной сфере представляется сверточная нейронная сеть, созданная специ- 
ально для решения проблемы визуализации внутренних дефектов с использованием методов ультразвукового 
неразрушающего контроля. Примечательной характеристикой данной модели является её относительно неболь- 
шой объем, составляющий 37,2 миллиона параметров. Этот факт обеспечивает высокую эффективность процесса 
обучения и расширяет область её применения в разнообразных смежных областях. Для более эффективного мо- 
делирования распространения ультразвуковых волн в данном исследовании был выбран метод конечных разно- 
стей во временной области. Этот выбор метода существенно увеличил скорость вычислений по сравнению с 
предыдущими подходами, применяемыми в ранее проведенных научных исследованиях. 

Результаты исследования. В данной работе рассматривались три состояния внутреннего дефекта. В случае А 
дефект моделирует полость, при этом параметры среды в области дефекта соответствуют параметрам воздуха. В 
случае В дефект моделирует жесткое включение, при этом параметры среды имеют параметры более жесткого 
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материала, чем материал полосы. В данном случае это карбид вольфрама. В случае С проводится вариация меха- 
нических параметров дефекта, которые соответствуют параметрам газов, жидкостей и твердых материалов. 

Результаты этого исследования, представленные на рис. 5, демонстрируют эффективность нейронных сетей 
в анализе и восстановлении информации о дефектах. Визуализированные фрагменты полосы в области дефекта 
отображаются на графике. Реальный дефект обозначен сплошной черной линией, тогда как форма и положение 
дефекта, предсказанные нейронной сетью, представлены пунктирной красной линией. Анализ рисунка позволяет 
утверждать, что нейросетевая модель способна восстановить размер дефекта, точно определить его местополо- 
жение, частично воспроизвести его форму и даже установить ориентацию. Эти результаты подтверждают потен- 
циал нейронных сетей в решении задач обработки и анализа изображений с применением высокой точности и 
надежности. 
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Рис. 5. Результаты визуализации дефектов. Реальный дефект выделен черным цветом. 
Прогноз нейронной сети выделен красным цветом. 


Валидация проводилась для каждого из случаев. Для оценки общей эффективности модели использована мет- 
рика Е1. В целом точность предложенного метода находится на высоком уровне. Среднее значение Е] для всего 
валидационного набора составляет 92 % для случая А, 90 % — для случая В и 85 % — для случая С. 

На рис. 6 показана зависимость точности визуализации дефекта от его размера. Размер дефекта — это макси- 
мальное расстояние между вершинами многоугольника, имитирующего дефект. Полученные данные показы- 
вают, что для успешной визуализации дефектов длина волны ультразвукового импульса должна быть в десятки 
раз меньше размера дефекта. На рисунке также показано влияние несбалансированного набора данных на точ- 
ность работы нейронной сети. Для дефектов наибольшего размера ошибка визуализации увеличивается. 
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Влияние размера дефекта на точность визуализации 
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Рис. 6. Зависимость точности визуализации дефектов от их размера для случаев А (слева) и В (справа). Черной штриховой 
линией показана длина волны зондирующего ультразвукового импульса 


Авторы также проанализировали влияние механических параметров дефекта на точность их визуализации. 
Для каждой нейронной сети из случаев А, В и С проведена валидация на наборе данных с различными механиче- 
скими параметрами дефектов. Данное исследование показывает, насколько хорошо обученные нейросетевые мо- 
дели способны проводить визуализацию различных типов дефектов, в том числе дефектов, которые отсутство- 
вали в обучающих наборах. Для простоты восприятия параметры плотности материала и скорости звука в мате- 
риалах представлены как относительные: 
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На рис. 7 показано, как точно нейросетевые модели А, Вис проводят визуализацию дефектов для широкого 
спектра параметров. В случае, когда нейронная сеть обучена только на наборе данных А, когда дефект моделирует 


полость, приемлемая точность достигается только для дефектов с С > 0,7. В случае В, когда дефект модели- 


тейайуе 
рует жесткое включение, точность работы нейронной сети существенно ухудшается при отклонении значения 
С 


Ге 


ре ОТ заданного в эксперименте. В случае С, когда нейронная сеть обучалась на наборе данных с различными 


параметрами дефектов, результаты работы устойчивы на всей области изменения параметров. Исключением явля- 
ются дефекты со значением С 


тейайуе 


—0. В этих случаях скорости звука в материале полосы и материале дефекта 


очень близки друг к другу. Также можно сделать вывод, что на точность предложенного метода в большей степени 
влияет изменение параметра С, чем р. Таким образом продемонстрированы границы применения различных 


нейросетевых моделей ДЛЯ идентификации дефектов с различными механическими параметрами. 
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Рис. 7. Влияние параметров дефекта на точность работы нейронных сетей. 
Нейронная сеть А слева, нейронная сеть В в центре, нейронная сеть С справа 
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Таким образом, можно заключить, что разработанная авторами модель нейронной сети представляет собой 
важный научный шаг в области неразрушающего контроля, обеспечивая эффективные инструменты для визуа- 
лизации внутренних дефектов и создавая новые перспективы для применения в индустрии и научных исследова- 
ниях [24—25]. 

Обсуждение и заключение. Полученные результаты предоставляют новый взгляд на область неразрушаю- 
щего контроля и его потенциал в обнаружении внутренних дефектов в различных материалах и структурах. 
В представленном исследовании предложен новаторский метод, основанный на применении сверточных нейрон- 
ных сетей, с целью экспрессной диагностики дефектов в структуре материалов. Ключевым аспектом данной 
научной работы является продемонстрированная потенциальная эффективность применения сверточных нейрон- 
ных сетей для повышения качества ультразвуковой визуализации. Авторы разработали модель, предназначенную 
для неразрушающего контроля, а также сформировали соответствующие наборы данных, используемых для обу- 
чения искусственных нейронных сетей. Эта сверточная нейросетевая модель способна не только точно опреде- 
лять форму, но и положение, и ориентацию внутренних дефектов в твердом материале. Полученные эксперимен- 
тальные данные демонстрируют информативность ультразвукового сигнала, а также его способность точно пе- 
редавать реальные характеристики внутренних дефектов. 

Авторы проанализировали несколько типов дефектов, в каждом из которых уделяется особое внимание: в пер- 
вом случае дефект классифицируется как А и моделируется в виде полости, во втором случае — как В и пред- 
ставляет собой жесткое включение, в третьем — как С, с различными механическими параметрами дефекта. Во 
всех сценариях исследования выявлено важное воздействие размера дефекта на точность ультразвуковой визуа- 
лизации. Кроме того, полученные результаты однозначно подтверждают, что выбор частоты зондирующего сиг- 
нала имеет критическое значение и должен быть оптимизирован так, чтобы длина волны была значительно 
меньше размеров внутренних дефектов. 

После проведения системной валидации сверточных нейронных сетей в контексте каждого из вышеупомяну- 
тых типов дефектов (А, В и С) для различных значений механических характеристик дефектов был установлен 
диапазон значений механических параметров дефекта, в пределах которого предложенный метод не обеспечи- 
вает достаточно качественную визуализацию. 

В заключение следует отметить, что разработанные в ходе данного исследования методики и технические 
решения имеют большое значение для дальнейших исследований в области дефектоскопии и решения обратных 
задач в широком контексте научных и практических приложений. 
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